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FOR THE INCREASE OF COMPLEX MECHANIZATION AND 
AUTOMATION OF PRODUCTION 


During May 12-16 of this year there was held an All-Union Conference of industrial workers on the question 
of developing complex mechanization and automation of productive processes, and of increasing the output of 
instruments and means of automation, At this conference, leaders of the party and of the government, and in- 
dustrial and scientific workers, spoke of the profound successes in mechanization and automation in our country, 
the blooming of creative force in Soviet science and technology, a striking demonstration of which was given 
by the launching of the third Soviet sputnik, outfitted with the most modern contrivances of measurement and 
automation technology. 


At this stage of its development, the economy of the Soviet Union has witnessed the maturing of a tech- 
nological revolution based on automata, atomic energy and radio-electronics, We are entering upon an epoch 
of creating qualitatively new forms of production which will allow productive work to be increased several 
times over. In these conditions, complex mechanization and automation of production in the national economy 
of the USSR is an extraordinarily important economic, political and technological problem, 


However, in this most important area there are still grave inadequacies, Despite the resolutions of the 
twentieth congress of the Communist Party (of the Soviet Union) on the broad automation of a number of branches 
of industry, this work is proceeding slowly, Theoretical investigations and preparations of technological processes 
for mechanization and automation, and also the modernization of technological apparatus in view of the require- 
ments of automation,are being performed on a completely inadequate scale, In the production lines of 
the machine industry there is an inadequate “specific gravity" of semi-automatic and automatic machine stations 
and automatic continuous-flow lines, many technological processes still requiring the expenditure of physical 
labor, There is a low level of mechanization of the warehousing and transportation processes, A unifying techno- 
logical policy is lacking for the cardinal questions of automation and the production of the technological means 
for mechanization and automation, There is flux and incoherence in tool manufacture, lagging in the areas of 
creating new methods and apparatus for automatic control and regulation, in particular, transducers for controlling 
physicochemical processes, electronic controlling machines and other automatic and remote control apparatus, 
and also a series of machines, mechanisms and other means for complex mechanization, 


The theoretical questions constituting the scientific stumbling blocks in automation have been inadequately 
worked out ;there has been insufficient training of the personnel for the tooling and automation of various branches 
of industry, There is lagging in economic science when it comes to estimating the economic effectiveness of 
automation; the attempts at automation have been inadequately generalized, and the state of technical informa- 

“tion in these areas is quite poor, 


Another serious inadequacy is the slow tempo with which the latest advances are incorporated in the areas 
of mechanization and automation of production, There is a large hiatus between the time of a scientific dis- 
covery and the moment of its realization, 


‘The most important problems in the domain of automata are the following: 


1, The creation of automatically working models of production complexes in various branches of the 
national economy on the basis of new technology and attainments in the techniques of control and regulation, 
plus the rapid dissemination of the results obtained from such models, 


2, The creation of one state-wide system of production of the technical means for automata and remote 
control, and such a development of industry producing these means which will guarantee the continuous satis- 
faction of the growing demand for them, 
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3. The development of research and development organizations and the effective allocation of them in 
the economic zones in accordance with their features, 


Before us the grandiose promises of automation open up: the creation of automated mines, automatic 
factories, tele- automated productive complexes, With each year automation will acquire ever greater value 
in every domain of our lives, beginning with the organization of our way of life and ending with space travel, 


The problem of our science consists in paving the way toward new solutions to these problems, and in 
providing ourselves with the means to realize these solutions. 


It is not to be doubted that the scientists working in the domains of automata and remote control will do 
all possible in carrying out the distinctive role which science plays in the matter of complex mechanization 
and automation, and that they will do their part in the building of Communist society. 
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THE DETERMINATION OF AN OPTIMAL SYSTEM BY SOME 
ARBITRARY CRITERION® 


V. S. Pugachev 


(Moscow) 


A method is expounded for the determination of an optimal system, by some arbitrary 
criterion of the Bayes type, out of the class of all functions to which this criterion can 
meaningfully be applied, very general assumptions being made as to signals and noise, The 
solution of the problem of determining the operator of the optimal system reduces to the 
finding of certain linear operators, and the minimum of some function or functional, 


1. Posing of the Problem 





The general problem of the optimal systems consists of this, that being given a random function Z(t), 
t€T, at the input, there must be developed at the output a certain random function W*(s), which would be as 
close as possible, in the sense of the criterion chosen, to the random function W(s), s€S, The random function 
Z is generally the sum of the effective signal plus noise, and the random function W is the result of some trans- 
formation of the effective signal,** For brevity, therefore, we shall henceforth call the random function W the 
effective signal, and the random function Z the observed random function, The random function W* , which 


the system must present at its output as the result of reproducing the random function W, we shall call the esti- 
mate of signal W, 


In various practical problems, different criteria of optimality are chosen, depending on the purpose 
of the proposed system, All criteria employed in practice can be presented as a minimum condition on the 
mathematical expectation of some function or functional of the effective signal and its estimate, Thus, the 
most general criterion possible for an optimum has the form 


p= M[r(W(s), W’(s))] = min, (1) 


where r_ is some function of the signal and its estimate, 


This function may depend only on the current value of the signal and its estimate, i.e,, be a function in 
the ordinary sense, or it may depend on all values of the effective signal and its estimate in some region S, i.e., 
be a functional, If the functions is given in advance, and its form does not depend on the manner in which 
the estimate W* of the effective signal is determined, then criterion (1) is a Bayes criterion, Of all the 
criteria known in practice, only the criterion of minimum loss of information does not fall into the class of 
Bayes criteria. Theoretically, the criterion of N, 1, Andreev [3] is not a Bayes criterion in all cases, However, 
in all cases in practice when it appears expeditious to use Andreev's criterion, this criterion does fall into the 





* The results of this work were presented on December 24, 1957, and January 7 and 21 and February 4, 1958, 
at the seminar on probabilistic methods of the theory of automatic control at the Institute of Automation and 
Remote Control of the AN SSSR, 

** We will consider here only real random functions and functionals, 
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class of Bayes criteria, It is thus possible to consider that all optimal criteria arising in practice, with the ex- 
ception of the criterion of minimum loss of information, are Bayes criteria, 


In certain cases, the function r contains undetermined parameters which depend on the manner of de- 
termining the estimate W* of the effective signal, it being then necessary so to choose these parameters that 
certain supplementary conditions are satisfied, Criterion (1), if it contains such a function r, might be called a 
conditional Bayes criterion, Among the conditional Bayes criteria is, for example, the well-known Neyman- 
Pearson criterion for revealing a signal in the presence of noise, The theory of N, I, Andreev shows that the 
criteria of the class investigated by Andreev are, in the general case, also conditional Bayes criteria, 


The correspondence between the observed random function Z and the estimated effective signal W* is 
established by the system's operator A: 


W’ (s) = AZ (t). (2) 


Consequently, the problem of determining the optimal system in accordance with an arbitrary criterion of 
type (1) reduces to finding an operator A guaranteeing the minimum mathematical expectation of the function 
r. 


In many cases, the effective signal presents itself as the sum of a linear combination of several known 
functions with linear coefficients plus a certain random function, Therefore, the observed random function Z 
can be given in the form 


N 
Z(t)= ¥ Unen()+ X (0), (3) 


h=] 


where @, are the given functions, Uj; are random variables and X is a random function representing the sum of 
the nonregular portions of the signal and the noise, 


The random function W, which is to be reproduced, frequently is the result of some linear transformation 
of the signal, In such a case it may be represented by the formula 


N 
W (s) = Si Unda (s) + ¥(s), (4) 


h=1 


where the functions ph are the result of subjecting the given functions g, to the linear transformation, and the 
random function Y is the result of applying the same linear transformation to the nonregular portion of the signal, 


We give here a method for determining the optimal operator A in accordance with an arbitrary Bayes (or 
conditional Bayes) criterion of type (1) for the case when the observed random function Z, and the signal to be 
reproduced W, are represented by Formulas (3) and (4), assuming that the random function vector [X(t), Y (s)] 
is distributed normally and is statistically independent of the random vector (U,,..., Uy). Without loss of 
generality, it may be considered here that the mathematical expectation of the random functions X and Y is 
identically zero (1, 2). 


The arguments t ands of the random functions will be assumed to be arbitrary scalar to vector variables 
which, in particular, can have certain discrete components, Therefore, all the theory presented here is applicable 
both to scalar and to vector random functions, To apply the theory to vector functions, it suffices to consider the 
components of each vector function as a scalar function of its argument and ordinal number, Thus, the method 
presented here for determining optimal systems is applicable to both one-dimensional and multidimensional 
systems, 


As is shown by the formula 
M [r(W, AZ)| = M [M[r (W, AZ) | Z)], (5) 
it suffices for the solution of the posed problem to find an operator A which guarantees a minimum conditional 
mathematical expectation of the function r with respect to the observed random function Z for each possible 


realization of the random function Z (with the possible exception of certain realizations having zero total pro- 
bability of occurring): 
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le 


M (r(W, AZ) | Z] = min. (6) 


Our problem, in its general form, is solved by the method of the canonical decomposition of random functions, 
which gives a sufficiently simple algorithm for finding the optimal operator, 


2. Solution of the Problem For the Criterion of Minimum Mathematical 





Expectation of a Function of the Current Values of the Effective Signal 





and its Estimate 








We present the vector random function [X (t), Y (s)] by some canonical decomposition (1): 
X (t) = z Viz, (lt), Y(s)= dS Voy (s), (7) 


where the V,, are noncorrelated random variables with zero mathematical expectation, As is well known, it is 
always possible to do this and, moreover, to do it in a nondenumerably infinite number ef ways(1Jj, After the 
decomposition of (7) has been constructed, we define a system of linear functionals, 9'” ) , satisfying, together 
with the function coordinates x,,, the condition of biorthogonality: 


os, i By q (8) 


This also is always possible to do, using, for example, the methods presented in [1] or in[2), Having done 
this we will also have satisfied the condition 


1 pag 
ty (t) = py, mK: (t, z), (9) 


where D,, are the dispersions of the random variables V,,, and Ky (t, u) is the correlation function of the random 
function X, The random coefficients of decomposition (7) are expressed by the formula 


V. = Qmyx (t). (10) 


We now determine the random variables Z,, = a'Y)z which, based on (3) and (10), are expressed by the 
formula 


N 
Ze = QZ (t) = D, > ann, + Vivo 


(11) 
h=1 
where 
i v 
Gyn = — OM op, (h=1,...,.N). (12) 
Then, taking (7) into account, we obtain 
N 
2.4% (1) = J Un SY) Den, (t) + X (2). (13) 
v h=1 v 


A comparison of Formulas (3) and (13) shows that in the case when all the functions gy, can be represented 
by the decomposition by the coordinate functions 


9 (t) = x Dyer, (t) h=1,...,N), (14) 
¥y 
then the following equality holds; 


> Z x, (t) = Z(t). (15) 
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Thus, in the case when all the functions g, can be represented by decomposition (14), the observed ran- 
dom function Z is completely defined by the given random variables Z,,, and, consequently, 


M [r(W, AZ) | Z| = M [r(W, AZ) | Z,). (16) 


We assume initially that in the canonical decomposition in (7) of the random function Y there are no 
random variables V, which do not also enter into the decomposition of the random function X, Then, re- 
placing the random variables V,, in the second formula of (7) by their expression from (11) and substituting 
the result in (4), we get 


N 
W (s) = a Une, (8) + dy (s), (17) 
h=1 v 


where 
on (8) = Yn (8) — Si Dytonys (8) = ts +N). (18) 


Thus, giving the joint distribution function of the random variables Z,, and Up, is equivalent to giving the 
joint distribution function of the random functions Z and W, In the case when the random function Y has, in 
decomposition (7), a finite number of random variables V,, which do not enter into the decomposition of the 
random function X, then in Formula (17) there appears a linear combination of such random variables V,,. 
Obviously, this linear combination may be joined with a linear combination of the random variables Up. Con- 
sequently, the case where the decomposition of the random function Y contains a finite number of random 
variables V,, which do not enter into the decomposition of the random function X leads to the previous result 
with only a simple increase in the number of terms entering into the sum of Formula (4), 


We note that a sufficient condition for convergence of the series in (14) and (18) is the convergence of all 
the series [1, 2): 


>) D2, (h=1,...,N). (19) 


On the basis of the results presented, the application to the observed random function Z of an arbitrary 
operator presents itself as a function of the independent variable s and the random variables Z,,, and, con- 
sequently, the optimal operator has the form 


AZ = 2 (s, Zy), (20) 


where p is some as yet unknown function, With this, the conditional mathematical expectation of the function 
r_ with respect to the observed random function is expressed by the formula * 


M[r (W, AZ) Z,) _ 


) ) N 
om { ee \ r (> Upp (S) + AE (s), u(s, z,)) 1 (un | Z,)du,...duy, (21) 


where f(up |z,,) is the conditional probability density of the random vector (Uy, . . . , Uy) with respect to the 
random variable Z,,, 


Based on well-known formulas of the theory of probability [1], the conditional probability density of the 
random variable U;, with respect to the random variable Z,, is expressed by the formula 


f (up) fr (2, | 4p) 
f (ur | 2) = >-— uy) fr (zy | Up : 
\ abe \ f (u,) fre, | up, du... duty 
* We denote random variables by capital letters, and the corresponding variables of integration and probability 
density arguments by lower-case letters. 





(22) 
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where f(u,) is the unconditional probability density of the random variable U;, and f(z,,| up) is the conditional 
probability density of the random variable Z,, with respect to the variable U;. In order to find the conditional 
probability density f4(z,| up), we note that the random variables V,, are distributed normally, since they are 
the results of a linear transformation (10) of a normally distributed random function X, Consequently, they are 
not only noncorrelated, but also independent, and their joint probability density equals the product of their prob- 
ability densities, It follows from this and from Formula (11) that the conditional probability function of the 
random variable Z,, with respect to the variable Up, will be normal, and the joint conditional probability density 
of all the random variables Z,, with respect to the variable U;, will equal the product of the corresponding 
conditional probability densities of the individual random variables Z,,, Taking into account that the random 
variable V,, has zero mathematical expectation and dispersion D,,, we obtain, on the basis of (11), for the con- 
ditional probability density of the random variables Z,,, with respect to the variable Up, the formula 


s Be (»- Dy M4 contin) fies 


Te, (Z | Un) = 





4 
V 2nD, 


For the functions f, in Formula (22) we substitute the products of Expression (23) corresponding to all values 
of the index v, obtaining the following expression for the conditional probability density of the random variables 
U}, with respect to the random variables Z,,: 


N N 
4 
f (un | 2) = p (2y) f (Un) exp bs Up SY) %vnZy or , UpUg > Dy} ; (24) 
A= v P.q=1 v 


where p(z,,) is some function of the variable z,,, the concrete representation of which is immaterial for our 
problem, 


Substituting Expression (24) in Formula (21), we are led to the conclusion that our problem reduces to de- 
termining the function pu for which the integral 


co ) N , 
t in \ ‘al \ r( un (8) + Dd) zy, (8), 1 (6%) f (an) x 





—co poe h=] . 
N , N 
x exp > Uh >) avnzy — > ape du, ... duy, (25) 
helo P.q~1 
where 
Ang = >) Dytvpi%vq (p,q =1,.-.,) (28) 
has the least possible value for all possible values of the variabless and z,, 
We consider the integral 
r) © N 
I (8, %. May +++ + WX) = \ si \ r( anc (6) + tor) (tn) % 
—o —o h=1 
N - N 
x exp | 3) um —-b > apap} diy. day. (27) 
awl P.q=1 


For fixed values of 8, No, 4» ++ » 2» the integral in (27) is a function of the variable x, Let X» be the 
value of the variable x for which the integral in (27) has the minimum value, It is obvious that x» depends on the 
values of the variables s, 19, Ny, 2+» 1N? 


Xo = Xo ($,Nos%as- ee hNn)- (28) 
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Comparing the integrals in (25) and (27), we see that the function p(s, z,) for which the integral in (25) 
assumes its minimum is defined by the formula 


p(s, 2) = Xo (s, Nor Nas+++s Nn), (29) 


N= D2 (Ss), "= >) Aunty (h=4,...,N). (30) 


If we replace the variable z,, in Formulas (29) and (30) by the random variable Z,,, we obtain, on the basis 
of Formula (20), the following representation for the optimal operator: 


AZ = x,(s, Hy, Hi, ..., Hy), (31) 
where the random variables Hy, Hy, . . . ,Hy are linear combinations of the random variables Z,,: 


H, = > Zvyy (s) = bY yy (s) QZ (2), (32) 
Aa = SenZy = YanQZ(t) (h=1,...,). (33) 


Thus, the optimal operator does not depend directly on the individual random variables Z,,, but only on 
N +1 linear combinations of them. 


Introducing the linear operators 


A® = Sy, (s) 2”, A” = Yana G&=1,...,, (34) 


v 
we can rewrite Formulas (32) and (33) in the form 
H, = AZ (s) @=0,1,..., 0). (35) 
Formula (31) for the optimal operator then takes the form 


AZ =y,(s,A%Z, AMZ, ... , AMZ). (36) 


Thus, if for any possible values of the variables s, No, 14,..., 1) there exists a value Xo of the variable 
x for which the integral I(s, 09, 13, . ~~ » 2)y.X ), considered as a function of x, attains a minimum, i.e., if for 
all values of s, No, Nas eee o MN X 


T(8, Nor Tar ++ » NN» Xo) Sl (8, Nor Na» --+ >» nr xX) (37) 


then the problem of determining the optimal operator reduces to finding linear operators ate) " aft), ae 
AW) and a function Xo, the value of which for given values of s , No, 14,...»%y is determined as the value 
of the variable x minimizing the integral I, defined by Formula (27), * 





* In practical problems it may frequently be convenient, in looking for the minima of the integral in (27) and 
the analogous integral in (58) , to make the following change of variable: 


N 
us > Up, @, (8) + No (A) 


h=1 
using u instead of one of the variable u,, for example, instead of uy (of course, it is necessary for this that 
wy 6) # 0, which is always possible to achieve by a proper renumbering of the random variable U;,). Then the 
integral over the variables y,..., Uy., can be computed and, setting 


(Footnote continued on bottom of following page) 
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Formulas (34) define the linear operators 4%), in the general case, in the form of infinite series, However, 


in some particular cases it is possible to find expressions for these operators in finite form, It is easily seen that 
the linear operators Ath) satisfy the equations 


AK, (t, u) ate K yx (s, u) (ué ge (38) 


AM K .(t,u) =on(u) (VET, h=4,..., N). (39) 


One may convince oneself of this by a direct substitution, taking Formula (14) into account (see also [1, 2)). 
Therefore, in all cases when it is possible to find solutions to Equations (38) and (39) in finite form, it is possible 
to avoid recourse to the definition of the operators ath) by means of the infinite series in (34), 


Formulas (18) and (26) determine the functions uy,(s) and the quantities a,, in the form of infinite series, 


These quantities may also be expressed in finite form by means of the operators A‘), if Formulas (12) and (34) 
are used: 


p(s) =n (s) — AM, (h=1,..., N), 


Agq = A, = AQ, (p, gi, :..', My). (40) 


Thus, although the method presented above for determining the optimal operator was obtained by means 
of the method of canonical decomposition which, in the most general case, gives an expression for the linear 
operators A“) in the form of the infinite series in (34), for the practical application of this method it is possible 
to avoid the use of canonical decompositions if one knows the solutions of Equations (38) and (39) in finite form, 


To estimate the quality of the optimal system, and to compare various real systems with it, it is necessary 
to be able to compute the mathematical expectation of the function r(W,W®*) for the optimal operator. Since 
Formula (31) defines the optimal operator as a definite function of a finite number of random variables, this pro- 
blem is, in principle, quite simply solved, Specifically 


N 
= M[r(W,AZ)}=M |? ( bs U neon (s) + Hy, x0 (8, Hy, Hy, --- , Hy)) = 
a= 
wt ig > a 
om \ oon \ } (Up) du,.. . duy \ ose \ (> (8) + ia Xo 6 Mo arses) 
—o —oo —oo —o =i 
x Io(Nor Ma» +--+ » NN | Ur) dyodn, . . . dnp, (41) 


where f2(n9, 4,...+»%N | Up) is the conditional probability density of the random variable H, with respect to 
the quantity U,. Since the random variable H, is a linear function of the random variables Z,, whose conditional 


probability distribution is normal, the conditional probability distribution of the random variable H; 1s also normal, 
and, consequently, 


q (U, Nor Tay +++ » Tyr 8) = 
@ @ N N 
° ‘ ; 
= \ eee \ f (uy) exp > uN, — > Zs oan du, ...dUy_,, = 
Tor) ire hea P.q=1 


where the variable uy should be replaced by its expression from Equation (A), we obtain the integral in (27) in 
the form 


T (8, Nov M+ ++ > Ay XD) = \ r (u, X) g (Us Nor May « «+ » Nyy $) du. (Cc) 
J 


**Ky x6, u) denotes the cross-correlation function of the random functions Y and X, 
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Ta(Mor Mas +++ » Qn | Ua) = 


‘ { 0 No — My... 1N — My 
= OX DP I -——— as MR, covsveeestecserccoce 
V @nyX*1) K| P 2| XK | No — Mo (42) 





where the conditional mathematical expectation m, andthe conditional correlation moment kpq of the ran- 
dom variables H, are defined, on the basis of the general theory of linear transformations of random functions, 
by the formulas [1] 


N N 
my = M {AZ | uy) = Sy und on = Samu (0 =0,4,..., 9), (43) 
h=1 h=t 


kpq = M ((A°Z — my) (AZ — mg) | up] = 
= AP AP K, (t, t’) (p,q = 0,1, eee, N), = 


and | K| is the determinant of the correlation matrix K = || Kpq||.* 


In the special case when r is a function of the current value of the difference between the effective 
signal and its estimate, and the random function Y is identically zero (i.e., when there are no nonregular random 
components in the effective signal), the results obtained above imply the previously known results obtained by 
Laning in a more complex heuristic way (6). 


The method presented gives an optimal operator only in cases when all the functions gp can be represented 
by a decomposition of type (14) by the coordinate functions x, . We now consider the case when certain of the 
functions yp, for example, 3, ... » Pk» are not represented by a decomposition of type (14). Then, setting 


Cn (t) = en (t) — S) Deen, (t) (h=t,..., (45) 
we obtain from Formula (13) 
a (6) = P.) Ungn (t) + X (t) — bs Uila(t) = Z(t)— > Onn (t). (46) 


Thus, the observed random function Z is completely determined in this case by giving the random variables 
Uy, . . . » Uy, and all the random variables Z,,, The random variables Z,,, in accordance with their definition 
(11) are, in their turn, completely determined by giving the function Z, We show that, in this case also, the 
random variables U,,..., U,, are completely defined by giving the function Z, For this we define the linear 
functions w{P) which, together with the functions ¢ ,, satisfy the biorthogonality condition: 


Yr bo (p.g=ts.--» A) (47) 


This can always be done, for example, by means of the method presented in[1, 2). After this we define 
new linear functionals 


o'?) — yl?) 4 b a” ty as 
p=i,..., & 
2 (48) 


such that they wili be orthogonal to all the coordinate functions x,,: 


@2,=0 (p=i,..., kh). (49) 


® It follows from Formulas (44) and Equation (39) that kp =4pq (Pp, d=1,..+,. ND. 








Since, as a consequence of (12), we have, for all vy, 


om, =0 (h = 1. oat Oh 


(50) 
the functionals @ P) will also satisfy the biorthogonality condition: 
oc, = bnq (p,@ wit, os. s k). (51) 
On the basis of (46), (49), and (51), we have 
k 
0= OZ (t)— SY UiO%, = OZ (t) — Up, (52) 
n=1 
whence 
Up = OZ (t) (p=t, ...) > (53) 


which was to be proved, 


Now, in the case being considered, giving the realization of the observed random function Z is completely 
equivalent to giving the random variables Uy,...,U, and Z,,, and 


M [r (W, AZ) | Z| = M[r(W, AZ) | Uy, ... , Un Bl = 
co co k N 
== \ > \ r( U 1p + > Un@p + Dd) Zvyvs (8, U,, _..Unas)) x 
anit anne he=l h=k+i v 
x J Ghen s+ 5 Oe 1 Oy, ins see Zy) duiz4,...duy, ata 


where f (uy,1,--+»UN | Use~--» Uk, Z,,) is the conditional probability density of the random variables 


Uxsyee+-,» Un with respect to Uy,..., Uy, Z,. This conditional probability density is defined by the 
formula 


f (up) fr (, | 4%) 





{ (Urq «++ » UN | Wy, .-. , te, WSS ™ . 
\ ii | Ftp) fa (yf ty) ity gy «dy (55) 


Thus, using Expression (23) for the conditional probability density of the random variables Z,,, as before, 
we obtain 


I (Ursa, coe 9 Un | Uy, eee 9 Up, Zy) == 
N 4 N 
= Pas + +4444) f(a)exp{ Sun Syeunts— 3 Sy east (s6) 
h=k+1 v Pq=i 


where p(uy, ... , Uy, Z,,) is a function of the variables uy, ... , Uj, Zy,, the concrete form of which is irrelevant 
to our problem, 


By substituting Expression (56) in Formula (54), we come to the conclusion that, to guarantee a minimum 
condition mathematical expectation of the function r, there must be found a function p(s, ty, «+» Uj» Zy) 
for which the integral 


a) oo N 
h= \ es \ r( Upn@, (Ss) + > zy (5), [2 (S, Uy, - ++ 9 Uns 2) x 
anti: —o h=1 v 
N gud 
x f (up) exp| > ue Dd) ney — > > Anata} duy4,-.- dun wn 
hek+1 ov P.a=1 


has the minimum value, 
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Let X9 be the value of the variable x for which the integral 


@ ros) N 
I, (8, %, Tk+is «+2 » Ny Uy, ++ Ux, X) = \ ae \ r( union + "or x) x 
—oa —o Awl 
N » 
x f (tp) exp} > Untn — = > Anas) du,4,...duy, (58) 
h=k+1 P.q=1 


considered as a function of X for fixed values of s, No, Mags +++ TN Ute- ++ » Uy, has a minimum value, 
Comparing Formulas (57) and (58), we see that the function y for which the integral in (57) attains a minimum 
is defined by the formula 


J (S, Wy, «~~ » Uny 2) =Xo(S) Nor Wry -++ » Uy Teta, ~~~ , Ny), (59) 


where the quantities n, are defined by Formula (30), If we now replace the variables uw, and n, by the random 
variables U}, and H;, and if we take (35) and (53) into account, we obtain the following formula for the optimal 


operator: 
AZ = y,(s, AZ, oz, Bitty ! oz, A*tDZ, ii AZ), 


(60) 
In the given case, the mathematical expectation of the function: is expressed by the formula 
co) © © Co) N 
p= Mir(W, AZ) = | ---\ f(un)duy...duy | --- { r( wn6n (8) + No» 
—c —@o —o@ —o@ Awl 
Yo (Ss Nor Way +++ Uk» Metis ees ww) Ta (Mor Nk+a «++ s4N | Un) @Hgdyx 41. +. dyn, (61) 


where fo (Mo, 1k43++++»%N | Up) is the conditional probability density of the random variables Hp, Hy+), 
. ++» Hy with respect to the quantity U;,, which is expressed by a formula analogous to (42), 


Since, in all the preceding calculations the variables played the role of a parameter, and the com- 
putations of the method presented must be carried out for each given value of s separately, then the theory 
presented is completely applicable to those cases when the quantities U}, are arbitrary random functions of the 
variables. Thanks to this, the theory we have presented is also applicable to those cases in which the canonical 
decomposition (7) of the random function Y contains an infinite set of random variables V,,, which do not 
enter into the canonical decomposition of the random function X. To apply the theory to such cases, it suffices 
to separate from the canonical decomposition of the random function Y the set of all terms containing variables 
V,, which do not enter into the decomposition of the random function X, This sum will be distributed normally 
and is independent of the sum of the remaining terms in the decomposition of the random function Y, It is 
therefore possible to denote it by Uj, (s) and thus to transfer it into the order of random variables Up. 


The method given for determining the optimal operator is quickly generalized to the case when the 
effective signal W is an arbitrary nonlinear function of the random variables U;, and the random function Y. 
For this it suffices, in all the preceding calculations, to replace the linear function (4) of the variables U;, and 
Y by an arbitrary nonlinear function of these variables: 


W (s) = (s,U,, ... , Uw, Y). (62) 


In the case when the signal W is a vector random function, and the function r depends on all the com- 
ponents of the effective signal and its estimate, the variable x and the function Xq in all the previous calcula~ 
tions are vectors of the same dimensionality, Formula (42) is correspondingly altered to give the proper normal 
conditional probability density of the random variables H,. For seeking the minima of the integrals in (27) or 
(58) in this case, it is possible to apply the method of steepest descent. 





al 


al 


3. Solution of the Problem for the Criterion of Minimum Mathematical 





Expectation of a Functional of the Signal and its Estimate 








In the previous section, the region T in which the argumentt varied, and in which the random function 
Z is observed, could depend on the value of the variable s, In that case, the optimal operator minimized the 
mathematical expectation of the function: for each given value of s, Here, we consider the case whenr isa 
functional on the effective signal W (s) and its estimate W*(s), defined over a certain region S$ of variation of 
s , assuming that the region T does not depend on the variable s, but may depend only on the region S, More~ 
over, we will consider that the random variables U}, do not depend on s, 


It is obvious that, if Y (s) = 0, then excluding from the computations of the previous section the variables 
1» and Hy, we can completely extend these computations to the case when r is an arbitrary functional on the 
effective signal and its estimate, In this case, however, the integrals in (27) and (58) will not depend ons, and 
will be functionals on the function X(s), In this case, the function xq is defined as the function for which the 
integral in (27) or (58) attains its minimum value, 


If the random function Y is not identically zero, then the method of the previous section is not applicable 
for arbitrary functionals r.Inorder to extend the method of the previous section to this case, we must limit the 
class of functionals r, To be precise, we shall consider here only those functionals r which are functionals on 
some function &(W, W*) of the current values of the effective signal and its estimate: 


r=r(&(W(s), W"(s))), = 
where the function € possesses the property 


E(x, y) =§ (x + A, A(y, A)) (64) 


for arbitrary x, y, and A. In other words, the function € is such that the surface z= € (x, y) is formed by ad- 
vancing a shifting curve, parallel to the xz plane, a certain point of which describes a curve parallel to the xy 
plane, 


In this case 


N N 
(2 Un, + No, U(S, +) =r ( (> UA, x)) ’ (65) 


h=1 h=1 
where the quantity x is defined by the equation 
H(S, 2) = A(x, mp). (66) 


Substituting Expression (65) in Formula (25), we are led to the conclusion that the problem of determining 
the function y (s, z,,) reduces to the determination of the function x (s) =X» G, 13,...>+ 1) for which the 
integral 


Fras <= os 0a) = Yoo J (8 (3 man) x6) ) x 


im * =a h=1 
N ie 
x J (Up) exp > Un, — > > anatst| du, ...duy (67) 
. hel Pq=1 


has a minimum value for the given fixed values of the variables ny,..., N- In solving this problem for all 
possible values of the variables ny, ... » yy» and replacing the variables No, ..., "yy by the random variables 
Ho, .. + » Hy, we determine, on the basis of (20) and (66), the optimal operator 


AZ = (yx, (8, APZ, ..., AMZ), AMZ). (68) 








Analogously to the case when the functions g,..., @, are not expressed by the decomposition of (14) 
by coordinate functions, we shall have, instead of the integral of (67), the integral 


aid co N 
Ty(uyy +. 5 Wey Mpa. 5 Iw, (8) = \ med (> (e/ Dd) une (5), x(s)) x 
—oo —o h=t 
N , N 
* 1 (un) exp > Untn — ~>- > a (69) 
hewk+1 P.q=1 


and the optimal operator is determined from the same Formula (68), in which the linear operators att Pay! iat 
A‘) must be replaced by the linear functionals of)... . a(®) 


The mathematical expectation of the functional r for the optimal operator is expressed, in this case, by 


N N ' 
the same Formulas (41) and (61), in which * ( >, UnOn + No» 10) is replaced by 7 (§ ( > UnMh; ts) 
h=1 h=1 


and, moreover, the variable , is excluded by integrating with respect to it, In accordance with this, the variable 
No is also excluded from Expression (42) for the conditional probability density of the random variables H,, 


The method presented is employed, in particular, when € is an arbitrary linear function of the current 
value of the effective signal and its estimate. 


To find the minima of the functionals in (67) or (69) one may use, for example, the method of steepest 
descent, 


We mention in conclusion that the method presented in this section is applicable only in those cases when, 
in the decomposition of the random function Y, (7), there are no random variables V,, which do not also enter 
into the decomposition of the random function X (or when there are a finite number of such variables). 


4. Determining Optimal Operators of Special Forms 





In certain cases, it is necessary to find an optimal operator, not in the class of all possible operators, but 
in a narrower class of operators, for example, in the class of linear operators, Since, as shown in the previous 
sections, the optimal among all possible operators depends on the random variables Z,, in depending on the H, 
which are linear combinations of them, it is then efficacious, in the given case, to search for the optimal operator 
of the special type in the form of a definite function of the variables H;, i.e., in the form 


AZ =p(s, AMZ, AZ, ... , A™Z, r)), (70) 


where yp (s, 7 dy ) is a definite function of the variables s and 7, and of the parameters ), . In this case, 
the problem of determining the optimal operator of the given special form reduces to finding the values of the 
parameters ) , for which the mathematical expectation of the function r has a minimum value. In particular, 


the optimal linear operator might be sought in the form 
AZ = AZ + ,AMZ +--+ WAZ, (11) 


where the coefficients ), ... A y must be so chosen that the mathematical expectation of the function r is 
a minimum. 


If we replace the function X» by the function » in Formula (41), we obtain the dependence of the mathe- 
matical expectation of the function r on the parameters ), for operators of the form (70): 


p= M [r(W, AZ)| = 


to) © oo os) _N 
= \ oe. \ f (Un) du, ...duy \ are \ r( D Unen + thor H (Ss Tv» »)) x 
—o —o —o —-~» “h=l 
Xfo(nor Mar +++ + Nw | Un) ayody,..dgn. (72) 





The problem of determining the optimal operator of the form of (70) then reduces to determining the 
values of the parameters ), for which the quantity p , considered as a function of the parameters ), , attains 
a minimum, This problem may be solved numerically by means, for example, of the method of steepest descent, 


In cases when the functions gy,..., 9, cannot be represented by a decomposition of type (14), the 
operators A)... | AW in Formula (70) are replaced by the linear functionals © ee gh), With this, 
the variables n},..., 1, and the integration over them in Formula (72) will be lacking. 


5. The Case When the Optimal Operator Can Be Random 





If, for arbitrary values of the variables No, 04,..., 9 (or No, Ug, - ++» Ups VMketeeees nN), there 
exists a value Xo of the variable x which minimizes the integral in (27) [or in (58)]) of Section 2, or a function 
Xo 4s) minimizing the integral in (67) [or in (69)) of Section 3, then this quantity X», or the function X9(s), de- 
termines, as we saw, the optimal operator A which possesses such properties that, for an arbitrary realization z 
of the observed random function Z, and for an arbitrary operator B, the following inequality holds: 


M [r(W, Az)|z] < M{[r(W, Bz)|z). (73) 


It follows from this that in the case when the optimal operator among all possible operators is determined 
also for arbitrary random operators, the magnitude of the mathematical expectation of the function: will not 
be less than for the operator A determined in Section 2 or Section 3, Thus, if the class of operators among which 
the optimal operator is to be found is extended by allowing all possible random operators to lie in it, no im- 
provement of the estimate of the effective signal is thereby brought about, 


Ifr is a function of the current values of the effective signal and its estimate, and the quantity X», mini~ 
mizing the integral in (27) or in (58), does not exist for all possible values of No, 24, . « . » Nyy (or, respectively, 
Of Noy Uy, » » » » Uke Mk 4ge » « » » MQ) but only for some of them, then the use of random operators might improve 
the estimate of the effective signal, 


The problem of determining the optimal random operator leads, in this case, to seeking the random prob- 
ability density 65( x 48 | Mo. %y..-+> 1») which minimizes, for given values of No, 14, . ++ » My, the mathe- 
matical expectation of the integral in (27): 


\ J (8, or Mas = +++ Ns X)8 (Xs | Mor Ma» «+» » MN) ay. (74) 


—co 


In this case, the optimal operator is defined by the formula 
AZ = X(s, A°Z, A™Z,..., AZ), (75) 


where X(s, Mo, 14,065 1») is a random function of the variables $, No, 94» + + « » Nye the one-dimensional 
probability density of which is §(x,8 | No. %.~.-~» My 


In the cases where the functions gy, ... , 9, cannot be represented by a decomposition of type (14), the 
variables 1y,...,}, in Formula (74) are replaced by the variables “ onde Ms and the linear operators at), 
omen Ak) are replaced in Formula (75) by the linear functionals o@ seoer® k). 


} 


6, A Supplement on the Determination of Optimal One-Dimensional and 








The general method presented for determining optimal operators has great generality and is applicable to 
very diverse problems, since the arguments t ands of the random functions can be arbitrary scalar or vector 
quantities, In particular, the method presented here is applicable to problems of interpolating and extrapolating 
scalar and vector random functions of several arguments which arise, for example, in connection with annual 
forecasts, The method given is also applicable to the design of one-dimensional and multidimensional automatic 
systems, 
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If the arguments t and s are different moments of time, then the linear functionals a ) and the linear 
operators A) are expressed, respectively, by the formulas [2] 


A”, = \a, (t) z(t) dt, 


(76) 
. 
AM: = \ git) a 
2 \¢ (s,t)z(t)dt (h=0,1,...,N), crn 
and Equations (38) and (39) become integral equations of the first type: 
\e (s,t) Kx(t, u)dt = Ky. (s,u) (weé7), (78) 
T 
h r —_ ° 
\« (8, t)Kx(t, u)dt = @n(u) (wT; h=1,...,). ae 


The solution of these integral equations, in the most general case, is given by Formulas (34), obtained by 
the method of canonical decomposition [1, 2), However, in particular cases, solutions of Equations (78) and (79) 
may be obtained in finite form by other methods, for example, by the method of integral canonical decomposition 
[4] or by other well-known methods, Determining the linear operators A‘h) and the function x4 minimizing the 
integral in (27) [(58), (67), (69)], we obtain the operator of the desired optimal system from Formula (36) which, 
in the given case, takes the form 


AZ = xo(s,{ g (s,t)Z(t)dt,..., ( g%) (s, t)Z (t) de). (80) 
T T 


Ift is a set of certain moments of time and the ordinal numbers of the components of the input vector 
function of the system, ands is a set of certain other moments of time and the ordinal numbers of the output 
vector function of the system, then the linear functionals a) and the linear operators A‘) are expressed by 
the formulas [1, 2): 


Os = FU ag (t) 29 (t) ae, (81) 
qui T 
(As), = >: | g0(s, t)tq(t)dt (p=i,...,m) 
q-1T 
n (82) 
A™; = > { gi) (s, t)zq(t)dt (h=4,..., N), 
q-1T 
and Equations (38) and (39) take the form of systems of integral equations of the first type: 
* T 
Dif elo) Kan(t, w)dt = Kpi(e, wy (“ENTE ) = 
q~-1T p=wi,...,m 
n 
’ oa 
bY Vas t) Kai (t,u) dt = @,, (u) («er + to) ; (g4) 
q=1 hk=1, ..., N 


In the most general case, the solutions of these integral equations are given by Formulas (34), obtained by 
means of the method of canonical decomposition [1, 2]. Other methods of solving the systems of integral equations 
in (83) and (84) are still not known, Although there was given, in [4], a general formula for this case, adduced by 
the method of integral canonical decomposition of random functions, its practical application is complicated by 
the fact that there still does not exist a method of obtaining the integral canonical decomposition of an arbitrary 



























random function, particularly a vector function, Determining the linear operators A) and the vector function 
Xo minimizing the integral in (27) [(58), (67), (69)], we find the operator of the desired optimal multidimen- 
sional system by Formula (36) which, in this case, takes the form 


(AZ}» = xoo( 8 3 aa (0,1) Zot) at, -... Df ahtn¢s t) q(t) at 


q-1T q=-1T 


(p=1,..., m). (85) 
Example 1, To determine the optimal system for the extrapolating linear function 
W (s) = U,+ Usgs (86) 


according to the criterion of a minimum probability that the error exceed a given magnitude a, with the con- 
ditions that the random variables U, and U3, and the noise X(t), which is independent of them, all be distributed 
normally, and that the correlation function of the noise be defined by the formula 


K , (t, u) =k, (¢—u) = e~ 4 | (87) 


In the given case, the functionr is defined by the formulas 


1 for |w—w*|>a, 
r(w,w*)= | (88) 
0 for |w— w* | <a, 
and the integral in (27) takes the form 
I (s, > Ne» Y=C \\ exp {urn + uN — . (ey; a ay) ut _ 
ju,-+-t498—X] >a (89) 


i 
—(Ci2 + 212) Uyug — > (C22 + 422) us } dujdu,, 


where C is some constant and the c,,, are elements of the matrix inverse to the correlation matrix of the random 
vector (Uy, U,). The integral equations of (79), defining the weight functions of the linear operators A@) and 
A(?) in the given case take the form 


T 
\ g™ (s, t) ec “Ide = 4, 
O<u<7). (90) 


0 
ry 
\ g (s, t) ede = u. 
o 


It is easy to convince oneself, by direct substitutiop, that the solutions of these integral equations are 
given by the formulas 


g™ (s, N= F+F8H4+380—7), 


(91) 
i 
B® (6, t) = Ft — 5 (8) — (1 + eT) Bt —T)]. 


It is obvious that the integral in (89) is a quantity proportional to the probability of a random point, governed 
by a two-dimensional normal distribution, falling in the region defined by the inequality underneath the integral 
signs, In order that this probability be a minimum, it is necessary so to choose the value of the variable x that 
the center of scattering lies in the middle of the strip, i.e., on the line 


Uy + ues — Xo = O (92) 





in (uy, ug) coordinates, Denoting by by and bz the coordinates of the center of scattering, i.e,, the values of the 
variables uy and uz for which the integrand in (89) is a maximum, we obtain 


Xo = bi + bys. (93) 


if we equate to zero the partial derivatives of the exponent in (89) for uy = by and ug = bg, we obtain the 
following linear algebraic equations for the determination of by and be: 


(Cry + @41) by + (Cig + @42) bg = Mm, 
(€x2 + @y2) by + (C22 + @y9) by = Np. (94) 


The quantities by and bz which are defined by these equations are linear functions of the variables 1, 
and Ng, Consequently, the quantity x» , defined by Formula (93), is also a linear function of ny and 19: 


Xo (S, May Ne) = Arm + Ao, (95) 


and Formula (80) for the optimal operator takes the form 


T 
AZ = \ [rg (s, t) + dog (8, 1) Z(t) de, (96) 


0 


where the weight functions g“) and g‘?) are defined by Formula (91). 


In order to compare the operator in (96) with the operator which would be optimal with respect to the 
criterion of minimum mean square error, we substitute in Formula (95) the expressions from (94) for the quantities 
11 and 72, and compare the coefficients for by and by with the corresponding coefficients in (93), As a result, we 
obtain the following linear algebraic equations for determining )4 and \2: 


(C1 + a1) Aa + (C12 + ie) Ag = 1, (Cre + aye) Ax + (C29 + aoe) Ag = 8. (97) 


These equations coincide completely with those which define )4 and ) 2 when the criterion of minimum 
mean square error is used[1, 2], Thus, the optimal operator in the given case coincides with that obtained by 
using the criterion of minimum mean square error, 


Example 2, With the conditions of the previous example, to find the optimal kg using the criterion 


-_ - e 
of minimum mathematical expectation of the function r(w, w*)=1—e w= w , assuming that the random 


variables Uy and U, are independent, where U, is distributed uniformly in the interval | w| <a, and Us is dis- 
tributed normally, 


In this case, the integral in (27) takes the form 
a uw 
I(s, ™m, Na x) = C \ du, \ [1 — e—*(rts—2)") 5, 
—a —@ 
. (98) 


U, 11+ Ughy— Ftp ag Qh tty — +(em+5) u, 
du,duz, 


xe 


where C is some constant and D is the dispersion of the random variable U,, To evaluate this integral we apply 
well-known formulas, obtaining 


nt gy 1 
Homma d—m—m| © (<vza-+ V2 (sux 5 "— sim )) + 


7” f2 1 32 
+0 aV 2a —/ re Box — 7 ua tm exp |— gee x? — (99) 


= Be m+ (1s am 2PaPr) xn} , 








where Kg and x, are functions of s, 14 and 1 which are not essential for our problem, 
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and  @ is the well-known function 
u 1 
i ¢--* 
(D = 2 dz. 
"75 \e * (101) 
Designating by v (p, q, v) the value of x for which the function 
w (x, p, q, v) =e P*-°* (@ (q + 2) + O(q — 2)] (102) 


assumes a maximum, we determine the value x 9 of the variable x for which Expression (99) attains a minimum: 


1 By v a Yor — BF — 4 Ja 
Xo = 28, m+ 5, mts oy fate aV 2a, g |/ .% [ron+@Brr4—Borsn ) ° (103) 


Then, Formula (80) for the optimal operator takes the form 


T -- T 





AZ a Vis t) Z(t) dt +e 3(S aV 2a, elt nzat), (104) 
where 
w) (s, t) = sa g')(s, t)4- g') (s, t), 
w') (s, t) = : / a [vog? (s, t)+-(281%9 — Bors) a (s, #)}. aan 


As we see, in this case the optimal operator is nonlinear, 


Example 3, To find the optimal operator for the general case, whenr is an arbitrary function of the 
current value of the difference between the effective signal and its estimate, and the random variables Uj), in 
(3) and (4) are distributed normally, 


In Example 1, a particular case of the problem under consideration, we found that the optimal operator was 
linear and coincided with the operator determined by using the criterion of minimum mean square error, If, in 
Example 2, we had assumed that the variables U, and U; had a normal joint distribution, then in Example 2, we 
would also have obtained a linear optimal operator, We now show that, in the general case, for an arbitrary 
function r, depending only on the current value of the difference between the effective signal and its estimate, 
the optimal operator is linear if only the random vector (Uy, . . . , Uyy) is distributed normally, 


Let (mj, ..., My) be the mathematical expectation of the random vector (Uy, . . . , Uy), K = || Koqll 
its correlation matrix, and C = || cpq || the inverse matrix, C = K™ 1. Then, the probability density of the random 
vector (U;, . . . , Ux) is given by the formula 


ToT" a 
f (u,) = VY or exp | > 2a Coq (u,, — my) (ug— ma) ‘ (106) 


P.q~i 





where |C| is the determinant of the matrix C, and the integral in (27) takes the form 
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Taking the argument of the function r as a new variable of integration instead of one of the variables u, 
+ + + » Uys WE May carry out the integraticn over the remaining variables up, by well-known formulas for in- 
tegrals of exponential functions with quadratic polynomials in the exponents.* As a result, the integral in (107) 
goes over into the form of a simple single integral. However, in the given case, this transformation is effected 
much more simply by using purely probabilistic methods, It suffices to note that the integral in (107) is pro- 
portional to the mathematical expectation of the random variable r(U), where the random variable U is a linear 
combination 


N 


U = p a W,,%), (8) + M9 — x (108) 
h=] 


of the normally distributed random variables W;, ..., Wy, the joint probability density of which is defined by 
the formula 


N N 


. 1 
f° (1, ... , Wy) = C exp > WY, — > > Gq" pa — 
h=1 P.a=1 
{ N 
i i DY Coq (Wy — My) (w, —™m,) : (109) 
P, q~l1 


where C is a nonessential constant depending on the quantities 1y,...,%). Since the variable U is also dis- 
tributed normally, and its mathematical expectation and dispersion are easily found from the known mathematical 
expectation and correlation matrix of the random vector (Wy,... , Wy) [1], it is not difficult to write the ex- 
pression for the probability density of the random variable U, and to express the mathematical expectation of the 
random variable r(U) by a simple integral. 


The mathematical expectations of the random variables W,, . . . , Wy can be determined as the values of 
the variables wy,..., Wy for which the probability density in (109) is a maximum, As the result, we obtain 


N 


N 
n, = M (W,I = p> lap (1, + > epqlq] (A a= i,..., WN), (110) 


p=1 q=1 


peri : Unp|| ts the correlation matrix of the random vector (Wy, . .. , WN), inverse to the matrix || a 
| The mathematical expectation and dispersion D of the random variable U, in accordance ory well- 
pie = cineen in the theory of probability, are defined by the formulas 


N N 
z= M{[U|= > ny, (8)+%—y7, D==D[U|= ys Lng p (8) @, (8). (111) 


h~t p.q~=) 


The mathematical expectation of the random variable r(U) and, consequently, the integral in (107), are 
nroportional to the function 
ed _ (u—x)* 
viz)= \ rile 2 du. (112) 








* See footnote on p, 518. 
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This function depends not only on the mathematical expectation of the random variable U, but also on 
the variable s, of which the dispersion D is a function, However, this dependence is not essential for us and, 
therefore, we shall not express it in explicit form. 


Let X» be the value of x for which the integral in (112) attains a minimum, Then, the value X, of the 
variable X for which the integral in (112) and, consequently, the integral in (107), attain minima, is determined, 
based on (111), by the formula 


N 
Xo = Not > AnMn + Hs 


rach (113) 
where, for brevity, we have set 
N 
. 
A» = Ln p™p (8) (A =n i,...,), 
p=1 (114) 
N N 
“= > Ap, > Chqi%q — Xo. (115) 
h=3) @=1 
In accordance with (36), the optimal operator is defined by the formula 
N 
— AM: (h) 
AZ=A 2+ dp A'S + (116) 
-1 


We note that since the matrix || Zpq || is inverse to the matrix || anq + Cpq ||, the quantities )),, defined 
by Formula (114), satisfy the ro of linear algebraic equations 


ee ey ee (p=1,...,N). (17) 


q=1 


These equations coincide exactly with the equations which define the quantities ,}, in the problem of find- 
ing the optimal linear operator in accordance with the criterion of minimum mean square error [1, 2], Thus, in 
the given case, the optimal operator is a nonhomogeneous linear operator coinciding exactly, up to a systematic 
displacement yp, with the linear operator which is optimal when the criterion of minimum mean square error is 
used with normally distributed effective signal and noise, 


We show further that, in the particular case whenr is a nondecreasing function of the absolute magnitude 
of the difference W- W*, the optimal operator we have found completely coincides with that given by using 
the criterion of minimum mean square error, For this we note that, in the case of an even function r, Formula 
(112) can be rewritten in the form 











) _ (u-xy* __ (u+x)* 
v(x) = \ r(ujfe "2 +e ™* Jadu. (118) 
0 
Differentiation of this formula gives 
‘ co _ (u-x)* (u+x)* 
v'(z)= D \ r (u) [(u — az)e => —(u+ z)e 2D ] du, (119) 
9 
from which we find 
v’ (0) = 0. (120) 


Differentiating Formula (119), we find the second derivative of the function y (x): 
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1 _- ——_—— 
v” (x) = 7 r (u) {[u — x)? — Die 2D +|(u + 2)? — Die 2D } du, (121) 
0 
whence 
v” (0) = ml (u?— D)r(uje 2” du = 
0 
o(¢t ae VD ma 
— i \ (u’— D)r(uje 7 du— \ (D—u®)r(u)e 2? au (122) 
VD 0 
Since 
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ak ee 


{ (u? aie Dye ? 
VD 0 


(123) 


and the function: is taken to be nondecreasing for u > 0, the first integral in (122) is always greater than the 
second and, consequently 


¥ (0) >0. (124) 


Thus, the function y has a minimum for x = 0 for any even function r(u), which is nondecreasing for 
u> 0, 


We have thus proved that if the effective signal and the noise are distributed normally, then, for an ar- 
bitrary function r depending only on the current value of the difference W— W*, the optimal, by criterion 
(1), operator of the class of all possible operators for which the criterion is meaningful is a nonhomogeneous 
linear operator, where only the systematic shift jy of Formula (116) depends on the form of the function r, If 
the function r is a nondecreasing function of the variables |W- W*| then the shift » does not depend on the 
form of the function r, and is always obtained the same as when the optimal operator is determined by the 
criterion of minimum mean square error, 
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DYNAMICS OF A RELAY-TYPE ELECTRIC SERVOMECHANISM WITH A LOAD 
VARYING PROPORTIONALLY TO THE MOTION 


N., S&S. Gorskaia 


(Moscow) 


The paper treats of the dynamics (combined oscillations) of a relay-type electric 
servomechanism, the load on which varies proportionally to the motion, 


The motion of the servomechanism considered here is described by a complete 
second-order differential equation,in the right member of which is a relay function with 
loops and a dead zone.* A complete solution to the problem is given by means of the 
method of point transformations [1]. 


1. Introduction 





There exists today a rather large number of works, on both the experimental and the theoretical side, in 
which the auto- oscillations of relay-type servomechanisms are considered, 


Among the works dealing with the dynamics of relay-type servomechanisms when the linear portion is 
described by complete second-order differential equations are works [2-6], However, these works are primarily 
given over to the study of the auto-oscillations in the systems considered, In none of them is there to be found 
a decomposition of parameter space into regions of qualitatively different system behavior, 


Moreover, the works cited dealt with system motions described by second-order differential equations the 
right sides of which were in the form of the simplest relay function, 


The present work investigates the dynamics of electric servomechanisms under certain assumptions which 
lead to the consideration of nonlinear systems, the motions of which are described by complete second-order 
differential equations in which the right sides are relay functions, the characteristics of which contain loops and- 
dead zones. The method of point transformations [1] gives a complete solution to the problems of: 


1) investigating the structure of the decomposition of phase space by trajectories; 


2) finding the decomposition of parameter space into regions of qualitatively different system behavior, 
i.e., into regions corresponding to the presence or absence of servomechanism auto-oscillation; 


3) obtaining analytic expressions for the critical relationships between the servomechanism parameters; 
4) investigating the stability of the periodic solutions, 


2. Posing the Problem, Equations of Motion of the System, 








In practical automatic control problems, widespread use has been made of electrical systems for positioning 
control by relay-type amplifiers, represented in Fig, 1, Used in this system are a dc motor with independent ex- 


* This problem was posed by V, V, Petrov, An approximate determination of the critical relationships between 
the parameters by means of degenerate limiting cycles was considered in the work of B, N, Petrov, V. V, Petrov, 
and N, S; Gorskaia, “Investigation of one- and two-stage relay-type servomechanisms with loads 

portionally to the motion," (Reports of the Institute of Automation and Telemechanics, AN SSSR, 1954), 











citation, 1, controlled by a type RP-5 polarized relay, 2, The output shaft is connected to potentiometer slide 
4; its deviation a, we shall consider equal to zero, Then, the angular discrepancy between the output and 
input shafts equals the angle of deviation a of the input shaft, taken with reversed sign, Potentiometer slide 
5 is connected to the input shaft, 
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The electric motor has a load varying proportionally to the angular rotation of the output shaft of reducer 
3, i.e., the load moment M, = K,.a. 


As is well known, for such a system of an electric motor controlled by a relay amplifier, when the in- 
duction of the motor’s armature winding can be neglected, the following system of equations is applicable: 


discrepancy equation 
_ AU =—U,,, 
(1a) 
relay element equation 
U =¥ (AU) = |U | F{AvU}, (1b) 
motor electric circuit equation 
U=IR+ E,, (1c) 
equations of motion 
dQ , 
Imig = Cul K, 2a Ka, (1d) 
Q er OQ; (le) 


K, 


The following notation was used in the equations: I is the armature current, R is the active impedance of 
the armature winding, U is the voltage impressed on the motor circuit, E, =C_ , is the counter-emf of the motor, 


C, and Cyy are motor constants, Q, and Q = “a are the angular rotational speeds of, respectively, the motor 


shaft and the reducer‘s output shaft, K, is the reduction coefficient, Up, =Kga is the feedback voltage, propor- 
tional to the angular rotation of the input shaft and J,,j is the moment of inertia of the reducer’s output shaft, 


The nonlinear function F(AU), the characteristic of the controlling relay, has the form shown in Fig, 2. 


Eliminating from the system of equations (1a) - (1d) all variables except Up), we obtain 


d*Up, dp 





o J_R KC, sie KR 
A= [OTR K Cc,’ whee |0| Kc,” — [0 | K,K,C,, 





or, 





We introduce the following dimensionless variables: t® = = x ae Upp, and also the dimensionless 


system parameters 
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Fig. 4, Servomechanism blockschematic: 1 is the 
Fic. 3 RP-5 controlling relay, 2 is the electric motor and 3 
8. is a rigid feedback connection, 
Equation (2), written with dimensionless quantities, takes the following form: 
d*z dz 
a? ta +e = —O@). (3) 


The nonlinear function @ (x) is shown in Fig. 3. This function can be written in the form 


1 for t>n—>, if O(z,) = 1, 
O(z)=) 0 for |z|<u+ >, if D(x) =0, 
—1 for z<—(u—F), if D(z) = —1, 


where Xp is the initial condition at the moment of time t® =+ 0, 


Figure 4 gives the block schematic of the electric servomechanism described by Equation (3), Thus, the 
system is characterized by three parameters, y, 1, and A, the load coefficient, the dead zone, and the loop 
width in the characteristic of the controlling element. 


3. Phase Space 





Setting k = y, we write the equations of motion (3) in the form 


z=y, yty+ye=—O(z). (4a) 


The (x, y) phase space, in accordance with Equations (4a) and the characteristic (x), consists of three 
plane sheets, lying one atop the other, 


On sheet I the system motion is described by the differential equations 
. . 4b 
c=y, yt+ytye=0. (a) 
Correspondingly, on sheet II 


z=y,ytytye=—t. (4c) 


The phase plot is symmetric with respect to the origin, since Equations (4b) and (4c) do not change when 
X andy are both replaced by their negatives, 














The singular points of the system of differential equations, (4b) and (4c), have the coordinates (0,0) and 


( 4, » 0} and are stable foci (in the case y> ry ) or stable vertices (in the case y< =) . For 


- o> a ie , the singular pon : 2.) does not lie in its region of definition. 
Y y We investigate the actual practical case when 
_— 
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feet wanes The motion of the representative point on the three- 
DY sheeted phase space occurs in the following way (Fig. 5). 
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In moving on sheet I [@ (x) = 0], the representative point, in 
o12)=1 accordance wigh the equation of motion (4a), moves during 
% a certain interval of time from the half-line 


A 
(z= —pet zy y>0) 
to the half-line 


(2 “h + t , y>0) }- 
It then moves from some point on this line 
(2 = Ye + > »y¥> 0) 
now on sheet II [@ (x) = 1), ats on half- line 
(=n, ¥<0) 
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at some moment of time, then once again moves to sheet I, 
' still does not reach the line 
| 
! 


(2=—w—F y<0), 
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transfers over to sheet II] [6 (x)= — 1], etc. 

It may happen that by its motion on the three~sheeted 
phase space the representative point describes a closed curve, a limiting cycle, mirroring the presence in the 
system of a periodic motion, 


Fig. 5 


4. Point Transformations 





It follows from our consideration of the structure of phase space and the motion of the representative point 
on it that the problem of discovering limiting cycles, and the investigation of the decomposition of phase space 
by trajectories, lead to the investigation of the point transformations of lines into lines, Each such transformation 
consists of two point transformations, or of the corresponding functions: of the function (A), which corresponds to 
the transfer of the representative point on sheet I of the phase space, determining the law transforming the half- 


A A : , 
line (x Oe >, y>0) into the half-line (2 =H+ +, y> 0) , and of the function (B) which 
corresponds to the transfer of the representative point on sheet Il of the phase space, determining the law trans- 
A RS : ; 
forming the half line (z= Ue +5 4 y>0) into the half line (x =—~-.— =, y<0). 
With such a representation, a fixed point of these half-transformations will, obviously, attest to the presence 
of a limiting cycle on the phase space, We find expressions for the corresponding functions, 


Sheet I[@ (x) =0] of Phase Space, At the moment of time t* =0, let the representative point (Fig. 5) 





start its motion from the point (2 =—Uet 5 ,_y= Yo) . With these initial conditions we determine the 
constants of integration in the solution of the system of differential equations (4a) for (x) =0. If, after the 
interval of time t* = = , Where w* = =. — y > 0, the representative point arrives at the point (x =p. + = 
y= ys) , we then obtain from the solution we have found” 
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mt babe Ftd (mt Smet emt Send 
n= te wit y,—(—u. + 4)(4—w*)]she+ wy oh. (5) 


Solving these equations with respect to wu, = “—y and y= -, we obtain a representation of the desired 
function (A) in the following form: 


(u, +5) ; + (u, +) (p sh t + ch *) 








== sht ‘ (6) 
(u—F)e* —(u + F)@psnr— en t) 
agers sht ; 


i 1 ° tT 
Here p=, w= 7—y> 0, t =, 
In the case y > 2 we find, analogously, 


nt fa tem {fy + +(- He + $)| sins + w,(— te + $) eos +}, 


(7) 





n= tem f—ty,— (—n.+4)(4—ut)]sine + wy cost}. 


Solving these equations for ug _—_ and v= ie , we find the expressions for the desired function (A) 
i 1 
in the following form: 





y= [(ue + se" + (ue — >) (p sin = + cos +), 


sin T (8) 
v =.= [(m— se (we + +) (psin« -- cos *)|. 
Here pP=>, wi = —i>0, ha. 


Sheet Il [® (x) =1] of Phase Space, We assume that at the moment of time t* =0, the representative 





A 
point (Fig. 5) lies on the half-line (t=u+ >, y>0) at the point (c=m+ 5, ¥=Ve, y>0), 


With these initial conditions we determine the constants of integration in the solution of the system of differential 
equations (4a) for @ () =1, 


If, after the interval of time t* =? » Where w*= 4. — y > 0, the representative point is found at the 


A 
point 7=—%e—Z, yY = ¥s = —|¥|> then from the solution we have obtained we find 
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eet +—)|sh “1 + u( Beth + —)ch a}, 
ly|= tet ~~ ety + (ue + > + 1) (4 —wt) sh  — wy oh 44) 


* ‘sh'= ‘sinh’; ‘ch’ = ‘cosh’; ‘ctg’ = ‘cot’; ‘tg’ = ‘tan’ — Publisher, 








Solving these equations for u, = dul and v —_ , We obtain a representation of the desired function 
(B) in the following form: 
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’ sh 7; [(u. 
1 1 . T1 
Here P=, W=7-y>0, =~. 
In the case y < 7 we find, analogously, 
ar 1 — 1 A 1\).. 
mot ta e™ nt ou + $+4)]sins t 
A | 
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1 —‘i1 ' A ; 
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Solving these equations for ug = dul and v= a » We find an expression for the desired function (8) 
1 1 
in the following form: 


aire [ (tHe + > + ~\e"™ + (ue — > = ~) (p sin t, — cos a)], 
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The curves given parametrically by Equations (6) [or (8)] and (10) [or (12)] are defined for all positive 
values of r and ry, However, because of the rationale for introducing the point transformations, the following 
rules of “selection” [8] must be satisfied by the branches to be retained out of the set of all curves which could 
be constructed in the first quadrant when r and ry vary from zero to infinity: 


1) for transformation (A), for each value u, > 0, that value of v > 0 is chosen which is defined by the 
least root r of the first equation of System (6) [or (8)]; 


2) for transformation (B), for each value ‘v > 0, that value of ug > 0 is chosen which is defined by the 
least root ry of the second equation of System (10) [or (12)). 


The curves thus singled out for (A) and (B) are given geometrically by the functions uy = uy (v) and u, = 
=ug(v), forming the corresponding point transformations (A) and (B), These curves, constructed in the Cartesian - 
coordinate system, where the quantity v > 0 is laid out along the axis of abscissas and the quantities u, > 0 and 
ug > 0 are laid out along the axis of ordinates, comprise the so-called (u,, uz, v)-diagram., 


An iteration of the point transformations (A) and (B) here corresponds to a broken line, consisting of line 
segments parallel to the coordinate axes, inscribed between the curves for (A) and (B), An intersection of the 
curves on the (ug, ug, v)-diagram bespeaks the presence of a fixed point of the half-transformation, i,e., the 
presence of a limiting cycle on the phase space, and the character of the intersection of these curves is indicative 
of the stability or instability of this limiting cycle.* 

* The criterion for stability (in the narrow sense) of a limiting cycle is that, for the point of intersection of the curves 
Ug = ug (v) and ug = ug(v) corresponding to this limiting cycle, the inequality duj/iv_ ey hold, With this, the 


du ;/dv 
angle of slope for the tangent may be either positive or negative, 
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Below we give various cases which can be treated with the transformations (A) and (8) which we have 
obtained, 


For the case y < 2 the function corresponding to (A) is found from (6) and the function corresponding 


to (B) is found from (10), We investigate the behavior of the curves uy = uy(r), v=v(r), and tg = ug(rq) 
v=v(r4), defined by these expressions [8], For this, we consider the qualitative behavior of Curve (6). The ex- 
pressions for the derivatives may be written in the form 





du; Px v pt dv u —pt 
dt” the © y eee (13) 
Consequently, 
du, 2prt 
oe. tt 
d*u 2 (p? — 1) e®P* A? A? 
re [(ui + 4-)sh pe + (ut—)pshe+ uch px]. ag 


From this it is easy to deduce the properties stemming from these expressions for the function and its de- 


rivatives, It follows from (14) that the derivatives “7 and r in the first quadrant are always positive as 


rT varies from zero to infinity and, consequently, the curve uy =u, (v) remains concave as r varies within the 
stated limits, For r—> 0, the derivative du,/dv tends to unity, The curve wy =u; (v) has the asymptote u, = 
=v+4pule. 


For small r we have 
He A : 
uy = 2— + wep + (ve +) [p> + P* 5 +...| 
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p= He A 2 
v= 2 — pep + (4. — +)|P 
and, fort — 0, the functionsu—> o andv~> o, 


d 
For vy =0, the derivative ~— equals zero, and in the interval 0 < vy <@ the inequality 1 > ya > 0 
holds, 


Thus, as r decreases from some value Tf, for which o =0, to zero, Curve (6) enters the first quadrant 


from some point on the axis of ordinates and thereafter approximates to the asymptote uy =v + 4py,, remaining 
always concave, 


We now consider the qualitative behavior of the curve ug = ug(v), defined by Expression (10), The expressions 
for the derivatives have the form 
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Whi’ %°’ dy chy * se 
Consequently, 
OMe ng tg re 
seal “lien, 
d™u 2 (p® — 1) e8P™ i\}, a 
a us [(u« + 7) | |sh ps + (16) 


+ [(we ++) —42]psha+(u++)dchpal . 


it follows from (16) that, in the fizst quadrant, the derivative ts. is always positive as ry varies from 





d 
Curve (10) has « point of inflection,* For small ry we have 


zero to infinity. The derivative os can be either negative or positive, depending on the value of 14, i.e., 


u, = - —Ap+(m+3 +4)[p 3 —p 3 +...] 


, 


v= —-> —Ap+(u—ZFt+—)[P FZ pa : ee 


and for ry; ~* po, the functions u,— (He —_ > + —)(p— 1), v— oo. 
For some value ry = 7, the function v =0 and the derivative = =0, For r;> 14, the function v > 
dug dv dug 
> 0 and the derivative 0, since ar > 0 and the numerator in dy increases with increasing rj. 
i 


d 
As ry tends to infinity, the derivative a tends io zero, For vy > 0, the derivative - <1, 


Uy, th Thus, the function ug=u,(v) > 0, for vy =0, continues to 
“ 


increase up to the limit (p — 1) (1. _— = + - with in- 


creasing X. 





The derivative = =0 for y=0 and for y=oo, Within 


d 
the interval 0 < v < o the inequality 1 > —_ > 0 holds, 





Thus, as Ty increases from zero to infinity, Curve (10) 
j enters into the first quadrant from some point on the axis of 
ordinates and thereafter approximates to the asymptote 








Fig. 6, Diagram of the point transforma- a 4 
tion, Point 1 corresponds to an unstable ua=(p—1)(p.—s +). 


limiting cycle, and point 2 to a stable 
limiting cycle, Figure 6 shows the diagram of the point transformation 


formed by the curves uy = uy (v) and ug = ug (v) which we have 

found, The qualitative investigation, carried out above, of the 
behavior of these curves shows that Curves (6) and (10) may have one or two points of intersection to within an 
even number of points (i,e., 1+ 2k or 2+ 2k, where k=1, 2, 3,...), It is also possible that there are no points 
of intersection, 


If these curves have two points of intersection, then on the phase space there exists one unstable and one 
stable limiting cycle. These two cycles may fuse into one semistable limiting cycle (Fig, 7, a) and thereafter 
disappear, The case is also possible when only one of the limiting cycles (the unstable one) disappears, and the 
other cycle persists (Fig. 7, b). 


The boundary in (y , A, u¢) space corresponding to this plot (a bifurcating parameter ratio) will separate 
the region of parameters for which two cycles exist from the region of parameters for which one cycle exists. 


The curves uy = uy, (v) and ug = ug (v) satisfy, respectively, the equations 


* Indeed, for small ry the derivative se is positive [Curve (10) is concave], since the expression in curly 


1 A 
brackets in Equation (16) is positive, having the form {(u, + —) A—> prt for small 14, 





du v onl 
= = — e apt . 
dr uy, 4 


At the point of contact of these curves, u(r) =ug(ry), v(r) =v (ry). 


Uy Uy 
























































Fig. 7, Diagram of point transformations: a) is for 


limiting 
cycle 
the case of critical parameter ratio, b) is for the 


case of a bifurcating parameter ratio, re ? 





Then, equating their derivatives, we find that the contact of the curves uy = uy (vy) and ug=ug(v) occur 
only on the line v = 0 (i.e., on yg=0, since v am. ), The boundary in (y, A, H, ) space, exactly defined by 


the analytic expression for the critical parameter ratio, uy (r) = ug(r4), v(r) =0, v(r) =0, divides the region of 
stability from the region of auto-oscillation, However, the construction of this boundary in (y, A, #, ) space is 
rendered difficult because of the transcendence of the equations, 


Figure 7, a gives the diagram of the point transformation in the case of a semistable limiting cycle, 


Thus, the semistable cycle on the (x, y) phase space will pass through the boundary of the dead zone 
(Fig. 8). 


For the case y > 4 the function corresponding to (A) is found from (8), and the function corresponding 


to (B) is found from (12), We investigate the behavior of the curves uy = ug(r), v= v(r) and ug = ug(7y), v= v (ry), 
defined by these expressions. For this, we consider the qualitative behavior of Curve (8), The expressions for the 
derivatives may be written in the form 
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and, consequently, 
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It follows from (18) that, in the first quadrant, the derivatives ~ and io are always positive as r 


varies from zero to infinity and, consequently, the curve uy = uy (v) remains concave as r varies within those 
limits, 


For small r > 0, the expansions of the functions uy = u,(r) and v = v(r) have the form 


2u, A 2u, 
y= “eit Map +(e +>)P>+ soe) V= ——— ep + (ie —$) p*F F.. 


ai. 


and, forr —> , the functionsu;—> mandy» o, The asymptotic curve of (8) is the line uy =v+ 4pue. 


We consider the qualitative behavior of Curve (12), It follows from (12) that, as ry—> 1, the numerators 
of ug(r4) and v(r4) tend, respectively, to the values 


pr 


[mete tjere $44) 


and 


($4) eet $44) 


and the denominators tend to the limits lim sin ry= + 0 for ry—> 170, Therefore, the functions ug(r,) 
and v(rq) tend to positive infinity in the case ry—* 1-0, The asymptotic curve of (12) is the line 


A i —pr A { 
= e~Ppr  <iee Ps mall enaite 
ume > + 2p[ (ue Ft+—)+e (u+t>+—)|. 
For small ry > 0, the expansions of the functions ug(1r,) and v(r4) have the form 
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The expressions for the first and second derivatives have the form 








ee co ee 

dv Ug ’ 
~ rs 2 (p? 3c f— [(v 4 + +4) sin rahe (19) 
+ [+P ]ouns. (ut taeda) 


It follows from (19) that, in the first quadrant, the derivative 1. is always positive as ry varies from 


zero to infinity, The derivative ye may be either negative or positive, depending on the value of rq, that 
is, Curve (12) has at least one point of inflection, * 


It thus follows from an analysis of the function uz = ug(v) that the curve Corresponding to (B) is the branch 
of Curve (12) for 0 < ry < ©, 





* Indeed, for small 1, the derivative fe is positive [Curve (12) is concave}, since the expression in curly 


A? 
brackets in Equation (19) is positive, having the form {(u, + =) A—> prs} » for small 4. 








The qualitative behavior, as investigated above, of the curves u, = uy(v) and ug =u(y), in the case y > 4. 
shows that contact of these curves, satisfying, respectively, the equations 
du, vy 2pr dug 


o a oe 


= —2pt, 
Ug sf 
also occurs just on the line v =0 (i,e,, on ys =0, since v = 72. ). 


The semistable limiting cycle for the case y > = is shown in Fig, 9, 


5. Critical Ratio of the System Parameters 





A, Case of a Large Load, y > v1 





We shall analyze the motion of the representative point on the semistable ‘imiting cycle shown in Fig. 9, 
1, Equations of the integral curves on sheet I 





The equation of motion of System (3) for sheet I will be +k +yx=0, The solution of this equation has 
the form 





f= Z 








1 
1+2—-- 
e ° cos (y 1-7" — arc tg ————* Bia 3} (20) 


where Xp and yo are the initial coordinates of the representative point on sheet I, 
We find the time t*, during which the representative point goes to the point x =x, y =y,=0. 
Using (20), we find 
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Here 








and, so that the equality yp=0 may hold, we must set * 
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i+2. oe AES 
Vi—-+ tin— arc tg ———*- vrai —arctg V4; —1 =>: 








* The equality yg =0 also holds when t*, =a, However, this case is of no interest, 



































Then, the time for the representative point to move to the point (xg, yg) is found from the expression 


Anme e 1-2 3 
ee Fe (21) 
tx.= - 1 
\ onli V 4y —1 





Substituting (21) in (20), and bearing in mind that 


125% 
eos | /y—4 ti, — arc tg ————* VR 


cos (F sang RSP: 


we obtain, after some simple transformations, the expression for x, 











1, 


z, 
443 — 
_=—-2 “+ Zt + Lex =A (5 y )}. 
0 Y Zz »| rather x + arc tg VES (22) 





2, Equations of the integral curves on sheet Ill 


The equation of motion for System (4) for sheet Il will be X +% + yx=1, The solution of this equation 
has the form 





i — ite 


Piss Bo. 4y ( a i 
( 7) Fi ¢ 2 cos V— — arc tg 7). (23) 


where Xz is the initial coordinate of the representative point on sheet III (x = xs, y = ys = 9). 
The final value of the coordinate x on sheet III is found from the expression 


% =(%—~)V fy ie — tes ( VW 1-7 tn — are tg =), (24) 


where t*sk is the time for the representative point to get from the point (xs, ys) to the point (xp, yo). 





(e386 4 e~ iB) 
2 
expand the exponential functions in powers of t*s,, limiting the series to three terms each, we obtain 





If we replace cos § in Equation (24) by where 6 =V Y — Mh ~ ate tan Ae » then 
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After substitution of the values for xp and xs, we get the following form for the expression thus found for 
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We get the expression for the coordinate y, by differentiating Equation (23) with respect to the time t*, 
and substituting the values for xy and t® =t*s; in the expression for the derivative. Taking into account that 


iy 4 RAY : r 4 1 
cos ( /4—4 ton are tg 7 — are tg V4; — 1) sin (tan V 1-4) . 
we can write the desired expression in the form 


een (m+ 5 Te? TES exp (-2) sin (t5q Vi- x). (26) 


3. Critical ratio of the parameters 











The semistable limiting cycle shown in Fig. 9 goes through sheets III and I of the phase surface through 
the points with the coordinates 


A A 
~ i » a a y; = 0, %=—w+sZ 
Yo,as determined from Equation (26), 
A : 
%=ht+s, y= 0. (27) 


Substituting (26) and (27) in (22), we get the critical ratio of the parameters 
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-t® a 
Oe expression thus obtained may be simplified by replacing in it the functions e -. e , and 


~t?sk 





e * by their values determined from Expression (24), in which xg and xs are replaced by their values 


determined from (27), Then, finally, we obtain the critical parameter ratio for (y, #¢, A), in the case 
Y> 4. » in the following form: 
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A= ez? An = 4t(me — F +5) +2(m—F)(u—F +4), 


As = tex VW t— 4, Ag=2(e — 5) (ue —F4+2)VA—E 


A, = 4, —1, 





and the parameter t*s, is defined by Expression (25), 















































Figure 10 gives, for the case y > =. the 


Fig. 10. ( ¥ ,#¢ , A) parameter space for the case 


decomposition of (y , 4 ¢, A) parameter space into 1 
regions corresponding to auto- oscillations of the Y> ~[ + Lis the surface F (y . He» A), II is the region 
servomechanism, and to their absence, of stability, III is the region of auto- oscillation. 


B, Case of a Small Load, y < 4. 





We shall now analyze the motion of the representative point on the semistable limiting cycle shownin 
Fig, 8. 


1, Equations of the integral curves on sheet I 
The expression for the coordinate x on sheet I as a function of time t*, and of the initial conditions Xx, 





and y, has the form 
°= Yot Gato e7aut® + Yo + %i% e~at* | (29) 
gz— M1 qi — Ya 
where 





—u=—>+ V 1-1, 


et a = 














The time t* yd, taken by the representative point to pass from the point (xp, yo) to the point (xy, yy = 0) 
is found by differentiating Equation (29) with respect to the time t*, and expanding the exponential functions 
in series limited to the first three terms of the expansions, 


As the result, we get 
{ 1 
#2 2 \2Yo + >%) e 


4y 
t — - / t —_— 9 = 0, 
id ~~ (— y) yo + 2yz, 39 GY) yo+ 272, 








from which we find the desired time 
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Substituting the value of t*4q from (30) into (29), we obtain the expression for the coordinate x, in terms 
of the coordinates xX» and yg: 


nal ta (s+ Vi— ve (4+ V 5-1) Sd 
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7 2, Equations of the integral curves on sheet III 















































= In the case y < > , the expression for the 
2h coordinate x, in terms of the initial conditions xs, 
1+~W) ‘ ys =0, and of the time t* sq for the representative 
ars point to move from the point (xg, ys) to the point 
(x9, Yo)» has the form 
IK | Rn (-- 4) 
! “ 2 the 
205] | Ni) [1 fig ot eee 
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ay Te ] (»— +) 
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Hil} . Hi f where qy and qg are the same as in (29), 

ay if! >|] tr | f Substituting in (32) the values of x, and x, 

1 Hi h H me h from (27), and of qy and gg from (29), and also ex- 
aeshatn MUI i j panding the exponential functions in series, keeping 
PP HH HHA I only the first three terms of the expansions, we get 
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a= . (33) 
Fig. 11. (y, Me , 4) parameter space for the case (4+) y¥+1 
Y< =. : 1 is the surface F (yu, A), II is the region 


of stability, Ill is the region of auto-oscillation, 
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By differentiating (20) with respect to time t*, and substituting the value t® =t*s4, we obtain the value 
for Yo 
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3, Critical ratio of the parameters 


(34) 


1 
oS 


By analogy with the case y> “4 + We obtain the critical parameter ratio for (y,u¢ , A ) in the case 
by substituting in (31) the value of y, from (34) and the values of xy and xz from (27), As the result, we 
obtain 
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(35) 
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n= 24(—we+F)VI—%, t= 4[(m + $)1 41], 0= VI, 


and the parameter t®sq is defined by Expression (33) 


Figure i1 gives, for the case y < 4 , the decomposition of (y , 4¢,A) parameter space into regions corre - 
sponding toauto-oscillations of the saleable and to their absence, 


SUMMARY 


1, In relay-type servomechanisms, the motions of which are described by complete second-order differential 
equations where the right members are in the form of relay functions whose characteristics contain loops and dead 








tial 





zones, there exist two different critical parameter ratios:* one for large loads, (28), and one for small loads, 
(35). 


2, It follows from the decomposition of (y,A, #_)** parameter space into regions of qualitatively different 
system behavior (Figs, 10 and 11) that an increase in the load, varying proportionally to the motion, aids the 
appearance of auto-oscillation in the system, 


The present work was carried out under the supervision of B, N, Petroy and V, V, Petrov, 
The author of this work wishes to express her gratitude to N, A, Fufaev for his advice and help in the portions 
devoted to obtaining sufficient conditions for stability. 
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. The critical ratio for the case of an average load, y = +. (corresponding to equal roots),is of no interest, 
since, in practice, the characteristic equation of the system has different, unequal roots, 


** Generally, lu +4 |<+ : 




































STABILITY OF PERIODIC CONDITIONS IN AUTOMATIC CONTROL SYSTEMS 
FOUND APPROXIMATELY ON THE BASIS OF A FILTER HYPOTHESIS 


V.A. Taft 


(Moscow) 


An approximate method of analyzing the stability of periodic conditions in 
nonlinear automatic control systems is outlined, 


A brief comparison between the above method and that based on the auto- 
resonance theory is given. 


1. Introduction 





When investigating nonlinear automatic control systems it is often expedient for determining oscillatory 
conditions to use a harmonic balance method [1], for instance in the form of Gol'dfarb's graphic representation 
(2). 


The determination of oscillatory conditions consists in the first place in finding periodic solutions to a 
system of differential equations which represent the process, and in the second place in separating the stable 
periodic solutions from the unstable ones, The necessity of solving both parts of the problem in finding oscillatory 
conditions, applicable to systems of automatic control, is connected with the recently discovered [3 and 4) con- 
tradiction between the basic assumptions on which the harmonic balance methods rest and the dynamic properties 
of automatic control systems, 


The harmonic balance method has been proposed and is valid for equations which contain the small para- 
meter #, so that when ¢ =0 the equation produces a harmonic periodic solution, In such a case, when y is 
sufficiently small, the periodic solution, if it exists at all, differs in its form but little from a harmonic one, a 
circumstance which permits one to find a periodic solution in the form of a harmonic function, In addition the 
setting-up process of a periodic solution, when disturbed but little, differs from the periodic solution itself only 
by slow changes in the amplitude and phase with respect to time. The latter condition validates the stability 
investigation methods used, for instance, in Gol'dfarb's construction [ 2], 


Equations representing the setting-up process do not, as a rule, contain the small parameter y», In order to 
justify the application of Gol'’dfarb'’s construction, other considerations are taken into account: it is assumed that 
the oscillations only contain frequencies whose harmonics lie outside the passband of the linear part of the system, 
this part being considered as a low-pass filter, These considerations suffice to justify the finding of periodic solutions 
in the form of harmonic functions, Hence, they justify the application of Gol'dfarb's construction for determining 
periodic solution in a certain frequency band,® 


These considerations, however, do not provide reasons to believe that in the transient state the amplitude 
and phase of deviations from the steady-state periodic movement change slowly; hence, it is not valid to use the 


e 
In the simplest case of an odd nonlinear characteristic this frequency band is equal to =. 
Wep is the linear filter cutoff frequency. 


<w< Wep» where 








harmonic balance method (in particular Gol'dfarb's construction) for investigating the stability of periodic con- 
ditions in such cases, 





Thus arises the problem raised in [6]: assuming that the system under investigation contains a linear 
filter, to find the conditions of stability for periodic solutions, obtained on the basis of these assumptions, 


An approximate solution of the problem is proposed below, 


2. Equation for Small Deviations and Reduction of the Problem to a 





Characteristic Equation of an Infinite Degree 













Let us examine a system of automatic control, consisting of an arbitrary linear part and a nonlinear 
feedback, The equations of the process will assume the form® 













D(p)X =—K (p)Y, Y =F (X), (1) 





where p= =, D and K are polynomials with constant coefficients and the power of K is lower than that of D, 


It is assumed that oy =0 when w > Wep» i.e,, that Wep is the boundary of the linear filter passband, 


Function F (X) is assured to be odd, but this limitation has been introduced to make the problem more concrete 
and it can be easily omitted, 


For determining stability, let us introduce increments, Let 


X=X (t) + z. (2) 
Inserting (2) in (1), we find 
D (p) X (t) + D(p)z = — K (p) F(X @ +2). 


By developing the nonlinear function into a power series in x, neglecting the nonlinear terms for sufficiently 
small values of x and considering that X(t) satisfies Equation (1) we obtain 


























ry 
J D(p)a+ K (p) Ax (t)z =0, (3) 
__ [dF (x) 
Ax) = ("2 |, xa 
is the given periodic function of time (with period r ), 
Let us represent periodic function (3) by its development into a Fourier series ** 
Az (t) = Ax + Azmyz cos (Qt + a2) + Axmys COS (2Qt + tig) +... (4) 
For the investigated periodic condition X(t), determined by Equation (1), to be stable it is necessary and 
a sufficient for the solution of the linear differential equation (3) with periodic coefficients to be stable, 
at 
2m, Let us find a general solution of Equation (3) in the following form: 
lutions ) 
ing z(t) = Re Dd) AxePt+iceat+onn , (5) 
k=—a@ 
. where ) are proper numbers to be determined and X, are the required amplitudes, 
ne 
* A more general case Y = F(X, X', X*) can also be examined by means of methods similar to those described 
below, . 
ere 


** It should be noted that in (4) ayy = koyy; hence, it is possible to get rid of quantity x by suitably choosing 
the origin of the time basis, 


By inserting (4) and (5) into (3) and taking, for simplicity, my, = m4, after simple transformations, we shall 
obtain 


AoK (p) | >} XyePetsecnseps 4. ME Sr KyeOHstesvart oa) 4 


k= —oo k=—oa 





+a 3 XpeOttsle—naH+on)) i ol 3 X pelt tilee+natton]) 
k=——oo : k=— oo 
mer = AL+sIGe- NAL x 
+—z— DD Xecormemreald +. .] + Dp) Sy XneMtiooreon a 0. 6) 
k=—e k=—oo 


By changing in all the index sums (k + r) the quantity k into (k =r) and introducing notation X, = 
= X;, elk ,» we can rewrite (6) in the form 


k=—co k =»—oo 


+4 AGK (p) mye-3% J Kegs 4... + AgK (p) meio x 


k=—oo 


x > X pape HO) 4 © ALK (p) myei*e >» | OP aioe ee 1 (7) 
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If the complex amplitudes at similar frequencies are equated to zero, Equation (7) transforms into a 
system of an infinite number of equations, each with an infinite number of terms 


{[4pK (p) + D (p)) e380) Xy + [5 AymeimK (p) +900] Xp, + 
+ [5 amie HK (p)e+9001] Xess +... + [F ApmreierK (p) 0400] Xy_, + 
+ [+ Aym,e-3#r Kp) eO-inayt] heat ...08 Gaon... ~4;6, 1,...,00. () 

By equating to zero the determinant of this infinite system of equations we obtain a characteristic equation 
of an infinite power with respect to ). Its roots are the required values of . For the investigated periodic con- 
dition to be stable it is necessary and sufficient that all these roots have a negative real part, Obviously,in this 
form the conditions of stability cannot be applied in practice, But we have not yet used the filter conditions, 
which can now be applied, 


3. Filter Conditions and Limitation of the Power of the Characteristic Equation 








If it is taken into account that the cutoff frequency is finite and the complex amplitudes of components 
whose frequencies exceed the cutoff, are equal to zero, the infinite system of equations can be converted into 
a finite one, In fact all complex amplitudes, corresponding to the value k which satisfies the condition 


| Im Ay + 4Q| Dd erep, (9) 
will be identically equal to zero, 


Taking into consideration that the \; are unknown, it is possible to substitute condition (9), if the sign 
in front of kQ is appropriately chosen, by the approximate condition 


(10) 
k im® > Gp: 














Hence, the limiting expression of k = k};,,, which corresponds to the complex amplitude equal to zero, will 
be determined as 


@ 
him >a (11) 
If, for instance, Wep= 2Q,we shall obtain a finite system of five equations instead of an infinite system 
(8). 


If we denote the diagonal term coefficient (with a complex amplitude X) by &, ; and the mutual co ~ 
efficient at the crossing of the k line with the (k + r) column by ay yy 7, and take Wep = 2 Q, we shall obtain 
from (8) a finite system of equations which in a developed form will be written down as 


G_4.-sX— + 6,-4.2 + G_29Xo+ a_94Xy+ a_o2Xs = 0, 


@_4,-gX 4 + G4, X_y $+ ay Xo + a4 X1 + O49 Xy = 0, 
Mp, aX + Ao,yX_4 + AooX 9+ ay, X + Go 9X, = 0, (12) 
@s,-gX 9 + 64,-.X—3 + Gig Xe + Oy3Xs + Grghe = 0, 
p,—gX 9 + Gg,-2X_1 + Og,oXo + G9,1X 1+ a9X_ = 0. 


According to previous statements the coefficients in (12) will be, respectively, equal to 


Gre = [AgK (p) + D(p)} eo tO 


{ (13) 
Ox,k xr = ~~ Agmye**** K (p) et +00" 


dF (X) 


If the expression for the coefficient ax 


, expanded into a Fourier series, is limited to a 
X=X () 


finite number of terms equal to y, all the ax, terms with r> y will become equal to zero, Thus, with 
some of the terms becoming zero,Equation (12) will be simplified, 








In equating to zero the determinant of the final system of equations with complex amplitudes we shall 
obtain a polynomial in ), with complex coefficients, The roots of this polynomial are equal to the required 
values of )j;. 


For the periodic state under investigation to be stable, it is necessary and sufficient that all Re ,; > 0. 
When investigating whether all the ){ satisfy this condition,it is possible to use both algebraic and frequency 
criteria of stability, 


Taking into consideration that in this instance roots of a polynomial with complex coefficients are being 
investigated, it is possible to use Shur's criterion when applying algebraic criteria of stability, 


When using frequency methods it is possible to apply, for instance, Mikhailoy's criterion (more precisely » 
the principle of his reasoning, since in this case the complex roots of the characteristic equation are not conjugate), 


4. Comparison with the Stability Conditions of the Harmonic Balance Method 


In conclusion let us compare.the conditions of stability derived on the basis of assuming the existence of a 
filtering action of the linear part of the system, with the known conditions proposed by Gol'dfarb [2] and based 
on the harmonic balance method, 


As M, A, Aizerman and I, M, Smirnova [4] have shown, the conditions of stability used in the harmonic 
balance method can be obtained if it is assumed that the system is only slightly nonlinear and that it is therefore 
possible to find small deviations from the steady periodic state in the form of x(t) = c(t) siniw t+ y ()], where 
c(t) and y (t) are slowly changing quantities, and that the following averaging over one period is permissible: 


1 ((c + Ac) sin (cot — yo — Ay)] © (co + Ae) g (cy + Ac) sin (cat — to — Ay) + 
+ (€y + Ac) b (cy + Ac) cos (at — ty — Ay), 





where Ac and A y are small quantities, 


Let us find on the basis of the filter hypothesis the approximate stability conditions for the case when the 
cutoff frequency satisfies condition 0 < Wep = Q. 


In this case, according to (12) and (13), 
[ApK (p) + D(p)} eI", 5 AymyeiniK (p) e100 
+ Aymye-imK (p)e", [AK (p) + D(p)] 


= 0. (14) 


By solving the determinant we shall obtain a polynomial in). In evaluating the roots of this polynomial 
it becomes possible to answer the question of the periodic movement stability, If Re , < 0, the movement will 


. be stable, 





It is easily seen that these conditions of stability differ from those obtained by the autoresonance hypothesis, 


Let us now examine the relation between the stability conditions obtained by the filter hypothesis and those 
obtained by the autoresonance method, 


Let us make assumptions similar to those adopted in the autoresonance method, namely , let us assume 
that the system is only slightly nonlinear, i.e,, that it is possible to neglect terms containing m* and that owing 
to the slow change in the amplitude it is possible to neglect terms containing powers of ) greater than one. 
Under these conditions determinant (14) will assume the form 


T=] 


A, {a + | Siay_s7 gay] ht bm -+ 


+ | Son—s s (fQ)— h 4. (Gn—y\ + Gn) + (bm—yA + bn) “te 


8=—] 


— | Agmt{an + San_-r (ay), | =o. (as) 


T=] 


By solving the determinant, neglecting terms containing m? or powers of i greater than one, we shall see 
that the stability condition will then be reduced to the condition that the real part of the coefficient of ) must be 
positive, When expanded this condition will be written down in the form 


Re| Aoan | ets + bm + Sans? (jQ)'— + Somes vay-| + 


r=1 8=] 


+ Agbm | Avan + bmat y ayn—rT (7Q)*-1 + ) bin—s8 gay-||> 0. (16) 
T=] 


8-1 


Above condition ,.ven with'an assumed narrow passband of the linear part of the system [only the harmonic 
of frequency (Im )4— jQ) was taken into account], differs from the stability condition obtained by the auto- 
resonance method, 
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A METHOD OF ANALYZING AND CALCULATING TRANSIENT PROCESSES IN 
AUTOMATIC CONTROL OF GENERATOR EXCITATION BY MEANS OF A 
MAGNETIC AMPLIFIER 
V. R. Kulikov 
(Khar*kov) 

The author proposes a method of calculating the equivalent time constant of a 


unit consisting of the operating winding of a magnetic amplifier and a single-cycle 
bridge-type rectifier with an inductive load, simulating a generator excitation winding. 


1, Peculiarities of Transient and Stable-State Processes in 





Inductive Circuits With Rectifiers 





Transient processes in nonlinear circuits, containing semiconductor rectifiers are,as a rule, represented by 
periodic piece-wise-continuous functions, For the mathematical representation of such processes it is more 
common at present to use the method of successive solutions of differential equations, mostly transcendental or 
nonlinear, The solutions are obtained for a large number of discrete parts of the process with a subsequent 
matching of end and starting conditions of neighboring parts, 


The existing methods of investigating stable-state processes in circuits with semiconductor rectifiers, in- 
cluding the method used in[1], are based on the application of the harmonic balance method, Such treatment 
is possible because the stable-state processes in these circuits are represented by periodic piece- wise~continuous 
functions which can be expanded into a harmonic series, The existing methods of investigating stable-state 
processes in above circuits do not follow from the consideration of transient processes, These circumstances pre- 
vent the application of such methods for investigating the dynamic operation of semiconductor rectifiers, in 
particular single-phase bridge-type rectifiers in circuits for the automatic control of electric generator excitation, 


On the other hand, the existing methods of investigating transients in magnetic amplifiers, as applied for 
instance in the paper of L, A, Bessonov and M, A, Rozenblat, include the operation of magnetic amplifiers with 
a load on the ac side only, i,e., do not take into account the special transient conditions of operating jointly 
magnetic amplifiers and rectifiers, 


Therefore, the results of investigation of transients in magnetic amplifiers with an ac load cannot be 
directly applied to investigating the dynamics of the automatic control of electrical machines in which the 
rectifier, feeding the highly inductive excitation winding, is an indispensable component, 


All this shows how essential it is to develop methods of analyzing and calculating transient processes in 
inductive circuits with rectifiers, as applied to electrical machine automatic control systems which contain mag- 
netic amplifiers, 


Article [2] only deals with the operation of a magnetic amplifier without a feedback to an inductive load 
through a bridge-type rectifier, Such amplifiers are not normally used in automatic control systems, 


In addition, in the above article it is assumed that the core of the magnetic amplifier has an ideal rectangular 
magnetization characteristic, which does not correspond to the actual condition in the circuit, In electrical machine 
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automatic control systems the core of the high-powered output amplifier, whichis directly connected to the 
excitation winding of the machine, is made of electrical steel, 


The publication of certain results of paper [3] in this article serves as an attempt to fill the gap in the 
technical literature on this subject, 


The main peculiarity of inductive circuits containing rectifiers consists in sharp changes in their character 
istics, due to switching, i.e,, in rapid, practically instantaneous changes in circuit parameters, These changes 
are often accompanied by a change in the path taken by the current, For examining such processes it was found 
expedient to use as the basic quantity in calculations the definite integral of tension with respect to time p, = 


t 


= \u(z)dr, known as the tension pulse [4], The tension pulse has the dimensions of magnetic flux (in sec) 


ty 
and can be measured with a fluxmeter, 


An inductive linear or nonlinear circuit, with tension u(t) impressed on it, can be represented by an 
equation of electrical balance 


ad¥ : 
a a Ri => u(t). 


By integrating this equation term by term in the interval of u(t) operation y <t < ts, we shall obtain an 
equation for pulses 
ts 
AY + \ Ride = Bes 


ty 


In the pulse equation, as applied to the control of electrical machines excitation, A¥ represents the useful 
ty 
part of the pulse, which forms the flux, and {Ride the tension pulse losses, 
ty 


The transient process in inductive circuits with rectifiers is characterized by the effect of a large number 
of tension pulses which are relatively short compared with the circuit time constant, In order to represent this 
process it is necessary to find the relation between the pulse equation terms in neighboring time intervals. 


The basis of the investigation method consists in a precise solution of the transient process problem which 
arises when an ideal single-phase bridge-type rectifier circuit with two inductors lL, and L, is connected to an 
alternating tension, which does not contain a direct component and is,in the general casenonsinusoidal [3], In- 
ductor Ly is connected in the ac side and L, in the rectified tension side (Fig. 1). 


The circuit is made ideal in so far as it is assumed that the active resistance of inductors and the forward 
resistance of rectifiers is zero. and the reverse resistance of rectifiers is infinite, 


In [3] it was shown that the transient process (Fig, 2) consists of alternating excitation intervals, during 
which the rectified current rises, and switching intervals, during which the alternating tension pulse is spent in 
changing the magnitude and direction of the alternating current, but does not affect the magnitude of the rectified 
current, which remains unchanged jn that interval, The duration of the intervals is changed in the course of the 
process, The part of the pulse spent in the ac circuit, without reaching the rectified current circuit, was called, 
in [3] ,the tension pulse switching loss, 


This loss in the above circuit for (n+ 1) half-cycles can be determined from the formula 


Px (™%) = 2L,i (m), (1) 


where i(n,) is the instantaneous value of the rectified current at the beginning of the (n + 1)-th switching (all 
switchings begin at the moment the tension passes through zero), 
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In the above case Equation (1) represents the basic law of the transient process. 


The current at the beginning of the (n + 1)-th switching is equal to the current rectified during the switching 
i (n,) = iy = pgs + Sy, where iyy» and iyy* are the currents in circuits I and II (Figs. 1 and 2), 
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Fig. 1 
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During the (n + 1)-th half-cycle for a periodic 
function, or during the (n + 1} th interval between two 
adjacent tension zeros in a general case, inductor L, 
accumulates an interlinkage increment A¥2(n,), 
whose value can be determined from expression 














L , 
AY, (m) = rT [Pu — 2L,i (n,)]. (2) 


Fig, 2 


Expression (2) can be easily interpreted, The 
bracket includes the difference between the full tension 
pulse p, forthe (n +1)+h half-cycle and the switching loss of the pulse during the same period, The factor 
in front of the bracket determines the fraction of that difference which reaches the rectified current circuit. 


The discrete time argument of difference A ¥2 has an ordinal number a unit lower than that of the (n+ 1} th 
half-cycle under consideration, since A¥, is considered to be an increment of function ¥, at the end of the pre- 
ceding half-cycle, 


The value of the current i(n,) =i, does not change during switching. 





¥ 
By substituting in (2) atv) for i(n,),we obtain 
L L 
AY, (n:) = EET [P. —2 i Ys (n)| ° (3) 


Expression (3) represents a finite difference equation, By substituting the first difference by its definition 
A ¥2(n,) = ¥2(n, + 1) — ¥2(n,), we obtain a normal form of a difference equation 


lis — I. Ls 
VY, (nm of 1) aa yee Ee: T. (m2) = eT i Pu- (4) 
Since the phase and the moment ty of connecting the alternating tension to the circuit under consideration 
(Figs, 1 and 2) is arbitrary, the value of the initial pulse py is also arbitrary in the limits 0 <pyo < py. The 
initial pulse pyg determines the initial condition in the difference equation (4). 


| 
¥, (0) = Pe Sy Puo- 


The solution of Equation (4), as is known, is obtained similarly to that of a linear differential equation 
with constant coefficients and represents the sum of a general solution of a homogeneous equation 


Ig — I. 
YY. (mm 4+ 1 ) _— Ete i (72) = 0 











} th 


tion 
e 








and a particular solution of the heterogeneous equation (4). 
The general solution of the homogeneous equation is obtained in the form 


Wag (m) = cr". 
Subscript "f* in ¥g¢ stands for “free” and “c” for constrained, 


Quantity 4 represents the root of the characteristic equation, 


L.—ly 


oan T+; 


=0. 





Thus, 


einen 


The particular solution of the heterogeneous equation (4) is found in the form of a constant =A, By 
inserting this constant in (4) we obtain 


A= VY, = 7 ~ Pus 
The general solution of Equation (4) is 


Yalu) = gr at os (EE) 


After inserting the initial condition, we obtain 


Pu Puo ) 
a= — hls TFT) 
The solution of Equation (4) has the following form: 


Ys (m1) = 37 Pu La (sr — ren (ES2) od 





By dividing both sides of Equation (5) by Lz we obtain an expression for the instantaneous value of the 
rectified current ij; at the instants corresponding to the beginning of switching (zeros of tension) 





. Py Pu Puo I,—Iy "t (6 
re inn (m) = 37" — (37> — rt) EFL) : 
Denoting 
Ia— hn _ equiv (1) 
LF Li ; 


where A rt = . and Te aujy 18 @ constant, and w = Qnf is the angular repetition frequency of function u(t) 


zero values, 


Let us call quantity Tyqyjy the equivalent time constant of an ideal bridge-type rectifying circuit with 
two inductors, By inserting (7) into (5), we obtain 


—n Ot 


P “FT 
Ya(m) = yp Pu— La( gr — Zee )e 7 eau " 


The continuous exponential function 


\2, +l, 


ne 
falt) = 37> Pu— L(G Pos a Fe (9) 








with the values of its argument t = = By = ,3 =. eee, (k#1) 5 has the same values as function ¥,(n,) 
Ww w Ww 


with its argument n,= 0,1,2,...,k. ‘Hence, the continuous function of time f,(t), which has continuous de- 
rivatives, can be used as an interpolating function for ¥2(n,) in the intervals between the points common to both 
functions, 


R For dealing with actual circuits an equivalent circuit 
is proposed in [3], in which the single-phase bridge-type 
circuit with semiconductor rectifiers is replaced by a re- 

u presentation of a real rectifier in the form of three com- 


pO . ponents — a constant resistance, a source of constant back- 
u, emf and an ideal rectifier (Fig. 3)—and conditions of 
cham | switching are examined, On the basis of analyzing switching 

i conditions and observing by means of an oscillograph, it was 
| tan &=R shown in practice that the instant at which switching in- 
0 shy b fr tervals begin in transient and steady-state processes coincides 
Fig. 3 exactly with that at which the tension passes through zero, 

if in front of the rectifier (on the ac side) there is a choke 

coil with a saturating core, and there is a considerable 


linear inductance on the dc side, In this case the saturating choke simulates the magnetic amplifier and the 
considerable linear inductance the excitation winding of a machine. 














2. Combined Operation of a Magnetic Amplifier and Rectifier 





in an Electrical Machine Excitation Circuit 





The circuit under consideration consists of a self- magnetizing magnetic amplifier (MA) with an internal 
feedback, a single-phase bridge-type circuit with semiconductor rectifiers, an inductive load and a source of 
direct emf on the rectified current side, An equivalent circuit is shown in Fig. 4. The load parameters simulate 
the generator excitation winding and the resistance of the rectifiers and the operating winding MA, whereas the 
direct emf on the rectified current side represents a design value, determined by the value required to flatten out 
the excitation and rectifier characteristics obtained from curves ug, in Fig. 5 and e,, in Fig, 3 [3]. 


In determining the equivalent active resistance R 
in the rectifying circuit it is necessary to keep in mind 



































two essentially different cases. 
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Case 1, Resistance R, of the magnetic amplifier ac winding is small compared with resistance R of the 
excitation winding, This is the normal case. 


Moreover, for low-voltage excitation the number of rectifier elements per mesh is small (one or two) and 
their resistance R,,. = tan a (Fig, 3) is approximately equal to 0.5 ohm, which is small compared both with the 
active resistance R of the load and with the reactance of MA. In this case it becomes possible to neglect the 
active resistance in the ac circuit and take the equivalent active resistance R, as the sum of the excitation winding 
resistance R and the dc bridge resistance nR,o, Rz =R + MRro. 
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The equivalent emf E, is determined in this case by the expression Ez = E— 2neyg, where E is the actual 
direct load emf or an equivalent quantity, measured by a length of curve Ugo in Fig. 5, n is the number of series 
elements per mesh, ¢,, is the equivalent back-emf of a rectifier element (Fig. 3), 


Case 2, The active resistance of the MA ac windings and the rectifiers in the ac circuit is comparable 
in magnitude with the active resistance of the rectified current circuit, 


In this case in order to find an equivalent circuit the following simplifications can be made: 


a) the resistance of the MA winding and the rectifiers can be transferred from the ac circuit to the bridge 
diagonal 2-4 (Fig. 4); 


b) the rectifier design back-emf can be transferred to the bridge diagonal, 


In this case the equivalent resistance is determined by the expression Rg = Ry + (2n + No.) Ryo + R,, where 
Ngoc is the number of rectifying feedback elements in meshes 5-1 and 1-5 (Fig, 4), and the equivalent emf is 
determined by the expression E,= E — (2n + ng¢) ero. 


When an equivalent circuit was compiled in [3],use was also made of the operating peculiarities of recti- 
fiers and magnetic amplifiers in the generator exciter excitation circuit, peculiarities conditioned by the exist- 
ence of a considerable nonoperative exciter excitation current required for the generator nominal tension on 
open circuit, These peculiarities include: 1) a large switching interval; 2) a practically complete saturation 
of one of the magnetic amplifier cores during the excitation interval, 


The latter peculiarity permits one to consider the magnetic amplifier dynamic inductance to be constant 
in the excitation interval, i.e., to consider the MA magnetizing characteristic to be linear above the saturation 
bend (Fig. 6). 


Under conditions determined by these peculiarities, 
8 & the magnetizing effect of the MA control circuit, as can 
Ate A i be seen from the experimental curves, is practically in- 
7 ar dependent of conditions in the MA operating winding or 
d / the inductive load which is connected in the dc side and 
/ simulates the exciter excitation winding. It is, therefore, 
L / possible to represent the MA control circuit as a linear 
Pith H 0 % My directional first-order section (aperiodic section), The 
linearity control coefficients can be obtained experi- 
Fig. 6 mentally as in[5] by means of investigating a model 
MA of the type required,of an arbitrary power and size, 
but working under actual load conditions, For magnetic amplifiers of the same type but a different power and 
size paper [3] provides design formulas similar to those obtaining in the technical literature, 





























The basic aim of this paper is to investigate the transient and steady-state processes in the self- magnetizing 
MA inductive load, which is connected in the dc side of a single-phase bridge-type rectifier and which simulates 
the exciter excitation winding (Fig, 4). 


Transients in the MA are caused by changes in the control current with a constant supply tension, This 
condition is typical for dc generators with magnetically controlled power and tension, when the controlling de- 
vice is supplied from an independent ac source, It is also characteristic for the normal operation of synchronous 
generator magnetic regulators if the degree of static changes is small, when in the course of regulation supply~ 
tension variations can be neglected, being small in comparison with changes in control current which is supplied 
from a highly sensitive measuring element, 


The MA control circuit current deviations from the original value, which corresponds to the initial steady 
state, play in the above excitation circuit the role of a disturbing effect, In automatic regulation the excitation 
circuit is a sluggish link connected to the MA control circuit, In investigating the possibility of making the ex- 
citation circuit linear, both in the case of its analysis and for experimental work, the transient caused by a sudden 
change in the control current was taken as the basic form under consideration, When the circuit under investi- 
gation is fed by an alternating tension,such a transient process is superimposed on the preceding stable-state con- 
dition of the system (curve ig in Fig, 7). 
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The quantitative representation of transient and stable-state processes in the circuit under investigation 
(Fig. 4) was obtained by solving a linear difference equation with constant coefficients, The method of com- 
piling a difference equation was similar to the one used for transients inthe ideal circuit case. 
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Fig. 7 
We obtain the following difference heterogeneous linear equation with constant coefficients [3]: 


ale — tRe 


YP, (m: + 1) _— @ (Ly + Le) ¥.(m) = 
: (10) 


In Equation (10) ¥,(n, +1) and ¥, (n,) are instantaneous values of the magnetic interlinkage in inductor 
L, on the rectified current side, R, is the active resistance of the circuit, referred to the dc side, pa and i, are 
parameters of the MA magnetic system, related to coordinates B, and Ha of the core magnetizing curve (the 
beginning of the linear characteristic beyond the saturation bend) by the expression 


; l 
Pa=w,S; Ba, ia = 3 Ha, 
where wy is the number of turns in the MA operating winding, S,; is the cross section of the steel core, / ,, 
is the length of the core magnetic path , L, is the constant inductance of the MA operating winding with a 
saturated core, p...[n;] is the equivalent control pulse for the n-th half-cycle, determined by the expression [3] 


O.2mcSsrwiwe . 
Pec (™] = 7 — < i. (rm). (11) 
t 





In Expression (11) w, is the number of turns in the MA control winding, ic(n;]) is the control current at 


the beginning of the n-th half-cycle, p¢= yi is the permeability for the constant component of the 


lst 
0.4 © Sgr 
flux, y 4 is the coefficient determined from the statistical characteristics LA $9 = f(i,) which were obtained by 
means of a ballistic galvanometer or a fluxmeter for a given type magnetic amplifier with a magnetic circuit 
of length 1,, and cross section S,; [5], 


Current i, in the control circuit is determined by the linear differential equation with constant coefficients 


[5] 


. di _ 
E.= Riot we rr 10°. (12) 


When an initial rectified current is followed by a sudden reduction in the control current from a certain 
value to zero, the solution of the difference equation (10) assumes the following form: 


(13) 








wT ,. ; 

Pu + rae Fe + Lig Pec [0] al, — rk, ny 

Y2(™m) = Ly = +L, m ( ol; + ols ) , 
-— Ri +1 = eth st 











where E, is the equivalent emf on the rectified current side and is determined when the characteristic of the 
exciter and rectifier is made linear; it is included in the initial conditions of the process, 


When the supply tension is sinusoidal the expression for instantaneous values of the rectified current at the 
instants the tension passes through zero takes the following form: 


is (n,) = 








i 
0,9 ueget+ Es — 2f (p, — Lig) 2/Pe I) ( la— as 
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The continuous exponential function 


; ail 
0,9 uefft+ Ey — 2f (py — Trig) Pe I ; T equiv 
R; + 2/1; Ry + 2fly 


i, (t) = 
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(16) 


Lla— FF R, 


coincides with function ig(n,) at points t =0, ~ 2 = ooo, 0 i and can serve as an interpolating function 


in the intervals between these points, Expression (15) can be represented more conveniently in the form 


t 
, 0,9 ueget B2— 2f (Pp, —Lrig)+2fPe I FP, 191 ~ Tequ 
7 ef 2 Pa A P P i 
a(t) = Ra + 2fly Ry + 2fly lara "2. we 





The first term of Expression (17) represents the stable-state current before the application of the signal, 
The process which arises when the signal is applied (when the MA control circuit is opened) is represented by 
the second right-hand-side term, similar in form to the function which represents the transient process when a 
negative direct tension u is connected to a linear L, R circuit: 


t 
i=—+-(i—e 7) 


Let us examine more closely the second term of Expression (17), 


When an analysis is made of the process due to breaking the control circuit in which a demagnetizing 
current is flowing, it is necessary to take into account the negative sign of the “signal” pe_ (0), i.e., examine 
a positive function 


i=psat—e equiv) (18) 


Switching on a positive constant magnetization is equivalent to switching off a negative magnetization, 
Hence, Expression (18) also represents the process of applying a constant positive signal 2fPec (0) = const, 


Let us find a differential equation whose solution is function (18), Let us find it in the form 


2f Pe, [9] 
2+ 2jl,° 
Substituting current i, determined by Expression (18), and its derivative a in Equation (19) and equating 
the coefficients of the right and left sides of the equation we shall obtain B=1, A = BT gquiy = Tequiv- 


Equation (19) will then become 
di... 2fp.A0 
Cui +? = ee 


di , 








Taking into consideration (11) we can write 
di . 
Laminar +i = kyi,, (21) 
where 


0,4nu_ S. wyw f 
be c st c. 
i. (Hy + 2fL1) 





i,= 


The magnitude of current i,, as it has already been pointed out, is determined by Equation (12), 
Expressions (19) and (21) are equations of first-order sections (aperiodic), 


Thus, the circuit under investigation (Fig. 4) can be represented, when used as a part of an automatic 
control circuit, by a combination of two directional aperiodic sections. 


In the case when the automatic magnetic excitation regulator has no sluggish measuring element,* the 
regulator consists of the two above-mentioned inert sections, In automatic control the exciter is always connected 
to the very sluggish circuit of the generator excitation winding, which almost completely smoothes out the exciter 
tension pulsations, so that any further consideration of closed automatic excitation control circuits which in- 
clude a magnetic regulator can be carried outon the basis of known methods of the automatic control theory. 


Above analytical deductions, in particular Expression (15), were verified by a large number of experiments, 
The results of one of them are given on the oscillogram in Fig. 7, obtained with Tequiy = 0.0992 sec, Circles 
denote the calculated points of the theoretical curve, which in this case coincides with the experimental one. 
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* Such regulators include, for instance, automatic magnetic locomotive power regulators, which have a tachometer 
generator as a measuring element, 
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SYNCHRONOUS REACTIVE MOTOR SPEED REGULATION IN SYSTEMS 
OF PRECISE MAGNETIC RECORDING 


L. A. Pusset 


(Moscow) 


The author analyzes the stability of speed control by a synchronous reactive motor in 
a system of accurate magnetic recording, The stability conditions are obtained for the case 
when an electronic or electromechanical phase discriminator is used as the sensing element, 


The use of a nonperforated magnetic tape as a magnetic memory storage unit is complicated, in many 
cases, by the fact that the tape is not a perfectly elastic body and suffers changes in dimensions with the passage 
of time, The variable tension to which the tape is subjected during recording and subsequent playback pro- 
duces permanent deformations varying in magnitude along the length of the tape, This results in the recording 
frequency coinciding only occasionally with the frequency of the reproducible signal, and even then only for a 
very short length of time. Figure 1 shows an oscillogram of the beat note between a highly stable frequency, re- 
corded on tape, and the same signal during its playback from the tape, 











Fig. 1. 


It follows from Fig, 1 that the time scale experiences appreciable distortion during the reproduction of 
some information recorded on tape, This distortion in the time scale can be tolerated in sound reproduction, 
but for many forms of recording (phototelegraph and video signals, telemetered data, by means of frequency 
modulation) they are absolutely inadmissible, Consequently, in all known arrangements for accurate magnetic 
recording on a nonperforated tape, automatic control is applied to the mean tape speed, thereby achieving an 
almost complete compensation for all changes in the time scale which occur on the tape [1-5], 


The controlled element in most systems of accurate magnetic recording is a synchronous motor which 
keeps the tape in continuous motion through the tape winder[1-3, 5]. Figure 2 shows the block schematic of 
the automatic speed control system for the tape drive synchronous motor, A control signal, in the form of a 
constant frequency sinusoidal voltage is recorded, simultaneously with the information signal, on the tape, 
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During reproduction the phase of the reproducible signal is compared with the phase of a signal, of the same fre- 
quency, from the standard generator a, The output of the phase detector b is applied to a variable-frequency 
signal generator c, thereby changing its output frequency, The output of the signal generator is amplified and 
applied to the synchronous motor which drives the tape in the reproducer, Information given in foreign literature 
(1, 3, 5) indicates great accuracy in systems similar to the one described above. 
Thus, in the magnetic delay system installed in Singapore, during the XVI Olympic 
a Games in Melbourne, for the reception of phototelegraph signals, the error in 
maintaining the mean speed, during the time of transmission of one blank of 
phototelegraph information, did not exceed 107°, 




















b rt é In the present work, the stability of automatic speed control of a synchronous 
motor is analyzed and recommendations are made for the choice of the parameters 
of this system, Since magnetic recording systems incorporate almost exclusively 
synchronous reactive motors, all further discussions refer to such a type of motor, 
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Stability Conditions for the Control Process 








Fig. 2 The speed control of a synchronous motor must take place in such a way 


that the rotor is at no time out of synchronism, i.e., that the angle y, equal to 
twice the angle between the axes of the stator and rotor poles, never exceeds ymax » */2(6]. It follows from 
this that the rate of change of the frequency of the driving voltage can never exceed a known limit, determined 
by the rotor’s mechanical time constant, Leaving the question of the value of this limit open, let us suppose 
that changes in the time scale resulting from a nonuniform deformation of the tape take place so slowly that 
the motor does not go out of synchronism during the whole control process, 


Thea, after making all simplifying assumptions, which are usually made during the analysis of the transient 
processes in a synchronous reactive motor [6], the rotor’s equation of motion can be written in the form 


<e +r + apsine = o} sin (p+ 9). (1) 


Here ¢, is the stable value of the angle y, x is the varying component of the drive frequency, 


On the basis of what has been said, the quantity » is functionally related to the phase difference between 
two signals, one of which is the output voltage of the standard signal generator, and the other is taken from the 
control track of the magnetic tape, Let us suppose that the error signal is developed by a conventional phase 
discriminator, The rates of the stabilizing processes in such a discriminator are several orders higher than the 
rates of the transient processes in a synchronous motor. We will, therefore, assume that for every phase difference 
that appears at the input to the discriminator there appears instantaneously a corresponding voltage at its output. 
Corresponding instantaneously to this output voltage from the discriminator, there is some value of the varying- 


frequency component of the driving voltage. 


Assuming that the control system is linear and that the time scale variations occurring along the tape are 
sinusoidal, we can write 


d 
oY = acosMt + 9(p — q). (2) 


Here Q is the angular frequency of the sinusoidal variations in the time scale, g— gy» is the phase 
difference determined by the presence of transients in the motor, 


From (1) and (2), we obtain 
ats N%cos Qt + — + a a * rQ sin Qt + 


“hee (3) 
+ wf sin [; +e cost] = tsin [9 + 9]. 








Setting 7 = bu, we rewrite (3) in the form 


d? d2 
+ r st + wig sin [—$ iF + 90 — bu cos Mt] = 


= w2q Sin (po + ) — qby.Q? cos Qt — bugrQ sin Nt 


(4) 


or 


oy +r Fh —aigsin (9 + 4) + wig sin[ + St + 9] = 


(5) 
2 
= cos [so =a + 0] cos Qt — pqbQ? cos Qt — 





— pbgrQ sin Ot + p2F (3 hs 1) 


Here F (it ot Hu ) is a periodic function of time t with period = 


Taking » to be a small parameter, we replace Equation (5) by a series of first-order equations: 


dz, 


re | : 
ap = Ps — of sin 7 Iea+ Po + oq sin [2% + go] + 








ub 2 
+ a cos [- Z+9o| cos Qt — pgbQ? cos Qt — 


— pbgrQ sin Qt + pF (x, t, p), a = 2, (6) 


where X3=y. 
For 1 equal to zero, System (6) has the trivial solution x, = xg =x g=0, On the basis of this we shall seek 
a solution of Equations (6) in the form xj = yj. Expanding each of the functions sin + Xg+ »|. cos [pas ol 


and sin [Xg+ 9] into Maclaurin series about the point x; =0, we obtain 





d bq 
“7 wari fc (5 COS Poa + W29 COS Po¥s + —5 © cos gocos Qt— 
— gbQ? cos Qt — bgrQ sin Qt + pF, (ye, ys, t, w), (1) 
“ys 


d 
a YY Fr = Ye 


The following assertion is associated with systems of equations of the form of (7), If the equation 


Qy;—A Qy2 Gin 
a2) Gea—A don 

Ape =0, (8) 
ani ang Anan 


does not have any roots with their real parts equal to zero, then the system of the form (7) permits, for sufficiently 
small 4, one and only one almost periodic solution which reduces, for y =0, to a generating one (i,¢,, to a 
solution of the system in (7) with » =0)(['7), In Equation (8) aj, are the constant coefficients associated with 

the first-degree unknowns y, in the right side of equations of the form of (7), In our case, the time-dependent 
right sides of Equations (7) appear as periodic functions of t with period 2/Q, Thus, on the basis of the assertion 
presented above, for the case of the absence in Equation (8) of roots with zero real parts, the system of equations 
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in (7) has only one periodic solution with period 21/Q , which changes, for 1 =0, to a generating one, This is 
equivalent to the assertion regarding the existence of a periodic solution for System (6), with the solution be- 
coming identically zero for 4 =0, Let us find the stability conditions for this solution, First of all, let us note that 
the right sides as well as the periodic solution of System (7) are analytic functions of yu. 


Let us form the following equations for System (7): 


d 
st = Spite (i, k = 1, 2,3). (9) 


It is not difficult to see that pj, are periodic functions having period 2x/ Q and are analytic with re- 
spect tow, and that, for—> 0 any of the Pi, tends to the corresponding coefficient of the firs~degree unknown 
in the equations of System (7), 


Liapunov’s Theorem states that if all the characteristic exponents of unperturbed periodic motion have 
negative real parts, then this motion is asymptotically stable [7], On the other hand, the coefficients of the 
characteristic equation of a system such as (9) will be analytic functions of y and will become, for » =0, the 
coefficients corresponding to the characteristic equation of a system such as (9), in which the pj,'s are replaced 
by their limiting values for 0, i,e., aj,. In this way, the condition for the existence and for the asymptotic 
stability of the periodic solution of System (7), for sufficiently small p, is the existence of negative real parts 
in all the roots of the fundamental equation (8), 


In our case, a4; =— Ty 842 =— whe COS Yo, 43 = we GCOS Poo 2g = 1, agg = 0, agg = 0, agg = 0, agg =1, and 
agg = 0, 


The corresponding Hurwitz determinants are equal to: 


A, =7, A, = rw cos go + @2¢ 60S, As = — wg cos GAs. 
It follows from here that the stability condition is wp» q cos gp < 0,and since ¢ < - then q < 0, De- 
noting |q | =p,we have for the second condition the inequality r> p, 


In most foreign equipment, designed for accurate mag - 
netic recording on nonperforated tape, the phase detector is 
e usually a system such as the one shown schematically in Fig. 3 
[1- 3], Here a and b are two synchronous motors with their 
rotors rigidly coupled to two shafts of the differential c. 
The differential's third shaft d, the angle of rotation of which 
is equal to the difference between the angles of rotation of the 
first two shafts, is rigidly coupled to the potentiometer e 
or to a system of selsyns which operate as transformers, 




















The stator winding of one motor is fed by a voltage 
taken from the output terminals of the standard oscillator, 
while that of the second motor is fed by the voltage from the 
control track of the magnetic tape, In the presence of a phase 
difference between these two voltages, the angles of rotation 
of the rotors are unequal and hence the shaft d turns through 
a certain angle and thereby changes the output voltage from 
the potentiometer e, This leads to a change in the driving 
voltage’s frequency and the speed of rotation of the motor 
driving the tape, The equation of motion of this motor's rotor 
can, as before, be written in the form 























SP 4 p SP 4 of sing = a sin (9, + 9). (10) 


The equation of motion for the rotor of the motor which is fed by the voltage derived from the control 
track will have the form 





als 
soli 





d , / ai 
ft 4, % +. wi sin y = w3sin (x, + 8), (11) 


where y and 8 are quantities of the same type as g and ». 


Since the motor driven by the standard voltage is not included in the control feedback loop then, assuming 
that its rotor is free from the influence of any disturbances, we can assume the rotor‘s rotation to be strictly uni- 
form, In other words, the rotating rotor of this motor appears as an original metering system, with respect to 
which the rotor connected to the other shaft of the differential rotates, The constant component of the relative 
angle of rotation of both rotors corresponds to some nominal value of the frequency from the driving signal 
generator which, in turn, is equal to the nominal frequency of the reproducible signal, Each change in the re- 
lative angle of rotation occurs with a corresponding additional turning of the rotor of the motor, the driving 
voltage for which is taken from the control track, Hence, we can write 7 =n(y — yo ), where 2 is the 
varying -frequency component of the voltage driving the tape- transport motor, 


On the other hand, 8 =a cos Qt+b(y— qm). Let us set f(y — yo) dt = €. Equations (10) and (11) can then 
be rewritten in the form 


3 ‘ . 
=a +r ri + be% sin [4 + %— - cost] — bo? sin[po+ne]= 


= —a 2*cos Ot — a Q sin Ni, 


(12) 


d%¢ ,d% ma. fd to ot 
ae +? 7a + o@sin| 5 + x0] = o7sin [xo + B]. 





Let us set a= yc and assume p to be a small parameter, Equation (12) can then be replaced by the system 
equations 


dz, 


= = — rt, — o? sin [% + xo] + ao sin [yo + 2], 
d d 
= = 7,5, = = %, 
— = — rz, — bw? sin (- + % — + cost2e] + be? sin [Go + nx] -— (13) 
— pc? cos QD t — wcQ sin Q, 
dz, 


a ~™ x, = B, wy = &. 


Let us rewrite Equation (13) in the form 





ce = — rs, — of sin [%,+ Xo] + 2 8in [Yo + %), 
d d. ‘d . 
ttn, ey, Ht = —ro,— bop sin [24 9 ] + 
pcus 
aa 9 ® cos| = + 9] sin Qt + beg sin [9 + nx] — (14) 
— peO?cos Dt — peQ sin Ot + y*F (2, t, p), 
dz, 
“ge Zar 


where F (x, t, 4) is a periodic function of t, with period = , and is analytic with respect to the parameter 1, 
For 4» =0, System (14) has the trivial solution x; =0, All the above discussions, concerning Systein (6) are 


also applicable to System (14), The condition for the existence and for the asymptotic stability of the per'odic 
solution for System (14) will be the absence of roots, with positive or zero real parts, in the fundamental equation 














Qy,— ayy 43 ang a5 
aq) Agg—A dos Ao, a2, 
a3) Az2 Agg—A Agy as, = 0. 
a4) AQ42 Q43 Ayye—h Ay 
G51 Aso as3 54 As, — A 
Here dix = (ast) , and X; are the right halves of the system of equations (14). 
[X,— 2X, X— xX, X,—0, p —O0] 


For the case under consideration: 4), = Qy3 = Oy, = 2;,=0, ay. =1, A, = 
= — G2 COS Py, dag = —T, Aes = bnw? COS Py, Ao = 0, ao, = 0, a1 =sa= G9 = 


= ds, = 0, ay, = 1, Ay = Ag = Qg3 = A= 0,7 a,, = 1, Oy) = Oy" COS Xo) Agq = 0» 
Ogg = 0, Ogg = — 2 COS Xo, O55 = — 7". 


The fundamental equation will be written in the form 


MB + Mr + 7’) +28 [rr’ + @? cos yo + @2 cos 9] + 
+? res? COS Xo + re COS 9] + hare? COS Pp COS ¥y — as) 
— w2w?bn COS 9p COS Xo = 0. 


From Hurwitz's conditions it follows, firstly, that — wi wi* bn cos gp cos xo > Oand, since @ < - : 
xo < : — , it must be that bn < 0, 


The second Hurwitz condition A; > 0 has the form 
(r +r’) [rr’ + w2 cos xo + «2 cos 9] > rw? cos Xo + 7’e?2 COS gp. (16) 


This condition is satisfied identically for r# 0, r° # 0. 
Hurwitz’s third condition 2, > 0 leads to the inequality 


[rw? cos x9 + r’w? cos 9] (7 + 7’) [77’ + @2 cos x + w2 C08 g] + 
+ (7 4- 7’) «20? cos 9 COS Xo > [72 COS Yo + 7’? cos Py}* + (17) 
+ (7 + 7’) wo? | bn | cos 9, COS Xo. 


Finally, Hurwitz’s fourth condition is written in the form 


[(r +r’) rr’ + rez cos 9, + r’a? cos xo] [702 cos xq + 7’w? COS Gy + 
+ bn(rr’ + ©)? COS Xo -+ @2 COS Py)] > (18) 
> wiw? cos Py cos x2 [(r + r’) + bn]?. 
Inequalities (17) and (18) impose definite conditions on the values of the coefficients b and n. To de- 
termine these conditions let us carry out a certain transformation of the inequalities (17) and (18), As a result 


of this transformation both inequalities, which must be satisfied by the absolute value of the product | bn| , re- 
duce to the inequalities 


A 
D|bn |? + (C — 2DB)|bn| + DB? -—-A<0. (19) 
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Here , ’ oe , ee 
A=[(r+r)rr’'+ rw’ COS Py + rw? cos ¥9)[rr’ +r ©? COS Yo-+ W2 COS Po], 


B=(r-+r’), 
C=[(r +7’) rr’ + raf cos 9 + 7’? cos Xo] [77’ 4+- @? COS Xo + w2 C08 gp), 
D = ww? COS Py COS Xo. 
Since forx-> — oo andx—> o the three-term equation Dx*+(C— 2DB)x +(DB*— A)—> +o, then 
for the second of the inequalities in (19) to be satisfied it is necessary that the following relation be satisfied: 


Ly <. | bn| <2, (20) 


where Xq and x, are the roots of the equation Dx* + (C — 2 DB) x + (DB*— A) =0, 


The use of the electromechanical phase-sensitive device discussed, instead of a simple electrical phase- 
discriminator circuit, has the following advantages: 


1) although the introduction of extra degrees of freedom complicates the control system, nevertheless, 
this makes it more flexible, i,e., it increases the freedom in the choice of values of the separate parameters for 
the provision of stability of the control process; 


2) any of the purely electrical phase-discriminator circuits impose very high requirements on the stability 
of the output voltage,in the absence of an error, since a given relative change in the magnitude of this voltage will 
bring about a corresponding relative change in the output frequency of the signal generator, 


In the case of the phase~sensitive electromechanical device, any change in the voltage driving each of the 
motors will bring about the appearance of a disturbing moment acting on the corresponding rotor, The hunting of 
this rotor will lead to a relative change in frequency produced by the signal generator. However, this change will 
be smaller relative to the voltage change, 
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INPUT CIRCUITS OF CONTACT-MODULATED AMPLIFIERS 


D, E, Polonnikoyv 


(Moscow) 


The main requirements for the input circuits of electronic amplifiers in automatic 
potentiomesers are set forth, A table of diagrams of the input circuits, containing contact 
modulators, is given, as well as the formulas for determining the parameters of these cir- 
cuits, The features of each circuit are set forth and recommendations are made as to their 
application in amplifiers of automatic potentiometers, 


Wide application has been found in automatic control equipment of amplifiers in which an incoming dc 
signal is converted to an alternating one at the input, The dc-to-ac conversion results in a considerable reduction 
of the zero-level drift, and thus, overcomes the basic shortcoming of dc amplifiers. Amplifiers with dc-to-ac 
conversion of the input signal find application in various systems of automatic control, tracking systems, auto- 
matic potentiometers, modulators, and other equipment where low-level dc voltages are required (of the order 
of a millivolt or less), In most cases the input circuits of such amplifiers contain contact modulators possessing 
the smallest drifts among those presently known [1-4]. 


Many circuits, each one unique in a different way, exist for connecting the contact modulator, The proper 
choice of the input circuit and its components is governed, to a great extent, by the amplifier’s characteristics and 
the system as a whole, 


In the literature of this country there have been no papers published on the comparison analysis of the in- 
dicated systems, In this article the results of experiments carried out on the more common input circuits, con 
taining contact modulators, are given, and recommendations are made regarding their application in amplifiers 
of automatic potentiometers, 


A series of special requirements is imposed on the input circuits of amplifiers used in automatic potentiom- 
eters. The more essential ones for the analysis and evaluation of the circuits being studied are listed below. 


The input circuit must convert the incoming dc signal to an ac signal at the network frequency, At the 
same time it is expedient to have the largest possible conversion factor (ratio of the effective magnitude of the 
fundamental at the output to the direct signal at the input), Insofar as the input signal contains, in most cases, 
a significantly large ac component of the network frequency, it is desirable for the input circuit to discriminate 
against, or significantly attenuate, the indicated alternating component, In order to prevent overloading of the 
source and yet obtain maximum sensitivity, the input resistance of the circuit must be at least several times 
larger than the output impedance of the measuring circuit, It is imperative that the disturbances, introduced by 
the input circuit (among them the zero drift), possess minimum values and do not extend into the regions of the 
equipment’s sensitivity, The input circuit must be as insensitive as possible to the action of external electric 
and magnetic fields, If the amplifier is designed for operation in fast-response systems, then the dynamic pro- 
perties of the input circuit must be such that the period of the transient process* is not greater than 3-4 periods 


* By period of the transient process is understood the time for the output to reach 98% of its steady-state value, 
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of the conversion frequency, i.e., 0,06-0,08 sec [1]. In those cases where, for some reason or other, the input to 
the amplifier cannot be grounded, the input unit must provide for the isolation of the input circuit, 


There are many varieties of input circuits containing contact modulators, Table 1 lists the more common 
circuits and the corresponding equations for their basic parameters: the conversion factor Kp» input resistance 
Rin» Mean input resistance Rin means 4nd the selectivity factor a, 


Equations for a are given for the case R; =0, In the table, g is the phase angle of the input signal re- 
lative to the phase shift of the contact modulator, L, is the inductance of the transformer secondary winding, 
R is the total transformer loss resistance referred to the secondary winding, Q is the quality factor of the trans- 
former, R; is the internal resistance of the signal source, 


The derivation of the equations for most of the parameters of the circuits in Table 1 is very lengthy (see 
[1J) and is not given here, Let us, therefore, limit ourselves to a few explanations, 


The pare meters of transformerless circuits, (1 to 8 in the table), were derived on the assumption that 
2fyCoR, > 1(f 9 is the modulation frequency) and the contact-travel time is zero, Moreover, the intercontact 
capacitances as well as the capacitances of the contacts to ground were ignored, this being permissible for input 
resistances up to 10 meg when fy = 50 cycles, 


Under the above assumptions one can assume that the output signal consists of rectangular pulses, To de- 
termine K,,, the pulse amplitudes were obtained, whereupon the output signal was expressed in terms of a Fourier 
series, and the fundamental component, together with K p> Was determined, 


Rin was derived on the basis of the condition that for Rj =Rjp the transfer coefficient Kp, is reduced by a 
factor of two, Knowing Rjp permits one to find Kp for a given source internal resistance Rj. 


The mean value of the input resistance Rin mean Was determined from the relationship between the input 
voltage and the input current, averaged over one period. Rinmean characterizes the loading of the source by 
the input circuit, The selectivity coefficient « was determined from the relationship between the conversion 
coefficients for dc signals (useful signal) and ac signals (conversion frequency), 


It is apparent from the table that circuits 6 and 8 possess the largest Kp, Rin, and Rinmean among the 
transformerless networks, Circuits 3, 4, and 7 possess slightly poorer parameters; however, their advantage rests 
in the fact that one modulator contact is free and can, in some cases, be utilized for other purposes, such as a 
demodulator, Circuits 1, 2, and 5, have one serious shortcoming: the absence of a decoupling capacitor results 
in the added dependence of the zero drift on the grid current of the input tube, This makes the application of 
these circuits in sensitive amplifiers highly undesirable, The transformerless circuits 1-6 possess one common 
disadvantage: their selectivity coefficient is very close to unity, As a result, considerable error will be in- 
troduced by the parasitic conversion-frequency input signal, If these circuits are utilized, the alternating com- 
ponent of the input signal must be suppressed to a level below the sensitivity threshold, The reduction of the 
alternating component necessitates the use of filters which have a very appreciable detrimental effect on the 
dynamic properties of the system, Therefore, in the presence of an appreciable alternating component in the 
input, it is expedient to use circuits 7 and 8 which are free of the above-mentioned shortcoming [6]. Because 
the parasitic alternating signal is delivered both to the grid and the cathode of the input tube, the fundamental 
component in the output is attenuated, If the circuit parameters are chosen so that the following conditions hold: 


for circuit 7 
| 4 1 
ees (ae + a) (1) 
Ri e+ Ra k 
and for circuit 8 


Ra=R (1+), (2) 


where : is the static amplification factor of the input tube, then the fundamental component is completely 
compensated and the parasitic signal is transformed to a spectrum of even harmonics, as can be seen from the 
oscillogram in Fig, 1, 
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TABLE 1 (Continued) 

























































































Circuit | i 
No. Circuit Basic Parameters 
K V2 2K, (Ri + Ba) 
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The conversion factor for the useful signal in circuits 7 and 8 is the same as in circuit 6, The input re~ 


sistance of circuits 8 and 6 is also equal, The input resistance of circuit 7 is considerably lower and thereby 
limits its applications, 


We can conclude from the comparison of the transformerless circuits that circuits 6 and 8 are the most 
suitable for application in autocompensator amplifiers. Let us establish other properties of these circuits, A 
more accurate value of Kp for circuit 6 can be determined from the formula 





K V2 4 “Ry+2R, 
po EE ~ FeR, RFR, (3) 


The phase shift of the fundamental frequency relative to the phase angle of the contact switching is de- 
termined from the expression 


z 


lh ot FR, ~ sg) " 





tang = 


Since the condition of constant phase angle [1] is met in autocompensator amplifiers, the dynamic pro- 
perties of circuit 6, relative to the envelope, can be characterized by the transfer function which, in operator 
form (for 2 fgCcR, > 1) has the form: 
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V2 R, Tr a—i (5) 
I (Pp) = = RR, (! + wk)! +a 
where p is the differentiation operator and a = jee ; 
I 
0.02 s 
I 


















































Fig. 1, Shape of the output voltage of circuit 8 (curve Il) for a sinusoidal 
input (curve D. 


The transient response of the circuit to a unit step function is shown in the oscillograms in Fig, 2, Ob- 
viously, to decrease the transient time it is necessary to decrease T, =C,R, and a. With proper choice of para- 
meters the transient time can be reduced to a value 2 to 3 times the carrier period, as is apparent from the 


oscillogram in Fig. 2,c. 


It should be noted that an accurate value for the input resistance of circuit 6, for R, > 10” ohms, cannot 
be obtained from the formula given in the table because, under these conditions, the conductivity associated 


with the overcharging of the contact capacitances and the capacitance of the reed to ground [1] has an appreciable 
value, , 


The dynamic properties of circuit 8 are somewhat poorer than those of circuit 6, especially for large values 
of Rj (Fig. 3, a), which is involved in the overcharging of capacitor Cd. Therefore, it is expedient to use circuit 
6 when a high input resistance is necessary. 


To reduce the transient time, the following process for calculating circuit 8 is recommended : 
1, Choose the value of R;, for normal self- bias conditions. 


2, Knowing the required input resistance Rin, choose R, =a 


3, Knowing R;, and R,, it is necessary to choose a time constant T; = Cg [asm (+ AY), ene Cq 
is determined (T, determines the transient time). 


4, From the condition of equal phase shifts in the grid and cathode circuits,C,. is determined: 
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Fig. 2, The transient response of circuit 6 for three sets of values of its parameters: 
a) R=0, T,=0.2secz b) Ri =1.5x 10° ohms, T, =0,2 secs c) Ri =1.5 x 10° 
ohms, Te = 0,02 sec, 
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5, From the condition for the compensation of the fundamental component,Ry is determined: 
2 


Figure 3, b shows the transient for parameter values chosen in accordance with the above suggested method, 
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Fig. 3. Transient response of circuit 8 for two sets of para- 
meter values: a) Rj =50 kohm, Ry =51 kohm, R;, =50 kohm, 
R, = 240 kohm, C, = 0.5 uf, Cy = 10 wis b) Ry =50 kohm, 

Rg =10 kohm, R,, = 10 kohm, R, = 240 kohm, C, = 0,06 uf, 
and C;, =0.7 yf, 


Interference, caused by internal electric and magnetic fields, has a very negligible effect on transformerless 
input circuits, An electrostatic shield usually proves sufficient in a properly set-up circuit, It should be noted 
that the effect of electrostatic fields is reduced by using a smaller value of Rip. Circuits with a grounded moving 
contact have, in this respect, a significant advantage which, however, is only effective in the absence of screening 
of the excitation circuit in the contact modulator, The zero drift brought about by the variation of the thermo- 
electric potentials which exist at the junctions of dissimilar metals is almost completely independent of the ci~ 
cuit chosen, and is mainly determined by the selected materials and the differences in temperature between 
various points in the input circuit. | 


Let us consider the converter with transformer input (circuit 9). This circuit has found wide application, 
especially in autocompensator amplifiers, The main advantages of this circuit are: 1) isolation of the input is 
simple (1); 2) a large value of the ac selectivity coefficient can be easily achieved; 3) a large value of Kp 
is easily attainable, 








The last advantage is especially important because an increase of the signal level, prior to the first am- 
plifier stage, improves the signal-to-noise ratio at the input and decreases the errors for a given sensitivity or, 
for the same error, raises the amplifier's sensitivity [2]. 


Under ideal symmetrical operation of the contact modulator, the transformer input circuit discriminates 
against conversion frequency parasitics, In practice, a finite difference always exists between the closing time 
of the first and second contact, i.e., 


p= +0, (7 


where ty is the closing time of the first contact, tg is the closing time of the second contact, T is the operating 
period of the contact modulator, 


Determining the fundamental component of the output signal, for a given 6, and assuming that its value 
is below the sensitivity threshold, we obtain the condition connecting the permissible ratio of the alternating 
component (network frequency) at the input U, to the sensitivity threshold U,, and the assymmetry in the contact 





modulators operation a. ° 
heey y 
_— Oa\ti—ts/]° (8) 
. iy ty 
In the derivation of Equation (8) os was assumed to be small compared to unity. If it is assumed 


that the maximum permissible assymmetry is 10%, then it follows from Equation (8) that the ratio of — can 
t 


be no larger than 75, It should be noted that the parasitic alternating signal, besides the fundamental, results 
in an appreciable value of the second harmonic which may lead to saturation of the amplifier and a reduction 
in its sensitivity. To avoid this it is necessary to reduce the amplifier's bandwidth, A very simple means of 
suppressing the higher harmonics is the use of a tuned transformer and the introduction of an extra resistor Rj, 
in series with the input, which reduces the shunting effect of the source's 

Rin bohm input resistance and thus raises the transformer Q. It has been established 
experimentally that, for Rj = 8 kohm, the second harmonic is attenuated 
by a factor of 10; at the same time Kp is reduced by only 30%, The 
transformer used in the experiments had the following data associated 
with it: dimensions— 24x 32x 44; w,=800 x 2, PEL~ 0.1; w= 
= 20,000, PEL— 0,05; I1-shaped core, 0.8 cm? cross section, permalloy 
N79M5, The indicated method of attenuating the second harmonic 
is, to us, the most rational one, In the table, the expression for Rj, 
~Q2 -ai 0 a Qa? 46 in circuit 9 is given for the case of ideal operation of the contact- 

modulator where the contact travel time is equal to zero, 








ss Figure 4 shows the dependence of Rip on the relative, contact- 
travel time 7 . & < 0 is understood to be the contact flashover time. 
The dynamic properties of the transformer circuit can be approximated quite closely by the transfer 
function 
f(p) = 2V2 ARin 4 
P me Rint Rig4 2’ (9) 
' Tho 


where Q is the transformer quality factor, fo is the conversion frequency, n is the transformation factor, p is 
the differentiating operator. 


The following shortcomings of circuit 9 should be noted: 


1) appreciable sensitivity to internal, alternating, magnetic fields which necessitates the use of specially 
constructed transformers and thorough screening [1]; 




















2) difficulty of achieving a high input impedance; 
3) relatively high cost, 


In concluding, the following recommendations on the use of the described circuits with contact converters 
can be made: 


1, In those cases where isolation of the inputs is not required and the permissible sensitivity threshold is 
10 microvolts, the most rational choice is the transformerless circuit 8, 


2, If an isolated input is required and a sensitivity threshold of less than 10 microvolts is necessary, the 
advantages of the transformer circuit are obvious (circuit 9), 


3, The application of circuit 6 is expedient in those cases where a very large input impedance with high 
dynamic characteristics is required, as in the measurement of small currents, 


4, The application of the other tabulated circuits in autocompensator amplifiers is, in general, inefficient. 
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A PUNCHED CARD METHOD FOR SYNTHESIZING SEQUENTIAL RELAY SYSTEMS 


V. I, Shestakov 
(Moscow) 
A method is presented for implementing a vector-algebraic method of synthesizing 
sequential systems of two-position relays by means of special punched cards, This method 


is applicable for the synthesis of both autonomous and nonautonomous relay systems, 


1. Introduction 








For facilitating and speeding the process of synthesizing relay circuits, both relay contact and vacuum 
tube circuits, a number of punched card methods have been proposed, but in all the works [1-6] presenting these 
methods what was discussed was only the minimization of the expression of the Boolean functions describing the 
structure of the relay circuit to be synthesized, These methods are quite effective since they allow a significant 
portion of the work and time expended on simplifying the circuit to be synthesized to be saved, They are parti- 
cularly effective in the synthesis of one-cycle or, as they are still called, combinatory circuits, since the simpli- 
fication process is the most difficult step in the synthesis of such circuits, 


These card methods may also be used for the synthesis of sequential relay circuits, However, simpliiication 
is only the final step in the synthesis of sequential relay circuits, This step is overshadowed by the determination 
of the conditions for operation and release of each of the relays in the circuit to be synthesized, The determina- 
tion of these conditions is carried out from data, given in some manner (for example, by means of a switching 
table), on the sequence of functioning of all the relays in the desired system for given sequences of exteinal stimuli 
(signals), This step, specific to sequential relay systems, is completely independent of the concrete elements 
(electromechanical or ferrite relays, vacuum tube or crystal diodes, etc.) from which the circuit to be synthesized 
will be constructed, and of which the equivalent relay circuits of the system will be synthesized, Thus, we re- 
place the term “relay circuit” by the more general term “relay system," introduced by the author in [7]. 


A vector-algebraic method [8-10] is proposed for the synthesis of sequential relay systems, It is always 
possible, by using this method, to obtain Boolean functions formulating the necessary and sufficient conditions 
for the functioning of each relay in the system to be synthesized, given the sequence of functioning of all relays 
in the system to be synthesized and the sequence of external stimuli to this system, The computation of these 
functions is very simple and does not require much time if the number of independent variables (external signals) 
and the number of relays in the system to be synthesized are small, However, with an increase in the number of 
independent variables and relays in the system to be synthesized, there is a very rapid increase in the length of 
the sequences of external signals and of the sequences of relay states which are to be synthesized, Because of 
this, the computations necessary for the synthesis of sequential relay systems with an increase of independent signals 
and relays rapidly become ever more extensive and burdensome, with the result that the probability of error in 
the synthesis also increases, The requirement arises of facilitating and speeding, in some measure, the necessary 
computations by means of mechanization of at least the simplest of them, 


This paper presents one of the simplest methods for mechanizing the vector- algebraic method of synthesizing 
sequential relay systems, This method is based on the use of special cards, which gives it its name, a punched 
card method of synthesizing sequential relay systems. For the exposition of this method we shall assume that the 
reader is familiar with the vector-algebraic methods of synthesizing sequential relay systems, Hence, we shall 
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give here only the minimum details of this method as required for the understanding of the method, presented 
here, of using special cards for the mechanization of the vector-algebraic method ot synthesis. 


2. Short Exposition of the Vector-Algebraic Method of Synthesizing 





Sequential Relay Systems 





At any moment of time t, the state of the n-relay system Y,,..., Yn can be unambiguously defined 
by giving two vectors: the vector of the independent parameters x (t) =[x,(t),..., (0), t,....Xm@), 
controlling the given relay system, and the vector of the dependent parameters y (t) =[y,(t),...,y,(0,..-> 
Yn (0)], controlled by means of the relays of the given system, When the relays of the given system are electro- 
mechanical, the variables %,..., Xi»... » Xp, afe contacts closed by relays or buttons Xy,...,Xj,..+>» 
Xm» external to the given relay system, and the variables yy,...,Yke-++»Yn are contacts closed by operation 
of the corresponding relays Yy,...,Y,....»Y, of the given system, 


If all the relays in the given system have identical lags in operating and releasing, and if all the controlling 
and controlled parameters are varied in virtual synchronism, then the process of changing state of this system 
can be described by the following vector equation: 


y(t-+2) =F(x(t), ¥(o), = 
where f is a single-valued vector function, 


Determining the process y(t) by giving the function f and the process x (t) is called analysis, while com- 
puting the function f(x,y ) by giving the processes x(t) and y(t) is called the synthesis of the system described 
by Equation (1), 


Under the assumptions made above with respect to relay functioning in the given system, it is possible to 
replace the processes x(t) and y(t) by the sequences x(j) and y(j), where j=0, 1, 2,.... Then Equation (1) 
is replaced by the following equation: 

y(7 +1) =f(x(7), y)), (2) 


where the natural* variable j (j=0, 1, 2,...) will be called the cycle number of the sequential (multicycle) 
relay system described by Equation (2). 


A relay system is termed autonomous if each controlled parameter in it is controlled only by relays lying 
in the same system. Such a system is described by an equation of the form 


y (7 + 1) =f(y()). (3) 
A set of several sequences y @ is called complete if among the terms of this sequence all 2" values 
Yg (8 =0,1,..., 2°— 1) of the variable vector y occur. 


It follows from the single-valuedness of the function f that for the synthesis of eas oe — only 
such complete sets of sequences y(j) are suitable for which each pair of sequences. y * “(j) and y ($j) satisfies 
the condition 


Af y(t) — yO (7) then 4 (71) + 1) = y (J $4), (A) 
where r ands independently run through the values 1, 2,..., NN is the number of all the sequences y (j) 
in the given complete set) 


In contradistinction to autonomous systems there are nonautonomous relay systems. In particular, if a 
relay system does not contain an autonomous subsystem it is called completely nonautonomous, or a single-cycle 


system, 


An algorithm for synthesizing autonomous systems of two-position relays from a given complete set of 
sequences y (j), satisfying Condition (A), is given, as was shown in [8], by the formula 


* By a natural variable we mean a variable which runs through the series of natural numbers: 0, 1, 2, 3,... 





t+ cepa 6 





N Ny-1 
CY=D DY Pony U +4), (4) 
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ral jmo 


where £ is the sign for the Boolean sum, N is the total number of sequences in the given complete set, N, is the 


number of states in the r-th sequence and pg(y) are constituent units computed, as shown in [9], from the vector 
components by the formula 


Ps(y) = i (Yen @ yx @ 1), (5) 
-1 


where y,, and Yg,k re, respectively, the components of the variable vector y and its values. Yg° and @ is 
the sign for mod 2 addition, 


The algorithm for synthesizing nonautonomeus systems of two- position relays, i.e., relay systems described 
by Equation (2), is given, as shown in [10], by the formula 


N N;—1 


FaN= 2D > Pn Pw (v¥G + 1), (6) 


T=1 j-0 
where the p, (x) are constituent units, computed from the components of the variable vector x and its values, 
X,, by a formula analogous to Formula (5), Formula (6) is a generalization of Formula (4), 


Sequences x (j) and y(j) for which the function f(x, y) may be computed from Formula (6) must satisfy 
the condition 


if x(r) i) = x(8) (7)and y” Uj) = y” (7), (B) 
then y” (7 +1) = y” (j@ + 1), 


where r ands independently run through the values 0, 1, 2,..., N (N is the number of pairs of sequences 
x(j) and y(j) in the given complete set), Condition (B) is a generalization of Condition (A), 


If the given sequence does not initially satisfy Conditions (A) and (8) then, in many cases, it may be 
extended so that these conditions are satisfied, this extension being carried out perhaps by the method described 
in [9]. 


3. The Use of Special Cards for the Synthesis of Autonomous Relay Systems 





The values y, of the variable vector y , characterizing the state of the relay system Yy,..., Yk, 
. » +» Yp, May be numbered in different ways. As in earlier works [10, 8], we will number them by the 
numbers 8, obtained from the components yg , = by of the vector y,, by the rule expressed by the formula 


B= y b,2*—. (1) 


k=1 


Conversely, by using this rule it is easy, being given the number 6 of an arbitrary value y, of the 
variable vector y, to write all the components of this value, For this it suffices to write the given number g in 
the binary system, 8 = bnbp-4... by... bgby. The binary digits of this number, read from right to left, will 
be the components of the desired value y, =(by, by, ..., Dk» ..+» Dn=a» Dy). Precisely by this method the 
values of the vector y are numbered in the table. 


Using (5) for each value yg. of the vector y, we obtain the corresponding constituent units pg Cy). 
However, the expressions for the constituents. pg (y) are easily obtained from the corresponding values y, by 
the following rule: replace each component bj of the vector y by the factor yy if by =1, and by the factor 
y'; (ise., the inverse of the value of y,) if b, =0, The values of 32 constituents pg (y) are given in the third 
column of the Table, 
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Since the components b, of the vectory, can take only two values, 0 and 1, there is no difficulty in 
prepariag cards simulating the vectors y,. For this, it suffices to define the places (positions) on card number 8, 
simulating the vector y,, assigned to the components by of vector y,, and then to make notches or perforations 
only in those positions corresponding to values of by =1. The positions k, corresponding to the components by, 
of vector Yg» Are most expeditiously ordered linearly, placing them, for example, along one of the sides of the 
card in increasing order of k. 
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Fig. 1. Cards 0,1, 2,3,...,16,..., 30, 31, for the synthesis of 
autonomous relay systems. 


Since there is a one-to-one correspondence between the values y, of the vectory and the constituents 
Pg (y) having the same number 8, simulation of the vector y can be considered as simultaneously simulating the 
constituents p(y), With this, a notch or perforation in position k of card number 6 corresponds to the factor yy 
in the constituent p, (y), and the absence of a notch or perforation in this position corresponds to the factor yx 
of the same constituent, Now, if opposite each position k of card number 6 is written the designation Pg of 
the constituent p, (y), cotresponding to the given card, then the cards thus obtained may be used for the synthesis 
of autonomous relay systems in accordance with Formula (4), We immediately add that, technically, it is more 








effective to place, not the constituent designations pg(y) on the cards, but only their numbers 8, with this 
proviso, that the numbers be large enough (but not larger than the notches or perforations) so that when one card 
is laid on another, the number of the lower card might be read through the perforations of the upper card. 


Figure 1 shows the cards, prepared by this method, for the values 6 =0, 1, 2,3,..., 16,..., 30, 31. 


As is obvious from Fig. 1, card 0 (card number 0) has no notches, card 1 has a notch in the first position, 
card 2 in the second position, card 3 in the first and second positions, card 16 in the fifth position, card 30 in 
positions two, three, four, and five, and card 31 in all five positions, 


As is well known [7, 8], the process of changing the state of an autonomous system may be given algebrai- 
cally either by the sequence of values yg, j of the vector y(j), describing this system, or by the sequence of 


constituents pg j (y), corresponding to these values of the vector y, Therefore, the sequence in which the cards, 


8; are placed can simulate the process occurring inthe autonomous relay system, whether we consider these cards 
as simulating the values YB; of the vector y(j), or as simulating the corresponding constituents pg, {y- With 


this, the left sides of the cards 8, containing the notches or perforations, are best considered as models of the 
values Ys of vector y, and the right sides of the same cards, filled with the numbers 6, as the bearers of the 
constituents pg(y), i.e., the number 8 is considered as a shortened notation of the constituent pg (y) itself. Then, 
if card 8; is laid on card 6 ;'so that the left (i.e., perforated) part lies on the right (i.e,, the part bearing the 
number 8 ;) part of the first card,then through the notches of the second card will be seen the number 8; of 

the first card, If we consider the left side of card 6 5' as the model of vector yg } then the left side of card 6 ;', 


lying on the right side of card 8 j can be considered as modelling the product of the vector y, } by the constituent 
Pg {7 i.€,, as modelling the expression pg eh YB i+" If the positions of the cards are reversed, we obtain the 
model of the other expression, pg j “Y)¥p ; 


For example, let card 7 be placed oncard 3 in the manner described, This gives 2 model of the expression 
PsY7=Ps(1, 1, 1, 0, 0] =[ps, ps, Ps, 9, 0), shown in Fig. 2. 


Through the notches on upper card 7 are seen the numbers 3 on positions 
one, two, and three, while positions 4 and 5 of card 3 are covered by the non- 
perforated portion of card 7, This denotes that card 7, lying on card 3 as 
shown in Fig, 2, models the vector [pg, Ps, Ps, 9, 0), i.e., the expression pgyy. 
If cards 3 and 7 are interchanged,we obtain the model, shown in Fig. 3, of 
the expression pry = p7[1, 1, 0, 0, 0) =[py, py, 0, 0, 0), As shown in Fig. 3, 
through the notches of upper card 3 are seen the numbers 7 on positions one 
— and two, and on the remaining positions the numbers 7 are covered over by 
the nonperforated portion of card 3, This denotes that placing card 3 on 
card 7 gives a model of the vector [py, py, 0, 0, 0], i.e., of the expression 


Pry 3- 


Thus, the result obtained by placing one card on another card depends on the order or sequence in which 
e placement is effected, 
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) Let there be given a certain sequence y(j) =yz j of values yg j of the vector y characterizing the states 


of some autonomous relay system, This sequence will terminate, as known from [7, 9], in some periodic process 
or, in particular, in some stable state, i.e,, in the general case the sequence will have the form 


Yor Yores ss Yara (Yours +++ Yes tm—a)? (8) 


where the sequence in brackets is a periodic sequence, 
If cards 8 ; are chosen in the sequence corresponding to Sequence (8), i.e., in the sequence 


Bo» Bi, dine Bi,—1 Bi,» ee) By,4-m—1» B;. (9) 








and laid one on the other such that the right (i,e,, the constituent number- bearing ) part of each previous card 
6, will be overlapped by the left (i.e,, perforated) part of the following card in the sequence 8 then the 


1° 
interlacing of cards 8; and B j+, will model the expression Pa; WY e544 = Pa WY (j+1). j+ 


If the given sequence is just one of several in a given complete set of sequences then, using a superscript 
index r for the given sequence, we have that the card 6 j+,, lying, as described above, on card 8 ;, models the 


expression P ,,(Y) y' (7 + 1). 
] 


The set or, more accurately, the union of all such expressions, obtained from all N,+1 cards 8 je lying 
one atop the other in the fashion described above, in the Sequence (9) is the Boolean sum 


Ny-1 


> Pin) yi + 4), (10) 
j-o  *3 


contained in Formula (4), 


By carrying out such an operation for each of the sequences yD, and finding the union of all N Boolean 
sums (10) thus obtained, we get the right member of Formula (4) and, thus, the sought for function itself, f(y). 
The components f(y) of this function are the conditions for operating the corresponding relays Y, of the 
autonomous system to be synthesized. 


To illustrate the punched card method of synthesizing autonomous relay systems, we consider the following 
example, 


Let it be required to implement the autonomous relay system in which the following processes would occur 
(7, 8}: 


P (Yo) = (Yoo Ya» Ya» Yo» Yas» Yur» Yoo Yo): 


P (Ys) = Yas Yo» (Yaar Yous ++» Yis)s 
P(¥4) = Yas (Yor Yar +++» Yur)s 

P (Yq) = Yo Ys» Yor (Yar Yar ++ » Yo): 
P (Y12) = Yaa» (Yas Yous ++» Yo)s 


P (Y14) = Via» (Yoo Ye:-++» Yu)- 


All these sequences converge to one and the same periodic sequence P(y,) and form a complete set of 
sequences, since they contain all 16 possible states of the four relays Y,, Y2, Ys, and Y,4. 


To find the function f(y) describing the desired autonomous system, we first find the functions f (Wy) 
(r=1, 2,..., 6), in each of which there occurs just one of the given six processes: P(y,), P(y:),..., 
P(yy). The desired function f(y) will obviously be the Boolean sum of the six functions f ) y). 


We first find the autonomous system, the changes of state of which would be described by the sequence 
P (yy), where we require that this process be maximally stable, i.e,, that from all the initial states which are not 
included in the given periodic process, the system will immediately transfer to state yo, the initial phase of the 
given periodic process P (yp). 

By stacking the cards 8; according to the method described above, in accordance with the values yg j 


entering into the given sequence P (y,), we obtain the disposition of cards 0, 1, 3, 7, 15, 11, 9, 8, 0, shown in 
Fig. 4. By considering, in accordance with the conventions we have adopted, that the left side of card 6 ; models 
the vector YB; and the right side of the same card models the constituent Pp We see from Fig. 4 that 


Po¥1 = [Po, 9, 0, 0}, PrYs = (Pi, Pr, 0, 9), 





it 


fo 





Ps¥7 = (Ps, Pa» Ps» 0}, P2Y¥is = (Pr Pa» Po Pads 
PisYi1 = (Pas, Ps» 9, Pasl, PuYo = [Pu, 9, 9, os 
Ps¥s = (0, 0, 0, Pol, PsYo = [0, 0, 0, O] = 


Taking the union of all these expressions, i.e., finding their Boolean sum, we find the desired function: 


£™ (y) = poy + PrYs + Ps¥x + Pr¥1s + Pas¥ar + P1r¥o + Po¥s + PaYoo 


where "+" is the sign for Boolean addition, We note that the last term here is zero, since yg =0 , and hence 
it might have been omitted, 


If we now replace the terms found above by their values, and carry out the addition according to the rules 


for vector addition, we obtain the components f @) (y) of the vector function f @) (y): 
” (y) = Po + Pi t+ Pat Po+ Part Pass fa” (¥) = Pr + Pa + Po + Pass 
fs? (y) = Pa Pr» fa? (¥) = Pr + Po + Pur + Pas: 


By comparing the numbers of the constituents in the expression for the k-th component ft) (y) of the 


function f aad (y) with the digits which are visible in the k-th row of the cards distributed as in Fig. 4, we see 
that the digits which can be seen through the notches in the k-th 
row, distributed in the manner described above, are the 6 V5 of 


a * the constituents of the decompositions of the components he dy). 
5 
We thus arrive at the tremendously simple recipe for deter- 
715) 0) 9 mining the components ic ’ (y) of the function f a (y) describing 


the autonomous system, by the individual possible processes in which 
the states can vary only in a given sequence, Youn" Yon’ 08 sp 





Fig. 4 


> Your (¥ a(n" ces Yon) It suffices for this to choose 


the cards with the same pamber as the vectors in the given sequence, and to place them one atop the other in 


accordance with the sequence "Bt? Be”... , BI, BI,..., Bima so that the right side (containing 


notches or perforations) of each successive card 8 ~ , overlaps the left side [containing the numbers sy” of 


the constituent Pet (y) J of the previous card a ’ . Then the numbers which will be visible through the notches 


or perforations of the k- th row of all these cards will be the numbers of the constituents in the decomposition of 
the k-th component of the function f™® (y). 


By taking the Boolean sum of the constituents with the numbers read off the k-th row of the cards, we obtain 


(r) 


the desired component ra (y) of the function f“(y), By carrying out this operation for all values of k (k= 


=1, 2,..., 1), we obtain all the components of the function f ” (y) and, by its definition, we obtain thus the 
function itself, 


In order to obtain the function f (y) describing the autonomous system in which occur not only the pro- 


cess Ply =, but the complete set of processes given above, it is necessary to carry out a Boolean summation of 
all N functions found in this way, 


In case the given complete set of sequences of state variations in the autonomous system contains a 
sequence leading to some general sub-sequence, the work necessary for synthesizing this system is shortened, 
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since the decompositions into constituents found for the components of the functions corresponding to this general 
sub-sequence maybe utilized for determining the constituents into which the functions corresponding to the se- 
quences converging to thissub-sequence are decomposed, 


The sequences P (y,), P (y4), P(yg), P(yy), and P(y4,) of our example all leadtoone and the same pestodic 
ot rerenty P (yo) and, consequently, among the constituents of the components i? (2) y), 63? y), t+) y), 
£5) (y), and f{*2(y) are contained all the constituents, already found A in the decomposition of the com- 
poh f,' (y). Therefore, for the determination of the components fi) (y), where r= 2, 3, 4, 5, and 6, it is 
obviously not necessary to carry out the operations, described above, with the cards corresponding to all the terms 
of the sequences P (7 A §) . It suffices to carry out these operations only with those cards corresponding to the 


initial terms of these sequences, including also the first term of their periodic sub-sequence, In particular, for 
the determination of the components Ay) of the function 1™y) it is only necessary to use cards 2, 5, and 

















tTela 11, By placing these cards one atop the other by the method described earlier, i.e., as 
i shown in Fig. 5, we get from a consideration of the first row of these cards that, besides 
sya the constituents already contained in fj* (y), this component contains the constituents 
i " Pz and pr, i.e., 
5\ 
# AL (y) = Pa + Ps + 12 (y) = Pat Ps + (Po + Pi + Ps + Pr + Pur + Pus): 














a From a consideration of the remaining three rows of cards 2, 5, and 11, placed as 
Fig. 5 shown in Fig. 5, we obtain, in analogous fashion, that in addition to the constituents 


contained in f Q) yy), (2 ) {y ), and g y), the components My ), Py), and My) 


contain just one additional constituent, namely, they contain, respectively, ps, pz, and ps. Thus, we obtain 


#2 (y) = Ps + AP (y) = Ps + (Pi + Ps + Pr + Pus) 
3 (y) = Pot fs” (y) = Po + (ps + Pr): 
f? (y) = ps +72? (y) = Ps + (P; + Po + Pir + Pis)- 
By a similar use of cards 4 and 9, we find the components £° ) (y): 


PY=atPy), MYw=f \y), 
fs” (y)= #8” (y), 1? (y) = pat f? (y). 
By an analogous use of cards 6, 13, 10, and 1, we obtain the components aM ) yy): 


f(y) =Pet Pot ly), Ay) =ps+ f(y), 
fs? (y) = Pe +f” (y), t (y) = Pe + Pas + 7%” (y)- 


Using cards 12 and 8, we find the components Ay): 


Rw=—My) MM=2 yy), 
- (y)= ? (y), . (y) = Pu + fa a (y). 


Finally, using cards 14 and 9, we find the components ag y): 


Ay) = Put f(y), My) =f), 
f(y) =/2 y). i (y) = pu +2 (y)- 


Taking the union of the decomposition into constituents of the components A (y) for all values of r, 
we obtain the decomposition into constituents of the components f, (y) of the desired function f(y): 





of 


by 


st 


ul 


“~~ oc -*- -_ FY 


Ni (¥) = Po + Pi + Pa + Ps + Pa + Ps + Pa + Pr + Pio + Pur + Pus + Piss 
To (¥) = Pi + Ps + Ps + Pr + Pis + Pris» fs(y) = Po + Ps + Po + Po» 
Ta(y) = Pa + Ps + Po + Pr + Pot Pur + Piz + Pra + Pis- 


Comparing the values obtained for the components of the function f {y ) with the corresponding components 
of the function f (y) obtained in Example 3 of paper [8], we note that they coincide completely, Consequently, 


the desired relay system may be implemented by the same relay contact circuit as in the example cited, i.e., 
by the circuit 


(Y¥, + y,) Ywit+n (Ys + y,) Yo+ yoy¥s + (Ys + yrya) Ve. 


4, Carrying Out the Synthesis of Nonautonomous Relay Systems 





By Means of Special Cards 





Comparing Formulas (6) and (4), we note that the only difference between the first and the second 
formula is this, that in the (j+1}th term y )(j+1) of each of the N given sequences, the values of the vector 
y of the states of the nonautonomous relay system to be synthesized are multiplied, not only by the constituent 
Pax) (y ), constructed from the components of the vector y, but also by the constituent Pal) (x), constructed 

J 


from the components of the vector x of the external stimuli, 














It follows from this that, for the synthesis of 























: 0 O : O ; nonautonomous relay systems, it is necessary to 
0 1 2 iS ok ie have, besides cards 8, simulating the values y 8 
0 f 2 J of the vector y and the corresponding constituents, 
0 1 2 J Pg (y), cards a simulating the constituents p, (x), 
The cards q must contain, at each position, the 
number a appropriate to the given card, The cards 
7 Aa 3} a in no way differ from the right sides of the corres- 
6 Ow Os ponding cards 8, designed for the synthesis of auto- 
6 —_ Ow Ow nomous systems, Briefly put, if the left side were 
Ee 16 Ow Ow clipped from a card g designed for the synthesis of 
| autonomous systems, we would obtain the card a 





having the same number, i,e., the card for which 
Fig, 6, Cards 0,1, 2,3,...,16,..., 30, 31, a =8. 


for the synthesis of nonautonomous relay systems, The cards 6 designed for the synthesis of 


autonomous relay systems may have the same con- 
struction as the corresponding cards for the synthesis of autonomous systems, but in all cases they must be twice 
as wide as the latter, though with the same dimensions for the numbers 6 placed on them, The left sides of the 
cards g must be twice as wide as the left sides of the same cards designed for the synthesis of autonomous 
systems, The reason for this is that in the case of nonautonomous system synthesis there must be visible through 
the notches or perforations of card 6 ;,, two numbers: the number a; in the corresponding position of card a; 
and the number 8 j in the card g j- e right side of the same card, 6 j+1» Must be twice as wide because, in 
addition to the column containing the number 8 j43 in all positions, it must contain (for the synthesis of non- 
autonomous systems) a blank field on which the card a;,, must be placed, 


In order that the numbers a, and 8; which are seen through the notches on card 644, be easily dis- 
tinguished from each other, it is preferable to print the numbers on cards a; and gj either in a different font 
or in different colors, If this is too difficult to perform, and the digits of the numbers aj and 6 j are executed 
in the same font and the same color, then different forms should be used for the apertures on card g for the 
numbers « and the numbers g; for example, retaining rectangular notches for the numbers g, use circular 


apertures for the numbers a, Cards 8, prepared in this way and designed for the synthesis of nonautonomous 
relay systems, are shown in Fig, 6, 














Figure 7 shows one such card g, namely, card 3, on which has been placed one of the cards a, namely, 
the one for a = 5, 


We recall that the numbers § on card g designate the constituent Pg ¥ ), and the number placed on 
card a designates the constituent p,, (x), Therefore, the numbers which are visible on each position of card A, 
with those on the card a lying on it, properly designate the product py &)p, (y) of the constituents 
P q(%) and pg y ). In particular, the pair of numbers 3 and 6 which is seen at any of the positions of cards 
6 =3 and a =5, shown in Fig. 7, designates the product ps (x) psy ) of the constituents pg &) and psy). 


If the right side of card 6;, with card a; lying on it, is covered by the left (perforated) side of card 6 j44 
then the total (overlapped) portion of cards aj , 6; and 6 j44 gives the expression Pay (x )pg {PB sas » In 


particular, if cards §;=3 and aj = 5, placed as shown in Fig. 7,are covered by the method just described by 
card 6 j43 =, then the total portion of all these cards gives, as shown in Fig. 8, the expression ps(x)ps(y )y7 = 
= ps(x)ps(y ) (1, 1, 1, 0, 0) = [py (x)ps ly), psx) ps Cy), ps x)ps fy), 0, 0). 









































or 3X5) 7 116)71@ «4 
315 I) 7 116)7@ «6 
a|5 3\5) 7 316)7\@ « 
v15 7 7 \(2) & 
I|s 7 6 
Fig, 7 Fig. 8 Fig. 9 


If some card aj43 is now laid on the right half of card 6 ;,3, and these are covered by the left half of 
card B j,42 in the manner described above, then the overlapping parts of these three cards give the expression 


Pa; re (x) Po; 44 (y) ¥, jae’ and the union of the intersection of the cards a j41, Bj+1 42d Bj,, gives the 
Boolean sum Pp, (X) Pp, (¥) Yo... + Pa, ,, (*) Pos, (Y) ¥p,,,° 


In particular, if we lay card a;,4 = 12 on the right half of card 7 in Fig, 8, and then cover the right half 
of card 7 plus card 4,4 = 12 by the left half of card 6 j,2=15, then the overlapping portions of these three 
cards give the expression py (x) p;(y)y45, and the union of the intersection of cards 15, a4; = 12,and 7, and 
the intersection of cards 7, aj=5, and 3 gives, as is obvious from Fig. 9, the Boolean sum 


Ps (X) Ps (¥) ¥2 + Paz (X) Po (Y) x 18- 


The first three components of this vector surn, as is obvious from Fig. 9, are identical. and equal the sum 
Ps (x) ps Cy) + pas (x) py Gy), and the fourth component equals the expression py2(x ) p; (y). 


In order to obtain, by means of cards aj and § ;, the values of the Boolean sums 


Ny+1 


D, Peter (x) Py (y) ¥ i + 4), (a) 
j-o ? j 


which are parts of Formula (6), it is necessary to proceed analogously, It is necessary to take all the cards ay, 
and 64 corresponding to the terms * and ¥,(n) of the given r-th pair of sequences x )(j) and yp, and 
j 


to place them, one atop the other, in the sequence 
at”, af”, pt, al”, ..., 807, al, ..., Ka, alta, BW (12) 





ct. fF &® #4 


al 


ra 





(r) 
Oy 


so that the left (perforated) part of each successive card ye covers the right part of the previous card 8 


together with the card at » » lying on the blank field of the latter, i.e., analogously to that shown in Fig, 9, 


where a fragment of such a disposition of cards a; and 8 Ns shown, Then, the numbers which are visible through the 
the rectangular notches and circular apertures in ards Bit 1) will model the vectors Patt) (x) p a y™ G+), 


and the entire set of pairs of numbers ps and a? visible through the rectangular notches and edit apertures 
of all N, + 1 cards, positioned in the way described, models the Boolean sum of all these vectors, i,e., Sum (11), 


By taking the Boolean sum of the paired products p ap ™ )p aD of constituents, where the numbers 


B “ and a - are visible through the rectangular notches and circular apertures of the k-th row of cards By, 
placed according to the method described above and interwoven with cards a; in Sequence (12), we obtain 
the kth components of the vector sums in (11), 


If we carry out the operations described above on all N given pairs of sequences x ® (j) and yo, and 
find the union of the results obtained, we get the right side of Formula (6), i.e, the values of the desired function 
f(x,y). Each k-th component of this function is the Boolean sum of the k-th components of all N Boolean 
vector sums (11), i.e., the union of all paired products of the constituents Pag) “Po the numbers of 


which may be read through the apertures of the kh rows of cards rod in all the sequences of (12), By definition, 
this union is the set of conditions for operation of relay Y, in the nonautonomous system to be synthesized, 


We mention in conclusion that the synthesis of autonomous relay systems can be carried out with the aid 
of the same 6 cards designed for the synthesis of nonautonomous sequential relay systems, For this, the cards « 
are not necessary, and the cards g, designed for the synthesis of nonautonomous systems, are used in exactly the 
same manner as the cards designed for the synthesis of only autonomous systems, Here the circular apertures in 
the cards and the blank fields on them play no role, and no attention should be paid to them, 


For the synthesis of autonomous relay systems, it is necessary to have no less than one complete set of 
cards 8, containing all 2” cards, beginning with card 0 and ending with card 2"— 1, However, it is convenient 
to have two complete sets of these cards so that some _ periodic sub-sequence,one of which terminates every 
possible sequence of states in an autonomous system, be terminated by the same card ai? with which it began, 


For the synthesis of nonautonomous sequential systems, in addition to two complete sets of cards 6, designed for 
the synthesis of these systems, it is necessary to have several complete sets of the cards a. For convenience 
in changing cards a upon going from one sequence of type (12) and its corresponding disposition of cards 8 ; 
and a , to another, it is better to make the cards a somewhat longer than the cards § so that they will project 
beyond the edge of cards g, as shown in Fig. 9. 


Sets of cards designed for the synthesis of relay systems constructed of n relays and depending on m in- 
dependent variables may also be used for the synthesis of systems constructed of fewer than n relays and de- 
pendent on less than m variables, In the examples considered above, the cards that were designed for five re- 
lays and five variables were used for the synthesis of autonomous systems constructed with four relays, and for 
the synthesis of nonautonomous systems containing four relays and depending on four independent variables, 


The punched card method of synthesizing sequential relay systems may be used with manual placement 
of the cards and visual reading of the numbers through the apertures in the cards for the synthesis of relay systems 
where the total number of variables (dependent and independent) does not exceed 8 or 9, 


If mechanized, the punched card method could apparently accommodate twice the total number of relays 
and independent variables in the systems to be synthesized, 


To use the punched card method in practice, it is necessary to provide for the most efficient construction 
of the cases in which the cards @ and g are to be stored, the construction of the framework on which the 
assembly of sequences of type (12) is to be carried out, the construction of the cards themselves to guarantee 
rapid assembly and storage of the assembled sequences of cards, etc, 
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ON THE STABILITY OF PERIODIC REGIMENS IN NONLINEAR SYSTEMS 
WITH PIECEWISE LINEAR CHARACTERISTICS 


M, A. Aizerman and F, R, Gantmakher 


(Moscow) 


A way is described which allows one to find the equations of a linear approximation 
solving the problem of stability of a periodic solution in a system with a piecewise linear 
characteristic, 


1, As in works [2, 3], we consider an automatic control system described by the equations 


t= Dac; +2, +F,) (i =4,...,m), (1) 


j=1 


where a;; and ){ are given numbers (some of which are zero), F;(t) are periodic functions with common period 
r (in particular, constants, including zero) and f(z) is the given piecewise linear function. 


For concreteness, let the transitions from one line in the characteristic f(z) to another occur at the 
moments when z, first attains the values og (a=1,...,1). We assume that there has been found a periodic 
solution of Equations (1) with period r, for example, as cited in[2, 3), and that it is required to investigate its 
stability (in the sense of Liapunov), For this purpose, the theorem proved in [4] may be used, This theorem 
generalizes to systems of differential equations with discontinuous right members,the classical theory of stability 
by linear approximations, developed by Liapunov, applicable to systems of differential equations with analytic 
right members. 


2, Before getting back to System (1), we first consider the more general equations 


z =f, (21,..., zy, t) (émi,..., n), (2) 
where all the f; are functions periodic int with common period 1, or functions not depending upon t. 


Let (zy, ...» Zs t) space be sectioned by the surfaces F,, (74, ... » Zs t) = 0 into regions H,, (a= 
=1,2,...,). It is assumed that in each region H,,, the given curves are sufficiently smooth; upon a transition 
through a surface F.,, there may occur discontinuities both in the functions fj, themselves and in their partial 
derivatives, The integral curves of Equations (2) are continuous, however, but are subject to breaks (discontinuities 
in the derivatives) at the surfaces F,,. . 


Further, let z; = Z; (t) be the periodic solutions of System (2) intersecting the surfaces F,, =0 at the moments 
a 


Together with System (2), we consider the linear equations, the right members of which are given for all 
regions H,,: 


t=t 


n 
a, = Dy b4;2; (i=1,...,a), (3) 


j=l 
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of 
where },; == oe are periodic coefficients with common period r, defined for arbitrary t except t=tg. 


We add to Equations (3) linear relationships defining the discontinuities in x; for each moment t =t,: 





n n 
z, (t, +0) — 2, (t, —0) = &, 2 hr x, (t, —0) = &, dat 2, (t, +0), (4) 
k=l k=l 
where 
oF, 
Oz, 


jaF,\= 
(a) t=2 (tq) 
Equations (3) and Relationships (4) jointly determine the discontinuities (for t = t,,) in the integral curves, 
We call the set of Equations (3) and Relationships (4) the linear approximation of Equations (2) for their periodic 


solutions 2; = a (t), Then, with certain additional limitations* imposed on the functions fi and F,,, we have 
the following theorem *® [1]: 


If the null solution x;=0 of the linear approximation system(3) and (4), is asymptotically stable, then 
the periodic solutions zj =z, (t) of Equations (2) are also asymptotically stable. 


3, Returning to System (1), we note that if, in the interval t,,., <t <tg, the section of the characteristic 
has the equation 


f (21) = k,21 + 8, 


then Equations (3) for System (1) may be written in the form 


n 
By = DY) 045 2; + A 2 (t) 2, (3") 
jimi 
whens 1 (t) is the periodic piecewise constant function 


LQj=k, (t,<t<t,, a=4, 2,...). 


The equations for the discontinuity surfaces for Equations (1) we write in the form zy — 0g = 0. 
Therefore, Conditions (4) lead to the conditions 
x, (t, +0) — 2, (t, —0) = RF x (¢, —0), (4") 


‘ 
R=-a (so 
le 2% _. 


Here € is the magnitude of the discontinuity of the function f(z) for 21 =2Z4(t,), and ¢=1lore=-1 
depending upon the direction of intersection of the surface by the trajectory. 
Integrating Equations (3") for t,.4 <t <tq, and adjusting them at the limits of the period by taking (4") 


into account, we find the linéar relationships expressing the coordinate values at the end of the period, x; (ty+ 1), 
in terms of their values at the beginning of the period, xj (t)). 


‘We write these relationships in matrix form thusly: 


* These limitations are not discussed here since they are always satisfied for System (1) which is of interest 
us to, 
¢* The analogous theorem on instability is also proven (see [1]). 
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2 (to + T) = Ux (to), 


where U_ is the constant transformation matrix, 


Thus, for the stability of the periodic solution being considered, it suffices to require that the characteristic 
numbers of matrix U , i,e,, the roots of the characteristic equation det (U- AE) =0 be distributed in the unit 
circle, 


4, Certainly, the question of stability of the periodic solution of System (1) could be solved by direct ad- 
justment, For this, Equations (1) would have to be integrated over a period, with individual integrations necessary 
for each step ty., <t <tg, to determine in this way the transformation of the quantities 2; (ts) to z(t) + r) 
(i=1,..., 0), and to investigate the stability of the invariant points of the transformation corresponding to the 
periodic solutions, Such an investigation, in its turn, requires the expansion of the solution obtained in a series 
of powers of the variation of the initial deviation and of time ty, a discarding of the nonlinear terms, the elimi- 
nation of the variation in time tq, and, finally, the determination of the transformation matrix of the coordinate 
variation during the period. 


The method described above, using the theorem on the linear approximation for piecewise-smooth systems, 
allows this computational work to be shortened, since it allows the introduction of the variation of time ty, to 
be avoided, and allows the immediate writing of the necessary linear relationships 
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THE TRANSFER FUNCTION OF A DC MOTOR CONTROLLED BY 
VARYING THE EXCITATION VOLTAGE 


E,. L, Urman 


(Moscow) 


The transfer function is derived for a dc motor with a supplementary series excitation 


winding. 


In many systems of automatic regulation and control, dc motors controlled by varying the excitation voltage 
are employed, For computing and investigating the stability of such systems,an expression for the transfer function 
of the motor, well known in the literature (for example, [1]) is generally used, This expression, and the corres- 
ponding schematic design in which the motor takes the form of a network with standard links and rigid connections 
between them, are correct in the cases when the motor has one excitation winding. In many cases though, an 
additional series excitation winding is employed in order to guarantee the required starting characteristics, The 
transfer function for such cases cannot be assumed to be the same as for motors without the series. winding. 


Below, an expression is derived for the transfer function 
of a motor which, in addition to the electromechanical time 
constant and the time constants of the excitation winding circuit 
and the armature circuit, takes into account the action of the 
series excitation winding, and the corresponding schematic 
circuit is given, 


The circuit for switching on the motor is given in Fig. 1. 
The motor armature is connected to a dc voltage source of 
magnitude Ua. 





For the circuits of the basic excitation winding and the motor armature, we may write the following 
equations: 


do’ 
Up = phe + veg, (1) 
di do’ 
Un =e t+ ig(Rat Ry) + Lag + “oa (2) 


The equation for the motor shaft moment has the form. 


iG MM, (3) 
In these equations the following nomenclature was used: ug is the voltage on the terminals of the basic 
excitation winding, i, and ig are the currents in the armature and basic excitation winding, e is the counter- 
emf, Q is the nominal speed of armature rotation, Mj, and Mg are the turning moment of the motor and the 
impedance moment, wr and ws are the numbers of turns in the basic and series excitation windings, Ra, Rg, 
and Rg are the winding resistances of the armature, basic excitation winding and series winding, L, is the in- 
ductance of the armature winding, J is the armature's moment of inertia and @ is the total flux linked with the 
basic and series excitation windings. 
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For the turning moment and counteremf we may write 
c= kyOQ. (5) 


In Equations (4) and (5), @ is the effective resulting flux threading the armature winding, The total and 
the effective flux are related by the expression [2] 


©’ = O[1 + y,(e—1)}, (6) 


where g is the leakage coefficient of the poles and y is the coefficient of leakage flux yr with the ex- 
citation winding, Taking a=Ys (o— 1), we may write 


@’ = O(1 +a). (6") 


Replacing the magnetization curve of the motor by a rectilinear characteristic, and ignoring the reaction 
of the armature and the eddy currents in the massive parts of the system, we write the expression for the total 
flux in the form 


a Fe. 
where Lg and Ig are, respectively, the inductances of the basic and the series excitation windings, 


Since it is assumed that there is no leakage between the basic and the series excitation windings, we have 


2 
Ws | 
that Ls* Le We rn 


We write the differential equations for the motor in incremental form, expressing the variables ug, ig, 
ig, Q, and @ in terms of deviations from the steady-state values: 


Q= 2+ AQ, D= @, + AD. 


The incremental equations, in operator form, will be: 


Aug = Ag Ret wp(i + «) pA®, (8) 
0 = ka (DAN + NA) + Ai, (R a+ Ry + L ppAi. gt wg (1 + a) pA®, (9) 
JphQ =k, (DpAi, + ADi, ), (10) 
L 
(1 +a) A® = Eni 4 <8 ai . (11) 
E S 


From Equations (9) and (10) we obtain the differential equation relating the incremental armature current 
with the incremental speed 








[Fats + pit +2) wg + Ps Jaa — 
® R (12) 
= PF | hay + p(t +) ay — a ae a as Ai,. 


J(R R 
We introduce the notation: T,= re is the electromechanical time constant of the motor, 
1 
La 


TA" ek is the time constant of the armature winding, Ts = csittligieee is the time constant of the 
A 
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w T (R gt BR 
series winding, 7, = 0 t= a ¢= a i g is the relative voltage drop in the armature circuit, 
Using this notation we may write, instead of (12), the following equation: 


AQ 1 (pT, + 2)e Ai 4 
Tet eTylt+ +9) Py) go Te Ta +(8 +2) ay} 77 ae 





We also obtain from Equations (9) and (10) the equation relating the incremental effective flux with the 
incremental armature current: 


Aia i— A® 
[1+ pT, (1 + pT,)] Ta, = —{ +— = pT, {1 + (1 + «) pT || o, : (14) 





Finally, from Equations (8), (11), and (14) we obtain the relationship between the incremental excitation 
voltage and excitation current: 


{+ pT) it + PTY + PTA) + PT, [1 +. + @) PT) +% p=} Aig 


1 (15 
= Fg {tt PP a(t + PT + xT lt + 1+ @) Ty) +% 7} Aug 





On the basis of Equations (11), (13), (14), and (15), it is possible to set up the structural schematic for the 
motor (Fig. 2). 
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Fig, 2 


Since the action of the series winding was taken into account, the structural schematic for the motor is 
more complicated than the one obtained when only one excitation winding is present. The structural schematic 
expresses the fact that when the additional series winding is present in the motor, it is impossible to present the 
circuit as one consisting of standard links connected in series, The transfer functions of the links in the structural 
schematic shown in Fig. 2 have the following expressions: 
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The over-all transfer function of the motor, if we consider the input coordinate to be the voltage in the 
excitation circuit and the output coordinate to be the nominal speed, has the expression 


AQ 
W, Ciel cqcaeum inne 
mi?) Rup = 
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, (21) 
®, 
"(1 + pT g) [1 + PT y (1 + PTA + PT yt +1 +o) pTy]x +x 7— 


If the regimen of operation of the motor is varied, a host of quantities entering into (21) have their values 
changed, The quantities dependent on the mode of operation are the nominal speed {p, the effective flux 4», 
the electromechanical time constant Ty, and the relative voltage drop across the armature ¢, 


It is possible to detect the influence of a series excitation motor winding on system stability in the example 
of a system for stabilizing the frequency y of the output voltage of an electromagnetic transformer, consisting 
of a de motor and an ac generator, Figure 3 gives a block schematic of the electromagnetic transformer with 


frequency stabilization, 
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Fig. 3. 1) is the measuring organ; 2) is the amplifier; 3) is 
the motors; 4) is the generator and 5 ) is a differentiator, 


For setting up the differential equations for the system with frequency stabilization, we consider that the 
measuring organ is noninertial, but the amplifier is an inertial link with gain of k, and time constant T,, In the 
motor, for simplicity, we ignore the electromagnetic time constant of the armature (T, = 0) and the pole leakage 
(a =1), In this case, the expressions for the transfer functions of the system links take the following form: 


measuring organ 





Ws (p) = ks, 
amplifier 
a 
W (p= a 

motor 

i—2¢ 

To—< tPT s 

Wy (Pp) oe iy e’ 
(1+ pF ) (1 + pT) + pT, (1 + pT.) * +%7—; 

generator 


Wg (p) = kg, 








The transfer function of the system, opened at the amplifier input, has the form 








W (p) = Wi, (p) We (p) Wy (Pp) We (p) + Wa (p) Wy, (2) (22) 
. . Tonct rT 
ak Ks g i—es Pf x 
4 (1+ ple) (1 + Ty) + pT, (1 + pT) x +%7— 
ss (23) 
(1 + pT) (1 + pr) ° 





In the example considered, T,; = 0.0683 sec , T)y = 0,118 sec , Ty = 0.00044 sec , 


*=0.2 sec, T, = 0.06 sec , ¢ = 0.124, x = 1.57, 
volts radians/sec 
Wty = Sediangeed Mhka= 10, kM= Sop, kyg= 3. 
We get the expression for the amplitude- phase characteristic of the system by substituting p = jw in (23): 


144 + 0.074 jw 40 je 


W (je) = — TF — 0.0088 wo + 0.B72}@ ~ (40.00 jm) (14 0.2]e) ° ta 





In the case where there is no series excitation winding (ws =0, Tx =0, x =0), the expression for the 


amplitude- phase characteristic will be 


144 40 ja (25) 


W (jo) = — G + 0.0683 ja) (140.118 jo) ~~ (1 4+ 0.06 ja) (1 +0.27e) - 
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Fig. 4 





2. 





Figure 4 gives, on a logarithmic scale, the amplitude- phase characteristics computed on the basis of (24) 
and (25): a is with a series excitation winding present, b is without the series winding. The curves of Fig, 4 
show that the series excitation winding of the motor allows stable system functioning to be achieved, For the 
given values of the system parameters, the system will be unstable in the absence of the series winding. 
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ALL-UNION CONFERENCE ON MAGNETIC ELEMENTS IN AUTOMATION, 
REMOTE CONTROL AND COMPUTING TECHNOLOGY 


During November 25-30, 1957, in Moscow, there was held an All-Union Conference on magnetic elements 
in automation, remote control and computing technology, organized by the Institute for Automation and Remote 
Control of the Academy of Sciences, USSR, and by the Board (Komissiia) for magnetic amplifiers and contactless 
magnetic elements, The aim of the Conference was the discussion of the state of scientific problems in the 
theory, design and application of magnetic elements in automation, remote control and computing technology. 


More than 800 participants took part in the Conference, representing 246 organizations, including the 
institutes of the Academy of Sciences, USSR, the Academy of Science, Ukrainian SSR, the Azerbaidzhan SSR, 
the Latvian SSR, the Ural branch of the Academy of Sciences, USSR, research and educational institutes, building 
and planning organizations of various branches of industry, Representatives of the Czechoslovakian Academy of 
Science, the Polish Academy of Science, the Bulgarian Academy of Science and the Chinese Academy of Science 
participated in the Conference, 


Seventy papers and communications were presented, drawing lively interest from the Conference participants, 


More than forty people came forward in the discussions, The work of the Conference was carried out in two 
sections: magnetic amplifiers and contactless magnetic elements, 


B, S, Sotskov, President of the Board for magnetic amplifiers and contactless magnetic elements, mentioned 
in his opening words that magnetic amplifiers and contactless magnetic elements play a large role in the de- 
velopment of the technological means of automation, However, the scope and tempo of their production and use 
still do not satisfy the requirements of our rapidly developing national economy. 


B, S, Sotskov, in his paper, “Modern state and problems of the development of magnetic elements for 
automation, remote control, and computing technology,” further stated that magnetic amplifiers and semicon- 
ductor devices do not compete with one another, The appearance of crystal triodes not only did not limit the 
domain of application of magnetic elements and, in particular, of magnetic amplifiers but, on the contrary, 
created favorable conditions for their application, The speaker backed up his opinion with many examples of 
new uses for magnetic amplifiers and of their use in conjunction with semiconductors, In conclusion, the speaker 
voiced his surmise that, as a result of the Conference, there would be a significant expansion in the necessary 
work by the organization of extensive production of various types of magnetic elements, 


In his paper, “Dynamic characteristics of electrical circuit elements,” K, M, Polivanov presented a new 
method of representing dynamic characteristics in the form of hypersurfaces, The use of actual dynamic 
characteristics allows the computation of surface effects to be carried out approximately, taking into account 
viscosity, the influence of which in materials with approximately rectangular hysteresis loops is dominant in 
regions of strong fields, 


R, V. Telesnin gave a paper on the subject “On the influence of magnetic viscosity on the processes of 
core switching.” 


In his paper, “Certain factors influencing the static and dynamic characteristics of toroidal cores," M, A. 
Rozenblat showed the influence on the static characteristics of magnetic amplifiers of eddy currents and the geo- 
metric relationships of the cores, 


During November 26-29, the work of the Conference was carried out by sections, 
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A, Section on Magnetic Amplifiers 








The paper by N, P, Vasil'evaandO,A,Sedykh, “Unified method of designing magnetic power amplifiers,” 
was devoted to the design of magnetic amplifiers with output power from several milliwatts to tens of kilo- 
watts, As the basis of their calculations the authors placed the functional dependence of the minimum specific 
volume of the cores, i.e., the volume per unit power under optimum working conditions, on the strength of the 
magnetizing field, This dependence is determined solely by the cores’ magnetic characteristics and, for various 
designs of the magnetic amplifiers, can be expressed by equations differing only in their constant coefficients. 
The optimal geometry of the cores for magnetic amplifiers is determined from conditions for obtaining minimum 
weight, or dimensions, or cost of the amplifier, The simple engineering method presented allows magnetic 
amplifiers to be designed according to arbitrary schemes, including ac amplifiers and high-speed amplifiers, 


In his paper, "On the design of ac choked magnetic amplifier circuits," N, A, Kaluzhnikov dwelled on the 
question of designing choked magnetic amplifiers with active inductive and complex ac loading, and also for 
the various forms of loading connected to the rectifier output. The speaker presented graphic computations for 
various loads and also for various loading characteristics of the magnetic amplifier. 


N, N, Aleksandrov, in the paper, “On the geometric locus of points on the simultaneous magnetization 
curves determining the functioning of a two-stage magnetic amplifier with constant load,” extended the method 
of computing the characteristic of a single-stage amplifier by means of a loading ellipse to two-stage magnetic 
amplifiers, 


L, V, Safris gave a paper on the subject “Computation of magnetic amplifiers by linearized magnetization 
characteristics." The magnetization characteristics of saturated choke cores are, by this method, described by 
linearized equations in which the influence of the constant and the variable components of the magnetic field 
are taken into account, The equations obtained render possible the establishment of the connection between the 
controlling current and the current and emf in the ac winding, as well as the input inductance of the saturated 
choke, Conversion schemes were given for various types of magnetic amplifiers, The method of linearization 
allows uniform computation for the steady-state and the transient regimens in a number of basic circuits for 
magnetic amplifiers with cores of various properties, 


T. Kh, Stefanovich, in his paper, “On the design of self-saturating magnetic amplifiers,” presented a 
theory of functioning of magnetic amplifiers with internal feedback connections (amplistats) without taking into 
account the width of the hysteresis loop of the core material, and also gave a design method guaranteeing the 
functioning of the amplifier on the linear portion of its characteristic, The design is carried out using specimen 
experimental characteristics of an amplifier prepared from the given magnetic material according to the same 
design as the planned amplifier, 


In the paper of V, L, Benin and I, N, Senkevich on the subject “Special features in the design of magnetic 
amplifiers in reversing circuits for controlling short-circuited, three-phase asynchronous motors,” the method, 
developed by the authors, was presented for computing the parameters of magnetic amplifiers, taking into account 
all the fundamental peculiarities of their functioning in three-phase systems with the condition of minimum 
outlay of active materials, Based on the method of symmetric components, quantitative relationships were ob- 
tained characterizing the functioning of magnetic amplifiers in various circuits, 


The paper of R. Kh, Bal'ian, “Engineering design of chokes for magnetic amplifiers with feedback,” was 
devoted to the design of magnetic amplifiers with feedback for a given load, guaranteeing minimum dimensions, 
weight or cost, The design method is applicable both to single-stage and to two-stage amplifiers with arbitrary 
load characteristics, The fundamental regularities of optimal core geometry were given, and the peculiarities 
of the thermal regimen in choked magnetic amplifiers were considered, as well as the choice of the admissible 
current density in the windings, 


In the paper of M. A, Boiarchenkoy and N, P, Vasil'eva “Fast-acting magnetic amplifiers,” the functioning 
of the fundamental designs for fast-acting magnetic amplifiers was analyzed, and the conditions were derived 
for obtaining fast action in arbitrary amplifier designs; a series of designs of magnetic amplifiers illustrating the 
application of these conditions was presented, The authors produced experimental data and approximate com- 
putations for a fast-acting amplifier with a dc output, and also data for the conversion of ordinary designs with 
de output to fast-acting ones, 








V. I, Molotkov, in the paper “Investigation of the functioning of a magnetic amplifier with a capacitor 
filter at the output of a rectifier bridge,” using analytic methods for investigating the functioning of a magnetic 
amplifier in the transient and steady states, derived the equation for determining the boundaries of the possible 
regimens of functioning of the amplifier. 


N, N, Aleksandrov presented a communication entitled “Peculiarities of functioning of a magnetic ampli- 
fier with the load connected through a rectifier bridge and inductively connected to short-circuited coils." As 
such short-circuited coils appear the contours of the eddy currents in the massive magnetic circuits of electrical 
machines, their commutating sections, or the working contours of the windings in bridge magnetic amplifiers in 
successive cascades, 


The communication of A, S, Bogoslovskii,"Saturable reactors in circuits controlled by high- power three- 
phase semiconductor rectifiers,” gave the analysis of circuits of three-phase semiconductors with saturated 
chokes, and the results of experimental tests on them, 


In the paper of A, L, Pisarev, “Investigation of a magnetizable reactor with a rectifier load circuit, con- 
trolled from a dc source,” there were investigated the processes leading to the saturation of the choke by purely 
active loads with rectifiers in the load circuit, and the functioning of the choke was investigated for approximate 
magnetic characteristics of the material, considering it as being three segments, and as being a parallelogram. 
The author gave an analysis of the functioning of the saturated choke with complex loading, A method was 
given for computations, allowing a step-by-step construction of the transients for a falling load current, 


In the paper of A, A, Golovan, “Magnetic amplifiers with high-ohmic loads,” there were considered de- 
signs of magnetic amplifiers working effectively with loads on the order of tens and hundreds of thousands of 
ohms; with a source supplying voltage pulses alternating in sign, there was investigated a one-half- period non- 
reversing amplifier design, and results of experimental investigation of this design were given, 


V. 1, Minkevich, in his communication, “Certain questions in designing summing magnetic amplifiers 
with low sensitivity thresholds,” presented a method for reducing zero drift in a differential magnetic amplifier 
by introducing an active impedance into one of its arms. The speaker also considered the possibility of employ- 
ing one-half- period circuits in summing amplifiers. A method was presented for sorting cores according to their 
magnetic properties by the use of a series oscillograph, 


In the paper of A, M, Bamdas, V. A. Somov, and A, O, Shmidt, “Transformers and stabilizers, regulated 
by magnetizable shunts,” there was given a sketch of the principles of designing and building regulatable trans- 
formers by means of shunts, considered as elements in the theory of transformer regulation, and data were pro- 
vided for the carrying out of the construction. 


I, G, Gol'dreer and Iu, V, Afanas’ev presented the paper “The theory of even-harmonic (phase- pulse) 
magnetic amplifiers,” in which were considered the specific peculiarities of functioning of even-harmonic 
magnetic amplifiers from the point of view of the general theory of amplifiers, Presented in the paper were 
methods of selecting regimens of functi- ning, and basic parameters of magnetic amplifiers. 


V. A. Oleinikov, in the paper “Investigation of the work of magnetic amplifiers in sampled-data control 
systems,” employed a graphico-analytic method of computation for the construction of the transient response, 
and dealt with the possibility of obtaining forced transient responses, The pulse elements in the magnetic 
amplifier are here considered as equivalent to inertial links, The author showed that the stability and speed of 
action of a servosystem, working in a pulse regimen, are increased in comparison with ordinary servosystems. 


F, 1, Kerbnikov and M, A, Rozenblat, in the paper “Magnetic modulators with lateral excitation,” gave 
a method for designing modulators with mutually perpendicular fields, working either on the principle of fre- 
quency doubling or with the output at the fundamental frequency of the supply source, Such modulators have 
high sensitivity and stability, the transient response in them lasting for one period of the source frequency, The 
authors made recommendations on the use of modulators in devices for automation, 


V. 8. Volodin and G, V, Subbotina, in the paper "Modern applications of semiconducting triodes and 
magnetic amplifiers,” gave a survey of the existing devices with the contemporary application of semicon- 
ducting triodes (transistors) and magnetic amplifiers, and gave an analysis of a design for phase-sensitive recti- 
fiers for the cases when certain types of magnetic amplifiers are used, The possibility was demonstrated of 
increasing output power and speed of action by the modern application of transistors and magnetic amplifiers. 
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R, A, Lipman, in the paper “Magnetic amplifier control by means of transistors," gave computed relation- 
ships in fast-acting magnetic amplifiers, controlled by transistors, and provided several circuits for cascade 
connections of transistors and magnetic amplifiers, 


B, A, Mitrofanov and N, I, Chicherin, in the paper “Applications of semiconductor diodes and triodes in 
amplifiers of servosystems,” dealt with the most widely used elements of servosystems, their circuit connections, 
and also provided data on servosystems with type ADP motors with powers of 4 and of 50 watts, 


G, I, Shapkaits, in the paper “Experiments in the use of two-stage magnetic amplifiers in circuits with 
electron tubes and semiconducting devices," spoke of the modernization of magnetic amplifiers carried out for 
the purpose of using them in circuits for computing machines, and also of attempts to create ac bridge magnetic 
amplifiers with triode excitation and motors at the outputs, 


Ia, M, Kleiman and L, S, Koniaeva, in their paper, “Circuits with magnetic amplifiers for the control of 
the main drive of reversible roller motors," presented interim results of the work on creating circuits for con- 


trolling main drives of blooming machines with magnetic amplifiers as pre-amplifiers to electromagnetic ampli- 
fiers, 


The representative of the Chinese Academy of Science, Professor Ming Nai-ta, presented a paper devoted 
to the synthesis of electric circuits with dissipative elements, 


Ia, B, Rozman, in the paper “Magnetic amplifiers in regulated de electric drives for machines,” dealt 
with the development and installation of dc motor regulation with application made of magnetic amplifiers and 
semiconductor magnetic amplifiers, and of their advantages and disadvantages, and provided results of the testing 
of foreign regulated drives by magnetic amplifiers. 


The paper of la, N, Shtrafun, “The use of high- power chokes and transformers with magnetization in ex- 
citation systems of power turbogenerators, " was devoted to the development of a number of excitation systems 
using contactless magnetic elements, and also to high-power semiconductor rectifiers, The speaker related that 
the results of field investigation of the new turbogenerator excitation system at powers of 30,000 kilowatts 
showed a sufficiently high quality and demonstrated the efficacy of using it for series high- power turbogenerators, 


E, F, Stepura and V, V, Semenov, in the paper “Some applications of saturated chokes in electromechanical 
automatic control systems," dealt with work on the creation of contactless nonlinear and computing elements 
for use in electromechanical automatic control systems for the purpose of improving their dynamic characteristics. 
Such elements possess simplicity and high reliability, and their errors do npt exceed 2-5%, 


V. G, Komar and P, P, Mokhov gave a paper on "High-power static ac voltage stabilizers with magnetic 
amplifiers,” in which they dealt with developments of high-power three-phase voltage stabilizers at 300 and 500 
kilowatts for searchlight filament sources, In the paper a design was given for connecting saturated chokes, 
guaranteeing a high-power coefficient and an increased coefficient of efficiency. 


F, F, Sokolov, in the paper “A single-phase voltage stabilizer with saturated chokes for powers of 1.25 
kilowatts,” spoke of the selection of the sensing, amplifying and executive organs of the stabilizer, and gave a 
method for determining the computed power of various stabilizers in order to compare them, 


B. A, Mitrofanov and N, I, Chicherin gave a paper on "Servosystems, ac, using magnetic amplifiers and 
ac motors with powers of 0,5-1 kilowatt," devoted to the development of ac servo drives (two-phase and three- 
phase), and giving the basic results of the work cited, 


S. S, Roizen, in his paper, “Rolling mill electric drives with magnetic amplifiers,” dealt with certain cases 
of the use of magnetic amplifiers in rolling mill electric drives, The speaker cited the efficacy of use of three- 
phase two-stage designs for the output magnetic amplifier cascades, and also recommended the use of higher - 
frequency magnetic amplifier supplies for the improvement of the dynamic characteristics of the electric drives, 
An example was given in the paper of the use of magnetic amplifiers in the electric drives of a reduction turbo- 
rolling mill containing twenty stands, each with an individual dc motor, 


B, Section on Discrete Magnetic Elements 





A, 1, Bakhir, in the paper “Durable and operative memory devices of ferrites and semiconducting elements,” 
considered the over-all design for constructing a device in which the writing of signs is carried out by means of 
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leads drawn through ferrites according to a given rule, while for choosing the address and reading the information 
written are used decoders and shapers, implemented completely by ferrites and semiconducting elements, 


In the paper of Mr, Ovcharenko, “Basic logical circuits, executed by magnetic elements and transistors,” 
there were considered the different characteristics of computing cells, made up of ferrites and transistors, and 
various types of such logical circuits were given. 


la, A, Khetagurov, in the paper “On one design for filling devices made of ferrites," dealt with the design 
of a filling device which might be used in hysteresis machines, and with the peculiarities of developing such 
designs and its basic blocks, The speaker gave a detailed treatment of the frequency characteristics of the design, 
and gave results of tests of the design with various types of ferrites, and also gave data on a small-size memory 
device executed with semiconductors, 


In the paper of K, G, Mitiushkin and V, A, Zhozhikashvili, “The use of contactless magnetic elements 
in tele-control devices,” there was noted the superiority of contactless magnetic apparatus for the construction 
of tele-devices by the method of continuous transmission of communications, and two types of contactless de- 
vices for tele-control were characterized: time the distributing device works sporadically and continuously 
acting devices with a pulse distributor source from a common system with synchronous frequency. 


V. A. Ordynets, in the paper “Computing cells of magnetic elements,” spoke of a series of circuits con- 
structed of ferrite cores, in particular, of circuits of the choke type, and of circuits using the saturation effect of 
the individual portions of the magnetic conductors, and gave practical circuits for registers, adders, triggers, valves, 
etc, 


The paper of N, M, Brusentsov, “Digital elements of the type of magnetic amplifiers with a source of 
current pulses,” was devoted to the use of such amplifiers as elements of high-speed digital machines, The author 
considered the advantages, disadvantages and characteristics of elements of this type, gave a method for designing 
such elements, and gave examples of circuits for ternary algebraic adders, ternary inverters, ternary decoders and 
switches, 


In the paper of E, F, Berezhnii, “Magnetostrictive elements in computing technology devices," there were 
considered the design and construction of magnetostrictive delay lines for memory devices, The author provided 
results of an investigation into the influence of temperature on the stability of the delay time and also answered 
questions on compensation for the variation in delay time due to temperature, The characteristics were given 
for operative registers composed of magnetostrictive delay lines, 


In the paper of V, V. Bardizh , “Questions of pulse magnetization of magnetic cores," the working con- 
ditions were investigated of ferrite cores with rectangular hysteresis loops in matrix-type, magnetically operating 
memory devices, and the basic requirements for increasing the speed of action of devices of this type were 
formulated, 


E, I, Gurvich gave a paper on “Pulse magnetization of soft magnetic alloys with rectangular hysteresis loops," 
in which the dependence of magnetization time on pulse current amplitude was considered, and in which ex- 
perimental data were given on the switching (magnetization) time of cores prepared from strips of various thick- 
nesses, 


In his communication, “Experimental investigation of commutators implemented by magnetic elements 
with rectangular hysteresis loops,” V, N, Tutevich provided experimental data from the investigation of commutator 
used in remote control systems, and characterized the influence exerted by the form of the supply pulses on the 
effectiveness of commutator functioning, 


The communication of V, V, Verigin, N, S,Kartseva, and S, P, Maslov, "Matrix memory devices,” was 
devoted to the principles and peculiarities of performance of ternary memory devices, The authors dealt with 
the functioning of blocks of devices of this type, 


A. S, Fedorov gave a communication on, “Increasing the speed of action of ferrite core memory devices," 
in which he cited several methods of increasing both the speed of action and the reliability of matrix type memory 
devices, The author based his choices on the optimal relationships between speed and capacity of the memory 
devices, and he provided experimental results of memory device testing, 
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In the paper of M, P, Sycheva, “Certain questions pertaining to increasing the reliability of magnetically 
operating memory devices,” there was investigated the ferrite memory device used in the BESM computer, and 
comparisons were made between various magnetic materials used for memory device cores, Also given were 
results of experimental investigations of models of memory devices, 


V. A, Zhozhikashvili and K, G, Mitiushkin gave a paper “On a method of designing pulse distributors 
(commutators) made of magnetic elements with rectangular hysteresis loops,” The design is done by the method 
of obtaining the minimum power necessary for transmitting information from one element to another, taking into 
account the load power, The design method allows one to find the optimal distributor parameters for the most 
economic expenditure of energy in the distributor, 


In the paper “Relay phenomena in ring circuits containing magnetic cores with rectangular hysteresis 
loops,” the same authors considered the questions of stability and transient response in ring circuits, 


V. V. Kobelev and Iu, I, Vizun, in the paper “Program-controlled magnetic memory devices,” reported 
on a shift register circuit of ferrite-transistor cells, using plane crystal triodes, which made it possible to obtain 
currents sufficient for the reading and writing of information in the memory devices, 


lu. S, Volkov presented a paper, “Inductive transformer of angular shaft rotations into digits, expressed in 
the binary system in the form of pulse sequences," The transformer gives shaft rotations with an accuracy of 
three minutes of arc in the form of thirteen-place binary numbers, 


V. A, Ziminym, in the paper "Shaping transformers,” showed that the use of transformers for pulse shaping 
provides for stable machine functioning, and also for simplicity of the circuits, The author gave a circuit for a 
conversion transformer and a mathematical analysis of its functioning. Formulas were derived for designing the 
shaping transformer; experimental data were given, 


Ia, G, Koblents,in the paper “Construction of devices for the automatic control of an ATS (Automatic 
Telephone Station) based on contactless magnetic elements with rectangular hysteresis loops,” considered the 
requirements on commutation elements, the construction of elementary operational circuits, registers, random 
access circuits, subscriber assemblies, circuits for recognizing the number of the calling subscriber and transmitting 
information along communication lines, 


V. P. Sabadashev, in the paper “ Contactless tele-control devices with phase selecting indications,” gave 
an analysis of phase-sensitive circuits with saturated transformers, obtained their loading characteristics analytically, 
and gave experimental data and recommendations on the application of various remote control circuits with phase 
selecting indications. 


In the paper “On the construction of circuit elements not requiring materials with rectangular hysteresis 
loops," V. A, Leokene dealt with the use of semiconducting diodes with large nonconducting portions (insensitivity 
threshold) for the construction of economical and operative computing cells from materials without rectangular 
hysteresis loops, 


V. P, Sal'nikova, in the paper "Technological problems of the preparation of ferrite cores with rectangular 
hysteresis loops,” reported on experiments and research concerning the development of new types of ferrites and 
their application to the mass production of ferrite cores. 


The paper of Eng, lu. la, Dusavitskii, on the subject, “Applications of contactless magnetic relays for the 
automatic control of furnaces for metal heat-treating,” was devoted to the use of contactless magnetic relays 
as measuring organs in automatic current control circuits. The paper included information on the method of 
selecting the block schematic of the device to be regulated, and provided design and experimental data on the 
magnetic amplifiers used, 


In his paper, "On magnetic memory devices,” F, V, Maiorov considered several circuits for controlling 
magnetic memory devices by means of semiconducting elements and ferrite cores: circuits for reading and writing 
amplifiers using transistors, circuits for transistor diode matrices taking coded numbers from a simple diode de- 
coder for the coded address, writing and reading circuits of crystal triodes, and reading register circuits, also con- 
taining semiconducting triodes, The peculiarities of functioning, and basic data, were considered for the circuits 
cited, 
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In the paper of M, I, Sataev, "Transformation of shaft rotation to a number and the inverse transformation 
of a number to a shaft rotation,” the author gave the principles of transforming shaft rotations to binary numbers 
in a digital machine by means of miniature transformers constructed of semiconductors and ferrites, and also de- 
scribed miniature inverse transformers and gave the basic characteristics of these devices. 


L, 1, Gutenmakher presented a paper, "Work of the Electro- Modelling Laboratory (Laboratorii elektro- 
modelirovaniia) in the domain of creating contactless elements and systems," The basic subject of the 
Laboratory's work is the creation of information machines with large, fast-acting memories with the purpose 
of simulating complex processes of mental work, In the paper there was considered the design of an information 
machine and the specific properties of its elements, in particular, new types of economic address systems, The 
use of operational ferrite memory systems allows the machine dimensions to be decreased, There was developed 
in the Laboratory a universal digital machine without mechanical moving parts, the memory device of which 
was implemented by ferrites controlled by magnetic elements without recourse to vacuum tubes. 


In the paper of I, B, Negnevitskii and R. A. Lipman, "Transient responses in magnetic amplifiers in the 
relay regimen,” there was given the physical basis and a quantitative analysis of the delay phenomena in 
operating and releasing a contactless relay, It was shown that the lag in load current variation occurs until 
the core induction attains the saturation induction value; thereafter, the ordinary transient response occurs, in 
accordance with the relay type. Expressions were obtained for computing the total times of operation and re- 
lease of magnetic relays, 


N, M, Tishchenko, in the paper "Contactless relay times," dealt with the development of one type of 
contactless relay, and gave its design and experimental data. 


The work of the Electro- Modelling Laboratory in the application of ferrites with rectangular hysteresis 
loops was detailed in the paper of A, A, Kosarevy. The author noted that a disadvantage of existing ferrite-diode 
designs is the necessity of using diodes which were not specially designed to work with very small ferrites, 
as a result of which multi-turn windings are needed, thus reducing the speed of action of the design. Diode-less 
variations of the device were considered, in particular, transfluxors, However, as the author noted, ferrite- diode 
designs are today more efficient, and are recommended to a higher degree, than are ferrite-transistors. 


N,V, Korol’kov gave a short communication on a two-core trigger developed at the Electro- Modelling 
Laboratory. 


E, Dzhakov, Corresponding Member of the Bulgarian Academy of Science, considered, in his paper, the 
times of relays with magnetic amplifiers having delay times of 10 seconds and higher, The delay time of the 
relay is attributable to an RC network at the magnetic amplifier input, An analysis was given in the paper of 
the influence of the various circuit parameters on the delay time, and relationships were given for choosing the 
designs of the magnetic amplifiers used in the given relay. 


At the conclusion of the plenary session, E. T, Chernyshev, N, G, Chernysheva, and E. N, Chechurina gave 
a paper on the subject “Contemporary state of questions on the testing of magnetic materials in dynamic regi- 
mens,” The authors considered methods of obtaining magnetic characteristics with an estimate of the possible 
errors of measurement, 


A paper on the subject “Fundamental principles of constructing series toroidal cores for magnetic ampli- 
fiers and contactless elements of automation,” was presented by O, A.Sedykh and M, A, Rozenblat, Presented 
in the paper were the fundamental principles of constructing series toroidal cores for magnetic amplifiers and 
contactless magnetic elements for automation, and there was given a method for selecting the various core di- 
mensions, The authors gave a design for series toroidal cores recommended for line production, meeting the 
basic requirements of industry for toroidal cores for magnetic amplifiers and contactless magnetic elements 
by means of the least possible number of standard- dimension mass- produced cores, 


V. E, Pankin gave a communication on the results of development of series standard toroidal cores for 


magnetic amplifiers. 


The representative of the Czechoslovakian Academy of Science, lu, Gashkovets, in his paper, gave a de- 
tailed discussion of the development of magnetic amplifiers in Czechoslovakia, 








The representative of the Polish Academy of Science, S. Vendzhin, spoke of the use, in Poland, of 
magnetic amplifiers in industrial automation devices, 


In conclusion, the participants at the Conference adopted a resolution in which it was noted that magnetic 
elements are very promising elements of automation, remote control and computing technology, and a series of 
measures was adopted directed toward improving research work in the domain of magnetic elements and their 


most intensive injection into industry. M. A. Boiarchenkov, V. S, Volodin, F. 1. Kerbnikov 
. . . 7 a ° . ’ 
G. D. Kozlov, G, V. Subbotina, and I, &, Trefilova. 





